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Abstract 

This article is devoted to a study of the asymptotic dynamics of generic 
solutions of the Einstein vacuum equations toward a generic spacelike 
singularity. Starting from fundamental assumptions about the nature 
of generic spacelike singularities we derive in a step-by-step manner the 
cosmological billiard conjecture: we show that the generic asymptotic 
dynamics of solutions is represented by (randomized) sequences of het- 
eroclinic orbits on the 'billiard attractor'. Our analysis rests on two 
pillars: (i) a dynamical systems formulation based on the conformal 
Hubble-normalized orthonormal frame approach expressed in an Iwa- 
sawa frame; (ii) stochastic methods and the interplay between genericity 
and stochasticity. Our work generalizes and improves the level of rigor of 
previous work by Belinskii, Khalatnikov, and Lifshitz; furthermore, we 
establish that our approach and the Hamiltonian approach to 'cosmo- 
logical billiards', as elaborated by Damour, Hennaux, and Nicolai, can 
be viewed as yielding 'dual' representations of the asymptotic dynamics. 
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1 Introduction 

Remarkable developments took place in the sixties, seventies, and early eighties, that 
molded our understanding of singularities in General Relativity (GR). On the one hand, 
the singularity theorems of Hawking and Penrose [I] proved the inevitability of spacetime 
singularities under rather general conditions. On the other hand, in a series of papers, 
Lifshitz, Khalatnikov and Belinskii, [21 |3l S] and references therein, set out to give a de- 
scription of the actual nature of generic singularities — henceforth we will refer to these 
authors and their work as BKL. These authors performed a heuristic analysis that eventu- 
ally resulted in the claim that a generic spacelike singularity for the Einstein field equations 
with a perfect fluid with a radiation equation of state as the matter source is spacelike, 
local, vacuum dominated, and oscillatory. BKL obtained this picture by (i) using certain 
spatially homogeneous (SH) metrics, (ii) replacing constants with spatial functions, (iii) in- 
serting the resulting expressions into Einstein's field equations and making a perturbative 
expansion, (iv) checking if this ad hoc procedure yielded a consistent result and thereby 
completing the basic viability test. Furthermore, BKL employed synchronous coordinates, 
i.e., Gaussian normal coordinates, such that the singularity occurred simultaneously. The 
BKL approach led to the conjecture that the time evolution of a generic solution in the 
vicinity of a generic singularity is schematically described by a sequence of generalized 
Kasner solutions (i.e., vacuum Bianchi type I solutions where constants are replaced with 
spatially dependent functions), where the transitions are 'mediated' by generalized vac- 
uum Bianchi type II solutions through the so-called Kasner map. In subsequent work, 
this map was shown to be associated with chaotic behavior [UEKT]. 

The BKL picture obtained further heuristic support from the Hamiltonian approach devel- 
oped by Misner and Chitre [8l[9l[10], originally for asymptotic Bianchi type IX dynamics. 
In one variety of this approach the dynamics was described in terms of a free motion in an 
abstract flat Lorentzian (minisuper-) space surrounded by potential walls, or alternatively 
by means of a spatial projection leading to a free motion inside a potential well described 
by moving walls; for further developments of this picture, see [11] and [121 Chapter 10]. 
In another variety of the Hamiltonian approach an intrinsic time variable was introduced, 
which led to a description of the asymptotic dynamics in terms of a projected 'billiard' 
motion in a region of hyperbolic space bounded by infinitely high straight stationary walls, 
see e.g. |9]. This work on 'cosmological billiards' was later generalized in order to deal with 
general inhomogeneous cases, which culminated in the recent work by Damour, Hennaux, 
and Nicolai, see [131 114j and references therein. Apart from these studies, special inhomo- 
geneous spacetimes with non-oscillatory singularities have also been investigated by means 
of Hamiltonian methods — notably by Moncrief, Isenberg, Berger, and collaborators, who 
also obtained numerical support for the general basic BKL picture, see [15], and [16] for a 
review and additional references. 

Despite the ingenuity of the BKL and the Hamiltonian methods, there has been, unfortu- 
nately, little progress as regards a desired sharpening of the heuristic arguments in order 
to turn conjectures into rigorous mathematical statements The requirement of mathe- 
matical rigor is met in a different approach to cosmological singularities: the dynamical 
systems approach. During the last decade, based on dynamical systems formulations, 

^Note, however, the successful rigorous treatment of special non-oscillatory cases in |16l 117). which is 
partly based on Fuchsian methods. 
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there has been considerable progress in the SH case as regards theorems about dynamical 
behavior, largely in connection with the book "Dynamical Systems in Cosmology" |12j : 
notably, Ringstrom obtained the first mathematical theorems in the context of oscillatory 
behavior for Bianchi type VIII and, more substantially, type IX models |18t I19j. In an at- 
tempt to extend the dynamical systems approach to the inhomogeneous context, Uggla et 
al. [20] introduced a dynamical systems formulation for the Einstein field equations with- 
out any symmetries — in the following, we will refer to this work as UEWE. The results 
were: a detailed description of the generic attractor; concisely formulated conjectures 
about the asymptotic dynamic behavior toward a generic spacelike singularity; a basis 
for a numerical investigation of generic singularities — following UEWE numerical results 
yielded additional support for the expected generic picture as well as the discovery of new 
phenomena and subsequent refinements [211 [22l |23] . 

The purpose of the present paper is two-fold: first, we establish a link between the Hamil- 
tonian picture as described by Damour et al. |13^ [13] and the dynamical systems approach 
to inhomogeneous cosmologies as initiated in UEWE. In particular, we demonstrate that 
the 'Hamiltonian billiards' and the corresponding dynamical systems description can be 
viewed as yielding dual representations of the generic asymptotic dynamics ('configura- 
tion space description' vs. 'momentum space description'). To avoid excessive clutter 
that would obscure the main ideas we simplify our presentation by confining ourselves to 
the four-dimensional vacuum case. Inclusion of matter sources such as perfect fluids, see 
e.g. [20l [26], and to more general cases like those considered by Damour et al. [131 (Hj 
should — in principle — be straightforward. 

The second and main purpose of this paper is to derive and give rigor to some of the 
key conjectures formulated in the BKL and Hamiltonian approach, starting from 'first 
principles' connected with the full state space picture of our dynamical systems approach. 
This derivation does not constitute a mathematical proof of the conjectures; however, 
in many respects we go beyond what has been accomplished previously: we identify the 
'billiard attractor' as the attractor for the asymptotic dynamics toward a generic spacelike 
singularity; we give decay rates that describe the approach to the attractor; we show how 
asymptotic constants of the motion, which were obtained from the Hamiltonian billiard 
approach initially, arise. An important ingredient in our treatment is the use of stochastic 
methods: we emphasize the connection between genericity and stochasticity. The detailed 
statements and conjectures we obtain are accessible to numerical experiments and can be 
compared with results for special models. 

The present paper is to a large extent self-contained. First, we give a condensed review 
of the Hamiltonian approach of Damour et al. [13^ [H] in Section [2l Second, we present 
the conformal Hubble-normalized orthonormal frame approach leading to the dynamical 
systems formulation in Sections [3] and [H Endowed with the basic techniques and the 
interrelation between the two approaches, we then derive the billiard attractor, which 
then serves as the starting point for further discussions, developments, and concluding 
remarks. 
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2 The Hamiltonian billiard approach 

In |13j . Damour, Henneaux, and Nicolai used a Hamiltonian approach to study the dy- 
namics of the Einstein-dilaton-p-form system in the neighborhood of a generic spacehke 
singularity. The authors described the asymptotic behavior of the fields by a 'billiard' 
motion in a region of hyperbolic space bounded by straight 'walls'. The techniques used 
in [13] represent generalizations of methods that are due to Chitre and Misner [HI [9l [10] , 
which have been applied and extended by many authors, see the references in [131. In 
the following we give a brief review of the Hamiltonian 'billiard approach' restricted to 
the four-dimensional vacuum case in GR; the presentation is based on [13], but the nota- 
tion is tailored to our later purposes. In particular, for frame indices we use a,(3,. . ., as 
in [m [20], and for spatial coordinate indices we use as in [T3] and [T2[ [20]. 

In [13], the metric is written in 3-|-l form with vanishing shift vector, 



ds"^ = -iV2 + dx' dx^ . (1) 

Let g^^ denote the inverse of the spatial three- metric gij and g = det gij . We introduce 

gijkl g^( k _ ^^^^ ^^^^^ ^^^^ ^^^^ _ . 

gijki ^YiQ DeWitt metric [27J multiplied with g^^^"^ so that it becomes a tensor instead 
of a tensor density; Qij^^G"'^''^ = 6ij^^ = <5|f(5^.). 

The Lagrangian (density) associated with the reduced Einstein-Hilbert action reads 

c = IVgN-^ G''^' 9im + ^N^R, (3) 

where a dot refers to the partial derivative w.r.t. x^ and '^R denotes the spatial three- 
curvature associated with g^j. To obtain the Hamiltonian TC we introduce the conjugate 
momenta 

TT^^ =dC/dgi, = '^N-'g'^'''gki, 

where we define 

N = g-^/'^ N . (4) 

This leads to 

C = ir'^gij -n = Tr'^ij - N [g^.M tt'^tt''' -g^R . 

Variation of TC w.r.t. yields the Hamiltonian constraint TC = 0, while the momentum 
constraints are introduced separately, see [T^ . 

There exists a unique oriented orthonormal spatial coframe {lj" | a = 1 . . . 3}, a;" = dx*, 
such that = E/3^"/3-^^m where V is a diagonal matrix, which we choose to express as 
V = diag[exp(— 6-^),exp(— 6^),exp(— 6^)] , and where AA"j a unit upper triangular matrix, 

/I M\ 

AA-i) = I 1 AA^g 1 = I 1 ns I . (5) 
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This choice of frame leads to the so-called Iwasawa decomposition of the spatial metric 
9ij- 

a 

Existence and uniqueness of the frame {co"} is associated with the theorem on the unique- 
ness of the QR decomposition in linear algebra; A^"j can also be viewed as representing the 
Gram-Schmidt orthogonalization of the spatial coordinate coframe {dx^}. The Iwasawa 
decomposition corresponds to a Cholesky decomposition of a symmetric matrix A into a 
product BF R, where the diagonal elements of the triangular matrix R are factored out 
and parametrized as exp(— 6°). 

To avoid confusion with the Greek frame indices a, /?, . . ., we prefer to use h as the kernel 
letter for the diagonal degrees of freedom instead of (5, which was used by Damour et 
al. [13]. Note that the negative sign in the exponentials is in agreement with the conven- 
tions of [13], but contrary to the conventions of e.g. |1H I12j. For representations that are 
adapted to metric anisotropics, see Misner [8] and e.g. [3 [12]. Following the summation 
convention of [13j, summation of pairs of coordinate indices i, j, ... is understood, whereas 
sums over the frame indices a, /3, ... are written out explicitly in this section. 

The frame {cq,} that is dual to {cj"} is given by Cq = ead^i = exp{b'^)Af^adj:i, where the 
matrix {N^a) is the inverse of (M'^i), 

1 AA^2 AA^aX /l ni /l -ni nin^ - 

M\]=\Q 1 AA^g = 1 na = 1 -n^ j . (6) 

1 / \0 1 / \0 1 

Expressing the Lagrangian ([3]) in the Iwasawa frame variables 6" and N'^i yields 



-1 



Y.Qaph'^b^ + \ exp(2(6^ - 6")) (aA^ATV 

o,/3 o</3 



where we have introduced the reduced metric Qafi-, which is associated with the diagonal 
degrees of freedom, i.e., '6"-space': 



when '6"-space' is endowed with the metric Qaf3, it is isometric to a (2 + l)-dimensional 
Minkowski space, since the signature of Qap is ( — h +). The inverse metric Q°^^ of Qaj3 is 
given by 

Y Q'^^VaWp = Y ^"'^VaWp " | E ^""^P = E '"^'^'^ ~ I ( E ( E '^'^ 
a,(3 a, (3 a, (3 7 a (3 

The momenta tTq that are conjugate to 6" and the momenta V^^ {i > a) that are conjugate 
to AA"i read 

/3 f3 



-8- 



2 The Hamiltonian billiard approach 



while = when i < a; accordingly, P*^, can be viewed as a lower triangular matrix, 
whose components we for convenience call Vf. 



Inverting the relationships ([7]) yields 



V\ 



0> 
Pi 



(8) 



6" = ^-NY^g^P^^ , n, = Ne'ib'-''^ (n,Vi + {e'('' + nl) V2) , (9a) 

ni = iV e2(*'-*') {Vi + na P2) , ^3 = iV e2(^'-^') , (9b) 
which leads to 



-AT i5:a"Vvr^ + i5:e-2(''^-^'^)(nAA'^.)'-ffi? ,(10) 
where the Hamiltonian TC can be read off easily. 

We assume that the spacetime with metric ([T]) possesses a spacelike singularity in the past; 
then, as argued in [ISj, 6" is expected to be timelike in the vicinity of this singularity, 
i.e., /3 < 0. Based on these considerations, it is possible to replace the metric 

variables 6" by new variables: we introduce = — ^ Q^pb'^b^ and 'orthogonal' angular 
variables, collectively denoted by 7, i.e., is orthogonal to dp . In these variables the line 
element associated with 0^/3 can be written in the form 

da^ = g^pdh'^db'^ = -dp^ + p'^dnl , 

a,/3 

where dQf^ is the standard metric on hyperbolic space. Making a further variable change 
according to 

A = logp = ilog(- (11) 

a,l3 

yields da'^ = p'^{—d\^ + d^\) and thus in particular 



a,(3 



[-7^1 + TT, 



cf. ((TOj). 

Consequently, choosing the rescaled lapse according to N = p^ leads to the Hamiltonian 

W = i[-vr! + vr^]+p^5:^c^e-W7). (12) 

here, the sum is over a number of terms of the same type: ca are functions of spatial 
derivatives of the metric, off-diagonal metric variables, and momenta; WA{"f) denote linear 
forms of the variables 7°, i.e., ^^(7) = (^A)a 7^- 
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The Hamiltonian approach relies on the expectation that p — > +00 in the approach to the 
singularity. As a consequence of this assumption, each term eyiY>[—2pwA{l)\ becomes a 
sharp wall, i.e., an infinitely high potential which is described by an infinite step function 
0oo(a^) that vanishes for 2; < and is infinite when x > 0. In the GR vacuum case 
there exist three 'dominant' terms in the sum, which is the minimum number of terms 
required to define the 'billiard table'. The coefficients ca of the dominant terms are non- 
negative functions of the variables, but 'generically' they are positive, see jl3] . The three 
dominant terms /3^ca exp[— 2/)t(;A(7)] thus generate 'dominant' walls 6oo(— 2i(;a(7))- Only 
the dominant walls are assumed to be of importance for the generic asymptotic dynamics; 
the other, 'subdominant', terms of (jl2p . whose exponential 6"-dependence can be obtained 
by multiplying dominant terms, are dropped accordingly. We thus obtain an asymptotic 
Hamiltonian of the form 

3 

Woo = \ [-vrl + + ^ eU-^WAil)) , (13) 

A=l 

where the sum is over the three dominant terms. 

This limiting Hamiltonian is believed to describe generic asymptotic dynamics. It is 
independent of AA"j, P*^, and A, which suggests that V^^^ ^^"i are asymptotic 
constants of the motion. The remaining non-trivial dynamics resides in hyperbolic space, 
i.e., in the variables 7. The asymptotic dynamics can be described via as a geodesic 
motion in hyperbolic space, which is constrained by the existence of three sharp refiective 
walls, i.e., by a 'billiard motion'. The asymptotic dynamics is thus given by a 'cosmological 



bilhard,' see Figure 8(a) 



3 The dynamical systems approach 
Conformal Hubble-normalization 

Physics is a science of scales: in many problems there exists a variable scale of particular 
importance, which can be factored out of the equations in order to obtain a simpler 
mathematical description. In GR there exists a single dimensional unit which is chosen to 
be 'length' (or, equivalently, 'time'). The geometrical way of factoring out a scale in GR 
is by means of a conformal transformation: choose a conformal factor that depends on the 
variable scale so that the conformal factor carries the dimension and the conformal metric 
becomes dimensionless. Then, for dimensional reasons, the Einstein field equations split 
into decoupled equations for the conformal factor and a coupled system of dimensionless 
equations for quantities associated with the dimensionless conformal metricJl This reduced 
dimensionless system carries the essential information about the problem, since one can 
solve for the decoupled equations for the conformal factor once it has been understood, 
and thus recover the physical metric. 

^Logarithmic derivatives of the conformal factor occur in the equations; however, one of the evolution 
equations (which is the Raychaudhuri equation in the present Hubble-normalized case) can be used to 
algebraically solve for the logarithmic time derivative, and the constraints can be used to generically solve 
for the logarithmic spatial derivatives. 
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The problem of generic spacelike singularities and the dynamics of solutions in a neigh- 
borhood thereof is one example where one has a variable scale: it is provided by the affine 
parameter of inextendible causal geodesies, or the expansion along such geodesies. Here 
we consider non-rotating timelike congruences that may or may not be geodesic; the lat- 
ter is of no particular relevance, since the generic behavior toward a spacelike singularity 
is similar, as will be discussed later. One equation plays an essential role in the singu- 
larity theorems — the Raychaudhuri equation for the expansion. This makes a suitable 
function of the expansion an excellent candidate for a conformal factor. By factoring out 
the expansion, which blows up toward the singularity, we will obtain a system of regular 
dimensionless equations!! Let us begin by introducing the basic set-up and nomenclature: 

Consider a four-dimensional Lorentzian manifold (the 'physical spacetime') with a metric 
g of signature (— + -I-+) that satisfies the Einstein vacuum equations. By convention we set 
c = 1 and SttG = 1, which leaves 'length' (or, equivalently, 'time') as the single remaining 
unit. Let us consider a timelike non-rotating congruence in the physical spacetime whose 
tangential vector field we denote by u. We choose an orthonormal frame (a = ... 3) 
such that eo = u, and introduce a orthonormal coframe dual to which we denote by cj", 
i.e., cj'^(ef,) = 5°'ij. Our choice of coordinates is that of a coordinate 3 + 1 decomposition 
associated with u, where we set the shift vector to zero and let d^o = Nu, so that the flow 
lines of u are the timelines. Accordingly, 

eo = N~^d^o , Go = ejd^^ and = Ndx^ , = e°dx* 

with a = 1,2,3 and i = 1,2, 3. The metric thus reads 

g = riabu^"^' = -N\dxy + 5^pe%e^.dx'dx^ , 

where r]ab = diag(— 1, 1, 1, 1). The Hubble scalar of the vector field u = eg is given by 
H = IVqU"^ where V denotes the covariant derivative associated with g, and hence H is 
related to the expansion 6* of u by = \e. Note that R-^ carries dimension 'length' 
(see \14\ 125] and compare with discussions about FRW models where this is a frequently 
used fact). 

The natural way of keeping track of dimensions in a parametrized theory like GR is by 
means of using rigid measuring rods that carry the dimension, which corresponds to as- 
signing a dimension ('length') to the orthonormal one-forms. To obtain Hubble-normalized 
variables and associated Hubble-normalized equations, we factor out the Hubble scale out 
of the frame vectors and dual one-forms, i.e., we go over to a conformal Hubble-normalized 
orthonormal frame da and its dual coframe fi": 

rl 

The dimensionless 'unphysical metric' G that is associated with the conformal frame is 
defined as G = r]ah it is related to g through a conformal rescaling: 

G = r^ab n^f^' = H^Vab^^u^' = H\ . (14) 

The conformal Hubble-normalized orthonormal frame satisfies 

ao = M~^d^o , da = Ejd^. and fi" = Afdx^ , fi" = E^x' , (15) 

■^For historical reasons we will choose the Hubble variable instead of the expansion, but since they are 
proportional to each other this makes no essential difference. 
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where N is the conformal lapse, are the conformal spatial frame vector components, 
and i?" the associated dual vector components, i.e., 

M = HN, ^"' = ^. E^i = He'',. (16) 

Since fi" is dual to da, i.e., il°'{dp) = 5"^, the same is true for the components and 
i?", i.e., E^-^Ejj"^ = d^i^. The unphysical metric G can be written 

G = Tiab ^''^^ = -Af^dx^ + SapE'^iE^jdx'dx^ . 
For further details and discussions see |24j . 

The starting point for the conformal Hubble-normalized orthonormal frame approach to 
the Einstein equations is the decomposition of the commutators of da, 

[do,da] = ifado + Fjdp, Fj = qda'' -^a" -ea^yR"' , (17a) 

[da,dp] = {2A[a5f3]^ + eapsN'^)d^ , (17b) 

which can also be written in the alternative form 

= {da + Ua)do - (<5„^ do - Fj) dp , Fa'' = q 6 J - llj - ej^ W , (HZi) 
= dp , CoP = e,^^ {d^ - A^) - Nj , (dlbl) 

see [24J. In ()17ap . are the frame components of the acceleration of d^ w.r.t. the 
conformal metric G, i.e., U = V^p^o, where V is the covariant derivative associated with 
G. The variable q is the negative Hubble scalar of the vector field d^ in the unphysical 
spacetime, i.e., q = —\0, where © = '^ad^ is the expansion of do w.r.t. the metric G. 
The object is the (trace-less) shear, is the Fermi rotation, which describes how 
the frame rotates w.r.t. a Fermi propagated frame, associated with the vector do and 
the conformal metric G. In (|17bp . A^"^ and are spatial commutator functions that 
describe the three-curvature of G, see below; we also refer to [121 [2l] where the analogous 
non-normalized objects are described. 

Using the analog of (|17ap in the physical spacetime, i.e., [eo,eQ,] = iia^Q + fj^e^, where 
the trace of fj^ equals cf. [21], it is straightforward to derive the equation 

doH = -{l + q)H . (18a) 

We thus see that, in terms of the physical spacetime, q is interpreted as the deceleration 
parameter associated with eg = u and g. As we will describe below, q can be algebraically 
determined by the Raychaudhuri equation. Analogously, we define the spatial Hubble 
gradient by 

daH = -raH, (18b) 

where is generically determined by the Codazzi constraint, see below. For dimen- 
sional reasons, the equations ()18p for H decouple from the remaining dimensionless equa- 
tions ()19p - ()2ip . which are given next. 



We are now in a position to describe the dynamics of Einstein's vacuum field equations in 
terms of the conformal Hubble-normalized variables introduced above. From the Hubble- 
normalized commutator equations, the Einstein field equations, and the Jacobi identities. 
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see |24] . we obtain a dimensionless system of coupled equations; we split the system into 
gauge equations, evolution equations, and constraint equations: 



Gauge equation: 



Ua=dalnM. (19) 



Evolution equations: 

doEj = EjEf,\ (20a) 

do^ap = -(2 - q)^ap - 2e-y\^ S^)^ Rg - ^S^p + (/s)a/3 , (20b) 

^oA^"'^ = {3q - 2F^(°) iV^)^ + (/^)°^ , (20c) 

doA^ = EJ + {lA)a ■ (20d) 

Constraint equations: 

= Caf^ Ep' , (21a) 

= 1 - S2 - Ofc - (/g) , (21b) 

= -SA^S^'^ + ea^^J^fs^Ns^ + {Ic)a , (21c) 

= A^iV/ + (/j)„, (21d) 

= {Ij^)a . (21e) 

Here, (. . .) and (. . .) denote symmetrization and trace free symmetrization, respectively; 
let us briefly comment on this system: the equations (fT9|) . pOap . and piap are obtained 
from the commutator equations; ()20bp is the trace free Einstein evolution equations, while 
the remaining Einstein evolution equation, which can be written as the Raychaudhuri 
equation, yields the expression for q that is given below in (j22ap : (|21bp is the Gauss con- 
straint; (|21cp are the Codazzi constraint equations; (120cp . (|20dp . (I21dp . (121eP are obtained 
from the Jacobi identities. 

In these equations we have employed the following definitions: 

S2 := , g = 2S2 + (/,) , (22a) 

nk:=-l^n, = -6A^ - ^B^'' + {Ir) , (22b) 

S^f, = 2A^ Ns^^e^)"'^ + + {Is)ap , B^p = 2N^^ N^p - N^^ N^p ; (22c) 



3 



here and in the following the norm of a spatial vector is written as Va V°' = V'^ ; 
^TZ, ^Sap are the conformal scalar curvature and trace free Ricci curvature, respectively. 
Finally, the expressions for (/*).., which are zero in the symmetry adapted spatially ho- 
mogeneous (SH) case, are given by 

ils)ap = -e^\a Np^s - 2r7) + {d^^ + A^^){Up) + 2r^)) + f7(„f/^) - 2r(^r^) , (23a) 

(/g) = -i(a, + f/„-2^, + 2r„)(f7- + r"), (23b) 

[InT^ = (a^ + f77)e^'^("F/) , {lA)c. = \{dp + Up){?,q5/ -Ea^), (23c) 

{Ig) = i (25„ - 4^, + r,)r'^ , (/c)a = + (2,5,^ + S/)r^ , (23d) 

= -1(9(3 Na^ + eJ"'df3A^), {I Ju)a = C^" U p . (23e) 

[IsU = d{aAp)-d^Ns{^efs)^' , {lR)=Ado,A'^. (23f) 
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The above dimensionless coupled system of partial differential equations is associated with 
a state space described by the state vector 

X = {Ej,J:^p,A^,N^p) ; (24a) 

we make the split 

X = {Ej) e S , where S = (S«^, A^^N^p) . (24b) 

The Fermi rotation variables Ra, which describe how the chosen frame rotates w.r.t. a 
Fermi propagated frame, are gauge variables: thus Ra = corresponds to choosing a 
Fermi frame while Ra 7^ yields a rotating frame. The variable Ua is regarded as a 
gauge variable determined by M (or, equivalently, by N); the spatial Hubble gradient 
is generically determined by the Codazzi constraint (j21cp . however, for certain choices of 
M one can derive an evolution equation for and include Tq, into the state vector S, see 
UEWE; sometimes it is also useful to elevate q to a dependent variable, see [21] , 

Iwasawa frame variables 

In this paper we choose the conformal orthonormal frame da = Ea^d^i to be a conformally 
rescaled Iwasawa frame, i.e., Ea^ = H~^ea^, where e^* is a lower triangular matrijiEl with 
the diagonal entries written in an exponential form, see Section [2l 

In Appendix[X]we show that choosing a conformal Iwasawa frame in the conformal Hubble- 
normalized dynamical systems approach corresponds to choosing the gauge 

S23 = —Ri , S31 = R2 , S12 = —R3 , 

and setting 

iVs = iVaa = . 

This makes it natural to introduce the notation = Ti^a, and to replace the off-diagonal 
components of T,af3 by Ra according to the above equation, i.e., 

Sii II12 SisX / Si —R3 R2 \ 

S21 S22 S23 = -R3 S2 -Ri . (25) 
S31 S32 S33/ \ R2 —Ri S3 / 

By this choice of gauge we write the state vector S as 

S=ii:a,Ra,Aa,Nap) . (26) 

For further details we refer to Appendix [Aj where we give the relationships between 
the Hamiltonian Iwasawa variables and the Hubble-normalized variables of the conformal 
Hubble-normalized dynamical systems approach. 

As we will see in this paper, the Iwasawa gauge has non-trivial consequences for the de- 
scription of the generic asymptotic dynamics of solutions in a neighborhood of a generic 
spacelike singularity. However, next we describe some central concepts that are indepen- 
dent of the chosen frame. 

^Since e^^ = exp(fe")A/'*a, the components Ea' form a lower triangular matrix when JV'a is an upper 
triangular matrix, which is in accordance with the definitions of Section [5] For the case of the coframe, 
both E"i and Af^i are upper triangular matrices. 
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4 Asymptotic silence, locality, and the silent boundary 

Asymptotic silence and asymptotic locality 

A generic spacelike singularity is expected to be a scalar curvature singularit}{l associated 
with ultra strong gravity which increasingly focuses light in all directions as the singularity 
is approached. The resulting asymptotic collapse of the light cone^ causes the particle 
horizons to shrink to zero size toward the singularity along any timeline. As a consequence, 
communication between different timelines is prohibited in the asymptotic limit, and we 
therefore refer to this causal feature of shrinking particle horizons as asymptotic silence; 
the associated singularity is said to be asymptotically silent. 

Evidence for the conjecture that generic singularities are asymptotically silent is discussed 
in UEWE and [26]; we also refer to corroborative results from the context of spatially 
homogeneous cosmologies, see \12\ [T8[ \T9\ [25] and references therein. In Section 4 in 
UEWE it is shown that, generically, asymptotic silence is connected with the property 
that the Hubble-normalized spatial frame variables vanish asymptotically, i.e., 

toward the singularity; however, this is not sufficient for asymptotic silence, as is illus- 
trated by special examples in [231 Section 4]. We also refer to [23] for examples of var- 
ious non-generic asymptotically silent singularities and for examples of singularities for 
which asymptotic silence does not hold — a phenomenon referred to as asymptotic silence- 
breaking. 

We define the dynamics along a timeline toward a singularity to be asymptotically local if 

Ej^O, daiS,r0,Uf^)^O, (r„,[/„)^0. (27) 

toward the singularity. (Note that this definition slightly differs from the one used in [22| ) . 
Asymptotic silence and asymptotic local dynamics are closely related concepts. The 
shrinking of particle horizons and the associated loss of communication between time- 
lines suggests that inhomogeneities are irrelevant for the asymptotic dynamics, since they 
are shifted outside the shrinking horizons faster than they grow, i.e., E^^ goes to zero at 
a rate that is faster than the possible growth rate of diS, dilog{H), f?jlog(AA). This is 
supported for generic timelines toward generic asymptotically silent spacelike singularities 
by numerical experiments |21[ 122] . but it may not be the case for all timelines. This is 
an issue that will be discussed in the concluding remarks in connection with 'recurring 
spike formation' [22]; see also the discussion in [23]. However, apart from in the conclud- 
ing remarks, in the remainder of this paper we will be concerned with generic timelines 
approaching generic asymptotically silent singularities for which the dynamics is assumed 
to be asymptotically local. 

The silent boundary 

Consider the Einstein field equations in the conformal Hubble-normalized approach, i.e., 
the system of equations p^ - ([23]) introduced in Section [3| It follows from (|20ap . in cou- 
sin contrast, a non-scalar curvature singularity is expected to require fine-tuning of initial data. 
®We may hence refer to this as an 'anti-Newtonian' limit. 
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nection with ()21ap . that 

EJ = (28) 

defines an invariant subset of the state space; we refer to this invariant subset as the silent 
boundary, which is characterized by the state vector S. Since da = E^di, Equation (f28l) 
imphes daS = and daUp = 0, daTp = in p^ - ([23|) . so that the equations on the 
silent boundary reduce to a system of ordinary differential equations. Note, however, that 
we obtain a system of ODEs for each spatial point x*; hence the silent boundary can be 
visualized as an infinite set of copies, one for each spatial point, of a finite dimensional 
state space. The equations on the silent boundary are identical to the equations for 
the state vector S of spatially self-similar and SH models. This is a direct consequence 
from the symmetry assumptions, since spatial self-similarity implies that daS = Q and 
9a{Ua-,ra) = in a symmetry adapted frame, which in turn entails that the equations for 
Ea decouple, so that one obtains a reduced coupled system of ODEs for the variables S. 

The SH subset of the silent boundary (in brief: SH silent boundary) is defined as the 
invariant subset of the silent boundary given by 

EJ = Q and = , = , . (i28l) 

Since this yields (/*).. = 0, the equations on the SH silent boundary reduce to a system 
of ordinary differential equations, which is obtained by setting the quantities (/*).. to zero 
in ()19p -(|23p. This system is identical to the (reduced) system of equations for SH models, 
since spatial homogeneity implies daS = Q and f/a = = in a symmetry adapted 
frame, which subsequently leads to a decoupling of the quantities Ea\ see also [29]. Note 
again that the SH silent boundary is an infinite set of copies, parametrized by the spatial 
coordinates, of a finite dimensional state space with S as state vector. 

The particular importance of the SH silent boundary in the study of generic spacelike 
singularities stems from its connection with asymptotically local dynamics; compare (j27p 
and (j28ip . The reasoning is as follows: The Hubble-normalized conformal orthonormal 
frame approach produces field equations that are regular in the asymptotic limit toward 
a generic spacelike singularity. This allows us to extend the state space to also include 
the silent boundary in our analysis of the dynamical system. This is highly advantageous; 
in particular, based on our expectation that generic spacelike singularities are asymp- 
totically local, see (j27p . we are now able to conjecture that the asymptotic behavior of 
asymptotically local solutions is reflected in the dynamics on the SH silent boundary; see 
the discussion in Section [101 Therefore, the SH silent boundary and its neighborhood in 
the full state space will be the main object of investigation in this paper. 

In contrast to the spatially self-similar and SH models, the variables S on the SH silent 
boundary depend on the spatial coordinates, and the constants of integration are therefore 
spatially dependent functions. The dynamical systems formalism thus explains the first 
step in the heuristic ad hoc procedure by BKL, which consists in replacing the constants in 
certain SH models with spatially dependent functions. Moreover, the ODE structure on the 
silent boundary in turn induces ODE structures associated with perturbations thereof — in 
a series expansion spatial partial derivatives only act on spatial functions associated with 
previous lower order terms, thus only yielding temporal constants for the perturbative 
equations; this feature naturally captures the BKL statement that spatial derivatives only 
enters 'passively' into the equations (BKL [H pp.656-657]). Note, however, that the 
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dynamical systems approach allows one to derive the results of BKL and thus put them in 
the rigorous context of the regularized dimensionless state space picture of Einstein's field 
equations, a picture which also simplifies comparisons with analytical results for special 
cases as well as with numerical investigations. 

The asymptotic ODE structure induced by asymptotic local dynamics can be exploited in 
several ways; in particular, it allows us to reparametrize the individual timelines, for which 
the spatial coordinates are fixed, so that dof = —df/dr along a given timeline, where 
f{x^,x^) is any variable occurring in an ODE, and t{x^,x^) := — log{£/i) is a 'local time 
function' directed toward the singularity. In this context, i = g^/^ and g is the determinant 
of the physical spatial metric. In our convention, hatted objects refer to objects that are 
functions of the spatial coordinates alone; note also that we choose a time direction toward 
the singularity, which is in contrast to UEWE. To obtain the solution in the chosen time 
coordinate x^ , one subsequently integrates the relation dx^ = —J\f~^dT so that 

x"^ = 5P{x')- r M~^(T\x')dT' . (29) 

J To 

The reparametrization freedom allows us to introduce a dynamical systems treatment of 
the equations on the silent boundary and its neighborhood, in complete analogy with the 
treatment of the so-called silent cosmological models from a dynamical systems perspec- 
tive, see [301131] and [12^ Chapter 13] III From now on we will use the local reparametrization 
freedom induced by asymptotically local dynamics and write 

do = -dr , 

to study the asymptotic dynamics along timelines by means of finite dimensional dynamical 
systems techniques!! 



5 Kasner circle stability and the oscillatory subset 

It is suggestive that there exists an open set of vacuum models that possess asymptotically 
silent singularities and obey asymptotically local dynamics for a generic set of timelines. 
This follows from the analysis of UEWE |20j in combination with the results of [26] , where 
a dynamical systems formulation was used that added the null geodesic equations to the 
Einstein field equations; further support comes from the numerical analysis of [21] and 
from [22l[23]. We therefore conclude that at each spatial point the SH silent boundary — 
or rather a subset of the SH silent boundary — constitutes a local attractor for solutions in 
the full inhomogeneous state space associated with the Einstein equations; note, however, 
that it can only be a local attractor in the full state space since there exists, e.g., an open 
set of solutions without any singularity at all |^34j; moreover, there may exist an open set 
of solutions with weak null singularities, see [35] and references therein, or an open set 
of solutions with a singularity that consists of parts that are asymptotically silent and 

'^In the present context it is irrelevant that the only nontrivial spatially inhomogeneous non-rotating 
exact silent solutions without a cosmological constant turned out to be the Szekeres dust models; see [32] 
and, for similar topics and additional references, [33) . 

®As a simple illustration, we derive the generalized Kasner line element of BKL in Appendix [B] 
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spacelike and parts that are weak null singular; for special examples of singularities of this 
latter type, see [23] . 

The main aim of the present paper is to identify the subset of the SH silent boundary 
that acts as the attr actor for generic asymptotic dynamics. As suggested by the analysis 
of SH models, the subsets that represent Bianchi type I and II models will be of crucial 
importance. The common denominator of these models is the most fundamental structure 
of the SH silent boundary state space: the Kasner circle. 



The Kasner circle, stable variables, and the oscillatory subset 

The SH part of the silent boundary = is given by = Tq, = ; it consists of an 
infinite number of identical copies of the finite dimensional state space spanned by the 
state vector S = (Sq,^, ^q,, A^^^t/?); hereby, the spatial variables can be regarded as acting 
as an index set. The associated system of evolution equations and constraints is obtained 
when the quantities (/*).. are set to zero in ()20bp - ()20dp and ()21bp - ()21dp . respectively. 
This system of equations possesses a one-parameter set of equilibrium points which is 
fundamental for our analysis: the Kasner circle K°. It is determined by 

l-^^ = Aa = N^p = Ra = 0. (30) 



Since Ra = corresponds to Sq,^ = {a ^ (3), Sq,^ is diagonal for every point on K^, 
i.e., T,af3 = diagpi, ^is], where the numbers (q = 1 . . . 3) are constants. (All 
considerations apply to each of the (infinitely many) copies of the finite dimensional SH 
state space. Since the numbers Sq, (a = 1 ... 3) may differ between the individual copies of 
the SH state space, are in fact spatially dependent functions. In this context recall that 
we employ hatted objects to denote temporally constant spatially dependent functions.) 
It is standard to represent the Kasner circle in terms of the (generalized, when viewed as 
spatially dependent) Kasner exponents pa, 



diag[i;i,S2,S3] = diag[3pi - l,3p2 - l,3p3 - 1] , (31) 



where we have omitted the hats on top of pa in order to agree with standard notation. 
Since tr = and = 1, the Kasner exponents satisfy the Kasner relations 

Pi + P2 + P3 = 1 , Pi + P2 + P3 = 1 • 
In Appendix |B] we derive the generalized Kasner line element, see [1], 

from the fixed point solution given by (I30|) and ()31l) . 

The Kasner circle K^ can be divided into six sectors, where each sector is characterized 
by a typical ordered sequence of the Kasner exponents (pa)] hence, sector (123) is defined 
as the part of K° where pi < P2 < P3, whereas sector (312) is characterized by the 



order ps < pi < P2, etc.; see Figure 1(a) Evidently, the six sectors can be identified 
with each other through permutations of the spatial axes. Of particular interest are the 
six special points on the Kasner circle that are associated with solutions possessing an 
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additional symmetry, the so-called locally rotationally symmetric (LRS) solutions. The 
points Qa correspond to the three equivalent LRS solutions whose intrinsic geometry is 
non-flat; for Qi we have (Ei,S2,S3) = (-2,1,1) and {pi,P2,P3) = cyclic 
permutations yield the parameters for Q2 and Q^. Each point Tq, corresponds to a flat 
LRS solution: the Taub representation of Minkowski spacetime; the Taub points are 
given by (Si,S2,^3) = (2,-1,-1) or (pi,P2,P3) = (1,0,0), and cyclic permutations; see 



Figure 1(a) 




(231) 




(312) 



(123) 



Ri,R2,R3 Ni,Ri 

(321)\ / (132) 

Ri,R2 




(a) Kasner sectors and unstable variables. 



U 




Qi T2 Qi Ti Q2 T3 Qi 

(b) Variation of the Kasner parameter u along the Kasner 
circle. 

Figure 1: The Kasner circle of fixed points; sectors and unstable variables and relation to 
the Kasner parameter u. Note that n = 1 at the Qa points and that u = 00 at the Taub 
points Tq. 

It is convenient to parametrize the Kasner exponents in terms of the Kasner parameter u, 
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which can be defined frame invariantly through 

det(Ea/3) = 2 - ^j__^^_^^2)3 ' ^ l^' ' (^^) 

see Appendix [C] for details. Due to frame invariance, the Kasner parameter naturahy 



captures the equivalence of the six sectors of K^, see Figure 1(b) On sector (a, /3,7), 
where Pa < Pp < P'y, we set 



-a 1 + u u{l + u) 

^° = T"^ \ 2 ' Pl3 = 1~, \ 2 ' = \ 2 ^ ^ °° • (^^) 



It is easy to check that Pa + Pp + Pj = 1 and Pa+p'^+p'^ = 1- The parameter u = 1 
describes the points Qi,Q2iQ3j i-e., the three equivalent representations of the non-flat 
LRS Kasner solution, while u = oo defines the Taub points Ti, T2, T3 and thus the Taub 
solution Ifl 

We proceed by performing a local dynamical systems analysis in a neighborhood of the 
Kasner circle K^. Linearization of the dynamical system at an arbitrary point {pi,P2,P3) 
of K^ yields 

drEa = —3(1 — Pa)Ea (no sum ovcr a) (34a) 

for the conformal frame variables, and 

drAa = —3(1 —pa)-Aa (no sum over a) (34b) 
drNap = -6(1 - p^)N^p (a / /3 / 7 ^ a) , (34c) 

for the variables A^, Nap (a 7^ (3) on the SH silent boundary; furthermore, 

drNi = -QpiNi , drN2 = -6p2 N2 (34d) 

drRi = 3{p2 - P3)Ri , drR2 = -3(^3 " Pi)-R2 , ^^^3 = ^{pi " P2)R3 , (34e) 

where we set A''i = A''ii and = ^22 (recall that N3 = N^s = 0). We see that the 
variables Ea^, N^p [a 7^ 0) and Aa belong to the stable subspace of each fixed point of 
K^ (except for the Taub points). In contrast, the variables {Ri, R2, R3) and (iVi, A^2) are 
stable or unstable depending on the sector of K^ the point {pi,P2,P3) lies in. Finally, the 
variables = Sq-q, belong to the center subspace, i.e., they are constant to first order. 



The analysis of the stability of the Kasner circle K^ is summarized in Figure 1(a) , where 
the unstable variables are given for each sector of K^. 

We now decompose the state vector S = (Sa/3, Aa, Nap) of the SH silent boundary into a 
'stable' and an 'oscillatory' part: S = ^stable © -S'osc, where ^stable = {Nap,Aa) (a / (5) 
and Sosc = C^a, Ra, ^1, ^2)- The subset determined by S^stabie = on the SH silent 
boundary is an invariant subset, which we call the oscillatory subset O. 

We conjecture that there exists an open set of solutions whose behavior is governed by the 
Kasner states as r — > 00, i.e., we expect solutions to spend an increasing amount of time 



^For the BKL definition of u, the order of Pa is fixed according to pi < P2 < Pa; we liave found 
it convenient to use tlie above frame independent definition instead, and to permute tfie ordering of Pa 
according to the sector one considers when dealing with frame dependent matters. 
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close to K° as r — > cxd. Equations ()34bp and (j34cp suggest that the variables -Eq,* and the 
stable variables Sstabic = ^a) (a / /3) decay rapidly, and that {Ea, ^stable) ^ in 

the limit r oo; this motivates the notation 'stable'. In anticipation of results to come we 
note that the behavior of the variables Sqsc = C^a, Ra, ^i, ^^2) on the oscillatory subset, 
by contrast, can be best described as 'oscillatory', whence the chosen nomenclature. Since 
we conjecture (-E^*, -Sstabie) 0, it is the complicated dynamics of Sqsc on the oscillatory 
subset O that will play the decisive role in our description of the generic asymptotic 
dynamics — we thus begin by closely investigating the dynamical system on O. 



The oscillatory dynamical system 

Recall from Section [5] that the silent boundary contains an infinite number of identical 
copies of the SH state space and thus of the oscillatory state space O; hereby, the spatial 
variables act as the index set. Although it is essential for our ultimate aims to consider 
this infinite collection of spaces O, this is irrelevant for our proximate purposes: in the 
present context we may regard the oscillatory subset O as one finite dimensional state 
space. 

The evolution equations on the oscillatory subset O are given by setting the stable variables 
to zero in Eqs. (|20bP - (p0cj) . i.e., ^stable = (^^«/3,^«) = (a / (Recall that the 
quantities (/*).. are zero because O is a subset of the SH silent boundary determined by 
Ea^ = 0, [/q, = Tq, = 0.) Accordingly, the dynamical system on O takes the form 

(35a) 
(35b) 
(35c) 
(35d) 
(35e) 
(35f) 
(35g) 
(35h) 

where 

= 1 (S? + + + 2Rl + 2Rl + 2RI) 
^Sn = U^i- N2){2Ni + N2) , (35i) 

'522 = -U^l - N2){m + 2N2) , '533 = -U^l - N2f . 

The Gauss constraint (|21bp becomes 

l_s2-^(iVi-iV2)2 = o, (35j) 
while the Codazzi constraints (j21cp result in 

iVii?2 = 0, Arii?3 = o, iV2i?i = 0, N2Ri = ^. (35k) 
Note that — a priori — the Codazzi constraints read 

A^ii?2 = 0, iV2i?i = 0, (iVi - iV2)i?3 = ; (IMH) 





= 2(1 - s2)Si - 2(i?i 


+ Rl) + 3511 


dr^2 


= 2(1 - S2)S2 -2(/Z? 


— i?|) + ^^22 


dr^3 


= 2(l-S2)S3 + 2(i?? 


+ Rl) + 3533 


drRl 


= 2(1 - S2)i?i + (S2 - 


- ^3)Ri + 2R2R3 


drR2 


= 2(1 - S2)i?2 - (S3 - 


-^l)R2 


drRs 


= 2(1 - S2)ii3 + (Si - 


- T,2)R3 — 2R1R2 


drNi 


= -2(S2 + Si)iVi 




drN2 


= -2(S2 + S2)iV2 , 
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however, unless A'^i = N2 7^ on a finite r-interval, the constraints ()35kp ensue. Now 
assume that A^i = 7^ 0, so that 1 — = 0; then i?i = and R2 = 0, and the 
equations (j35ap - ()35hp imply Si = S2 (and 9t-Si = dr^2)i whereby = 0. Therefore, 
N1R3 = and N2R3 = 0, and (l35kl) holds in the case TVi = iV2 7^ as weh. 

The constraints determine the structure of the oscillatory state space O: it consists of a 
number of invariant subsets ('components') that are connected with each other by parts of 
their boundaries only. This partitioning of O is to be understood in terms of the Bianchi 
classification: the invariant subset given by (A^i = 0) A (A'^2 = 0) is the Bianchi type I state 
space, the components (A^i = 0) A (A'^2 7^ 0) and (A''i 7^ 0) A (A''2 = 0) describe Bianchi 
type II states, and the components N1N2 < and N1N2 > are the Bianchi type VIq and 
VIIq subsets, respectively. In other words, the equation system (|35]) . when restricted to 
the respective component, is identical to the equations for the Bianchi type I, II, VIq, VIIq 
vacuum models. In the following we more closely examine these 'Bianchi components' of 
O. 

The Bianchi type I component is defined by (A'^i = 0) and (A'^2 = 0), and hence {Ri, i?2; 
are arbitrary, but subject to the Gauss constraint. Every orbit on the Bianchi type I 
component describes a Kasner solution, since 

^72 = = and 1 - T,"^ = = N^^p iV"^ = (36) 

on this set; we therefore denote the Bianchi type I component alternatively as the silent 
Kasner subset K,. However, since {Ri, R2, R3) 7^ (0,0,0) (except for the fixed points on 
K^), every orbit on /C is a representation of a Kasner solution in a frame that rotates 
w.r.t. a Fermi propagated frame; this is a consequence of the Iwasawa frame choice; see 
AppendixOfor details. We distinguish several invariant Bianchi type I subsets which are 
defined by the vanishing of different combinations of the three variables (i?i, i?2) -^3)- We 
will use ;S as a kernel letter for an invariant subspace and a subscript that denotes the 
non-zero variables associated with the different subsets. Accordingly, Br^ is an invariant 
subset, where Ra 7^ 0, Rjs = 0, i?^ = (a 7^ /3 7^ 7 7^ a); for Br^r-^ we have Ri 7^ 0, 
i22 = 0, i?3 7^ 0; the set Br^r^r^ is characterized by R1R2 7^ or R2RZ 7^ and 
thus comprises the case of all three variables being non-zero. (It is easy to see that 
the equations ([35]) prevent the set (i?i = 0) and the set (i?3 = 0) from being invariant 
subsets.) The Kasner circle K^ is a special subset of the Bianchi type I component; here, 
(i?i, i?25 R-i) = (0, 0, 0). In the standard notation B^, for the closure of a set, we find that, 
e.g., Br, = Br, U KO, or, /C = Br.r^r^ = Br,r,r, U Br,r, U Br, U Br, U Br, U KO. 

The Bianchi type II component is defined by either A^i or A'^2 being zero. The Codazzi 
constraints (|35kp enforce N^Rp = (07^/?), i.e., the Bianchi type II component consists 
of four invariant subsets: Bn„ Bjy,, B]y,R,, and Bjy^R,. Solutions on the former subsets 
represent Bianchi type II models in a Fermi propagated frame, since Ra = for all a; 
solutions on Bn^r, , BN2R2 representations in a frame that rotates w.r.t. a Fermi frame. 

Finally, we denote the Bianchi type VIq and VIIq subsets by Bn,N2- ^"^^ ^NiN2+i i'^- 
spectively, where the subscript denotes the sign of N1N2. The Codazzi constraints (j35kp 
enforce Ra = for all a; hence all solutions on these components are represented in a 
Fermi propagated frame. 

In Figure [2] we present a diagram containing the subsets ^B* introduced above. We give the 
dimension of the different subsets and show how they are related to each other by setting 
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Figure 2: Subset contraction diagram of the oscillatory subset O. In our notation, the sub- 
scripts denote the non-zero variables, and D hence describes the dimension of the subsets 
B^:. Solid lines correspond to setting i?Q,-variables to zero (frame rotation contractions); 
dashed lines are associated with setting A^i or to zero (Bianchi type contractions). 
Note that, e.g., B]\j-^ includes the two possible representations (one for each possible sign 
of Ni) of the Fermi propagated subset; similar statements hold for the other subsets that 
involve A'^i or N2. Analogously, e.g., Br^ is the union of two subsets that are characterized 
by the sign of Ra- 



variables to zero; the figure thus represents a contraction diagram. 

We conclude this section by noting that models of Bianchi type VIII and IX are not 
described by the equations on the SH silent boundary we consider. This is due to the 
inherent incompatibility of the SH frames of the Bianchi type VIII and IX models with the 
Iwasawa frames we use throughout this paper; w.r.t. an Iwasawa frame Bianchi type VIII 
and IX solutions appear as inhomogeneous solutions and are thus associated with general 
solutions in the full interior state space X. In contrast, the SH frames of Bianchi type 
VI_i/9 models are Iwasawa compatible as we shall see in Section [I3l Since these models 
are expected to possess the same attractor as generic timelines of generic asymptotically 
silent models when expressed in an Iwasawa frame, it is rather these SH models that are 
of interest in the present context. 



6 Dynamics on the components of the oscillatory subset 

In the following we will describe the dynamics of the system ()35p on the oscillatory subset 
O. Since O decomposes into independent components, we will analyze each component of 
O separately. We will see that the solutions of (I35p are heteroclinic orbits, which, except 
for in the Bianchi type VIIq case, connect different fixed points on the Kasner circle and 
thus provide transitions between different Kasner states; accordingly we will henceforth 
refer to such heteroclinic orbits as transitions. 
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The silent Kasner subset (Bianchi type I subset) 



The silent Kasner subset, /C, is defined as the invariant component of O given by (A^i = 0) 
and (A''2 = 0) and thus by the conditions ([36]) . Consequently, solutions on fC represent 
Kasner states, though in general not in a Fermi frame, since (Ri, R2, R3) 7^ (0,0,0). In 
this way, a given Kasner state has several different representations on /C. 



We call the solutions of (j35p on /C frame transitions. Depending on how many of the 
variables {Ri, R2, R3) are non-zero, we distinguish between single frame transitions and 
multiple frame transitions. 

There exist three families of single frame transitions, denoted by Tr^, Trj, Tr,^, which 
are associated with the subsets Br-^ , j -Sfla • In the case of a Tr^ transition the equa- 
tions (j35ap - (j35cp yield that = const and that S^, are monotonically increasing or 
decreasing, see Figure [3l The initial and final states of the transitions are equilibrium 
points on the Kasner circle K^ that are related by (Sa)+ = (£„)_, (S/3)+ = (S^)-, 
and (S^)+ = (S/3)-, where the subscripts + and _ denote the final and initial states of 
the transition, i.e., (S^)-^ = limT-^-too ^^^(r). Consequently, a Tr^ orbit corresponds to 
a Kasner solution viewed in a frame that rotates around the a-axes; in their final state, 
the P- and 7-axes are rotated by tt/2 w.r.t. their initial position, which corresponds to an 
interchange of the /3-/7-directions. Since the Kasner parameter u is a frame invariant, it 
is necessarily invariant during a frame transition; in particular. 



U = Ua 



const ; 



(37) 



see Appendix O for details. Expressed in u we obtain for Tr^, 



|AS;3| = |(S;3)+-(S^). 



' 3(l + 2n)(l + n + 'u2)-i 
3(2u + n2)(l + n + u2)-i 
3(-l + n2)(l + u + n2)-i 



(38) 



depending on whether the transition starts from sector (/?, 7,0), (/3, 0,7), or (q,/3,7). A 
simple consequence is |AS/}| < 2\/3. 

For multiple frame transitions more than one of the variables {Ri, R2, R3) is non-zero. 
There exists one class of double frame transitions, Tr^r^ , which is associated with Br^r^ , 
and one class of triple frame transitions, Tr^r^r^, connected with Br^r^r^] for the latter, 
Ra ^ Va. The initial state of a multiple frame transition is an equilibrium point on sector 
(321) of K^, while the final state lies on (123); for details see Appendix |Dj The Kasner 
parameter u is invariant under multiple frame transitions, since it is a frame invariant; 
furthermore, we have | AS2I = and (Si)-|- = (S3)=p (as for Tr^ transitions). For examples 



of multiple frame transitions, see Figure 4(a 



The silent Bianchi type II subset 

Solutions of the oscillatory system ()35p that satisfy A^^ ^ for some a are located on 
the Bianchi type II, VIq, or VIIq component. These solutions are characterized by non- 
vanishing three-curvature, i.e., ^IZ = —NapN'^^ + ^{N^)"^, see ([22]) . We denote solutions 
of Bianchi type II and VIq as curvature transitions; for solutions of Bianchi type VIIq we 
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(d) Tjv2 (e) Tr^ 



Figure 3: Projections of single transitions for the Iwctsa-Wci frame onto dicigonal XjQ;-spa,ce. 
The first three transitions, T/Vu Tr^-, '^Rs-, are the transitions relevant for the billiard 
attractor in Section [TTl the two other ones, T/Vj , , are of relevance for the approach to 
the attractor, but asymptotically they are suppressed as shown in Section [TOl 

do not use this terminology for reasons that will become clear below. Curvature transitions 
(and Bianchi VIIq solutions) differ from frame transitions in one important respect: the 
Kasner parameter u changes under the transition. We begin by considering those curvature 
transitions that are associated with the Bianchi type II subset. 

The silent Bianchi type II subset is the subset of O determined by {Ni 7^ 0) A (A''2 = 0) 
and (A^i = 0) A (A''2 7^ 0); it is the union of the invariant subsets Bni, SniRi, and Bn2, 

^N2R2 ■ 

The simplest solutions on the silent Bianchi type II subset are the orbits on Bni or Bn2i 
denoted by T/Vi and T/V2 5 which we call single curvature transitions; the absence of a third 
family, T/V3 is due to the peculiarities of the Iwasawa frame. In the following we con- 
sider Tjvi transitions; the treatment of T/V2 transitions is analogous. For TjVi , the Gauss 
constraint p5jp yields Nl/12 = 1 — S^. The resulting system for (Si, S2, S3), Eqs. (j35aP " 
(|35cp . possesses solutions whose initial states are equilibrium points on sectors (123) or 
(132) of K^, see Figure [3l Without loss of generality we consider T/Vi transitions from 
sector (123) (as T/Vi orbits that originate from sector (132) are easily obtained by a permu- 
tation) . The orbits can be parametrized by u = tt_ , the Kasner parameter characterizing 
the initial state; in terms of an auxiliary function C(r) we obtain 

El = -4 + (1 + ^2) C, T.2 = 2-u^C. S3 = 2-C, (39a) 
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where the evolution equation for ( is given by 

9.C = 2(1-S2)C with (1-S2) = -^(C+-C)(C-C-) , (39b) 

where C± = 3/(1 =F + u'^)] by definition, < C- < and < C+ < 3. The function ((t) 
interpolates monotonically between C_ (as r — > — oo) and C+ (as r — > oo). Letting u = 1 



in (i39]l yields the orbit that connects Qi with Ti, see Figure 3(a) . u = oo corresponds to 
the point T3. 

Evaluation of (j39ap at r = =poo (corresponding to C^) yields the initial/final states of the 
T/Vi transition, and thus a map that connects the initial and final Kasner states, which 
are described hy u = U- and u+, respectively: 

f n_ - 1 if u_ > 2 

u+ = < . (40) 

\ (u-- 1)-^ if 1 < n_ < 2 

This formula was first obtained by BKL via different methods; we will refer to (j40p as the 
Kasner map. Note that this result holds for both T/v^ and T/V2 curvature transitions. 

For later purposes we define the growth factor g according tc0 



_C+_ l + u + u"^ _ l + u.+ul 



Since u = £ [1, 00) for TJvj orbits, we have 1 < g < 3. 

There exist two classes of mixed frame/ curvature transitions: orbits on BmiRi and orbits 
on Bn2R2^ which we denote by T^-^ji-^ and T^^ji^, respectively. The initial state of a T^^^r^ 
orbit is an equilibrium point on sector (132), and the final state lies on (213) or (231); for 
T/V2_R2 transitions, the numbers 1 and 2 are interchanged; see Appendix [Dl In analogy to 
single curvature transitions, the initial and the final Kasner state are related by the Kasner 
map ()40p . This fact relies on the frame invariance of the parameter u: T/v^ and Tj^^^i^ 
transitions are representations of the same Bianchi type II solutions — they are merely 
expressed in different frames (Fermi frame vs. rotating frame); the transformation law (j40p 
of the frame invariant u is thus unaffected. For more details on mixed frame/curvature 
transitions, see Appendix [Dl for examples, see Figure |4(b)[ Note that henceforth, for 
brevity, we will refer to mixed frame/curvature transitions as mixed curvature transitions. 



The silent Bianchi type VIq and VIIq subsets 

The silent Bianchi type VIq subset is given by N1N2 < and Ra = for all a; it is 
denoted by BniN2-- We call the orbits on BniN2- double curvature transitions: T/ViATj-, 
or simply T/v^jVj- In AppendixO we show that the a-limit of each T/ViAf2 transition is the 
Taub point T3; the w-limit is a fixed point on sector (312) or (321) of K^. Accordingly, 
in terms of the Kasner parameter n, we find 

TjViiVa (Bianchi VIo) 

U- = 00 > n+ G [l,ooj ; (42) 



^°The quantity g should not be confused with the determinant of the spatial metric — it should be clear 
from the context which one is meant. 
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Si 




(a) Double frame transitions Tr^e.-^^. 



(b) Mixed frame/curvature transitions TniRi ■ 



Figure 4: Figure (a) shows a projection of double frame transitions onto S^-space; note 
that single transitions form the boundary of this subset and that the final state of a 
double frame transition can be obtained in terms of a series of single transitions that are 
associated with the same initial Kasner state. Figure (b) shows mixed frame/curvature 
transitions for which analogous remarks hold. 

note that for each u £ [l,oo) there exist double curvature transitions that map the Taub 
state M_ = oo to exactly that final Kasner state, i.e., n = n-(_. It is evident that the behavior 
of double curvature transitions is fundamentally different from the behavior of Bianchi 
type II transitions (single and mixed curvature transitions) ; compare (|42p with the Kasner 
map (j40p . We observe that there do not exist double curvature transitions with arbitrarily 
small amplitudes: there exists e > such that maxr |A^i(t)| > e (and max,- |A'^2(7")| ^ e) 
uniformly for all double curvature transitions; for a proof see Appendix iDl 

Finally, let us now consider the silent Bianchi type VIIq subset B]\f-^M2+} which is defined 
by N1N2 > and Ra = for all a. The state space contains a line of fixed points, £3 , 
given by the conditions (Ei, E2, S3) = (1, 1, —2) and A''i = N2; see AppendixiDjfor details. 
The a-limit of each Bianchi VIIq orbit orbit is a point on C^, the w-limit is a fixed point 
on sector (312) or (321) of K°. Hence, since Bianchi type VIIq orbits do not connect two 
fixed points on K^, and we therefore refrain from calling them transitions. For this reason 
we will not introduce the nomenclature Tn^n2+- 

7 Sequences, eras, and phases 

In the previous section we have described the dynamics on the oscillatory subset by an- 
alyzing the solutions (transitions) on the different components of O. These transitions 
form pieces of a 'heteroclinic orbit puzzle' that governs asymptotically local dynamics; in 
the following we describe and characterize its key features. 
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Sequences of transitions and Kasner sequences 

We define a sequence of transitions as a heteroclinic sequence on the oscillatory subset 
O, i.e., as an infinite concatenation of transition orbits (heteroclinic orbits) on O. (Note 
that, by definition, Bianchi type VIIq orbits are excluded from these sequences; since the 
initial states of type VIIq orbits do not lie on K^, but on £3 , they cannot be concatenated 
with transition orbits.) Accordingly, a sequence of transitions, which we denote by Sq-, is 
described by St = {l^i^n = {Tq,Ti,T2, . . .)0 where each transition % is either a frame 
transition Tr^, Tji^, Tr^, Tr^r^, Tr^R2R3-, or a single or mixed curvature transition TJvi, 
'^2' '^NiRi, '^2^2) ^ double curvature transition T^^n^. It is understood that the 
w-limit of each transition % (which is a fixed point on K°) coincides with the a-limit of 
the transition T^+i. 

The initial state of each individual transition 7^ is a fixed point on K° (which coincides 
with the final state of Tl-i), and is thus associated with a particular value of the Kasner 
parameter u. Therefore, every sequence of transitions Sq- generates a series of Kasner 
parameters in a natural way. Since frame transitions do not induce any change in u, 
see (f37|l . the series of Kasner states is not associated with the sequence Sr directly, but 
rather with the sequence of curvature transitions of Sq-, which may be regarded as its 
supporting 'skeleton'. This makes it convenient to introduce a second running index that 
consecutively numbers curvature transitions: / = 0, 1, 2, ... . The index I can be viewed 
as a function / : N — > N such that l{i) = U, denotes the number of curvature transitions 
among the transitions {7o,7i, . . . ,7i_i}. 

Let Ul denote the initial Kasner state of the curvature transition; then the transition 
maps Ul to Ui+i, 

;th curvature transition 

Ul > n^+i . 




For single and mixed curvature transitions the Kasner map ()40p applies: 

■f (43) 
if 1 < < 2 . 

A double curvature transition, on the other hand, maps a Kasner state characterized by 
u = 00 to a Kasner state with u E (1, 00). 

Two cases occur: (i) The initial state uq is an irrational number, which is the generic case 
since Q is a set of measure zero in M. Then the sequence {ul)l(^^ is given through the 
recursion ([^3]) : in particular, for all /, ui G (l,oo). A sequence (tx/)/gN of this type we 
denote as a Kasner sequence. Note that a Kasner sequence is associated with a sequence 
of transitions Sr that does not contain any double curvature transitions Tjy-^jy^. (ii) The 
initial state uq is a rational number. Then there exists k < 00 such that u^-i = 1, whereby 
Uk = 00. At this point the recursion defined by the Kasner map ()43|) is interrupted. A 
double curvature transition Tn-^n^ must follow, which yields n^^i G [l,oo). Generically 
the Kasner parameter Uk+i is an irrational number, hence the series {ui)i^k is a Kasner 
sequence given by the recursion ([l3j) : by a shift of origin (i.e., by redefining uq := u^+i) 
we obtain again a standard Kasner sequence. If the Kasner parameter Uk+i is a rational 



^^In order to agree with the conventions of 3 we define the first transition of St to carry the index 
number 0. 



-28- 



7 Sequences, eras, and phases 



number, which is the non-generic case, we are back to square one, at the beginning of the 
loop. 

From this discussion we are led to classify possible sequences of transitions. Sequences 
with a finite number of double curvature transitions are innocuous: they are standard 
Kasner sequences by a simple shift of origin and thus do not require special treatment. 
Sequences that contain infinitely many double curvature transitions are qualitatively dif- 
ferent, but they are highly non-generic since they emerge from repeated fine-tuning: for 
such sequences the final state of each individual Ti^^n2 transition must be represented by a 
rational value of the Kasner parameter. This indicates that these sequences are not to be 
treated on an equal footing with generic sequences, i.e., sequences that are free from dou- 
ble curvature transitions. Henceforth, unless otherwise stated, a sequence of transitions 
Sr always denotes a sequence where double curvature transitions are excluded, i.e., it is 
an infinite concatenation of frame transitions and single and mixed curvature transitions. 
It is sequences of this type that are intrinsically tied to Kasner sequences (ti/)iGN which 
are defined through the recursion ()43p . 



Eras, large curvature phases, and small curvature phases 

Here we investigate qualitative properties of Kasner sequences {ui)i^^; in particular, we 
introduce a partition of the sequence into phases of small curvature and phases of large 
curvature; however, we begin by defining the concepts of an epoch and an era. 

Consider a Kasner sequence {ui)i^^ as given by the recursion (j43p . A Kasner epoch is 
simply defined as an individual Kasner state ui of the sequence. The sequence {ui)i^ii 
possesses a natural partition into eras: an era is a set [/in,^out] 3 ^ such that ui is mono- 
tonically decreasing from a maximal value u™ = ui.^ = {ui.^-i — to a minimal value 
1 < -u"*^* = ui^^^ < 2. The length of an era is given by the number of Kasner epochs it 
contains: L = /out — An + 1- Since = ui — 1 for all / of an era [/in, /out], we have 
L = int(n™), where int(x) is a function that gives the integer part of x, see [5j. Note that 
an era can be, and many will be, of length one. 

We now introduce a somewhat more flexible concept that comprises the concept of an era 
as a special case. Let M 9 r/^j > 1; we define a phase of small curvature — a small curvature 
phase — as a set [/in, /out] 9 I such that ui monotonically decreases from a maximal value 
= ui.^ to a minimal value = ui^^^ according to 

= ui,^ = {ui,^.i - >Vu>l i^l^^^ ry„ + 1 > = ui^^, > Vu • (44) 



When r]u is chosen to be equal to 1 the definition of a small curvature phase reduces to the 
definition of an era; however, the usefulness of the concept of a small curvature phase stems 
from the possibility of choosing r]u^ 1, and this is the choice we will typically make. In this 
case a small curvature phase contains only large values of u. Accordingly, the Kretschmann 
scalar Wi associated with the metric G is small, since Wi = 27ti^(l + n)^/(l + u + n-^)^ = 
27u~^(l + + 0(n~^)), cf. Appendix O Since this quantity determines the magnitude 



of the curvature during a Bianchi type II transition, see Section [T3] and Figure 9(c) , the 



terminology 'phase of small curvature' suggests itself. Finally, note that the length of a 
small curvature phase is defined in analogy to the length of an era as the number of Kasner 
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epochs (i.e., Kasner parameters) it comprises: 

L = lout -hn + l = int(n^"^ - r/„ + 1) « n"^ - 7]u ^ u'" , (45) 

where the last approximation holds only in the special case ^ rju- 

The complement of a phase of small curvature is a phase of large curvature. A set [li, l{] 3 I 
is called a large curvature phase if 

ui.^i > Vu , ui <r]u for all / G [h, k] , n^^+i > r]u ■ 

Accordingly, the length L of a large curvature phase is at least int(u;.) > int(?7u — 1), 
since ui. € {rju — l,rju)- In a large curvature phase, the values of the Kasner parameter 
are comparatively small, i.e., small in comparison with the u-values of a small curvature 
phase (j44p ; it follows that curvature is comparatively large in the case of a large curvature 
phase, see Section [T^ which motivates the nomenclature 'large curvature phase'. Note 
that the concept of a large curvature phase does not exist when rj^ is chosen to be equal 
to 1; recall that in this case the concepts era/small curvature phase are identical. 

Small curvature phases and large curvature phases occur alternately. A large curvature 
phase continues and ends the era that began simultaneously with the small curvature phase 
preceding the large curvature phase. When 1 < < 2, the Kasner parameter number 
{I + 1) begins a new era: if u^+i < rju, then that entire era belongs to the continuing large 
curvature phase; if, and only if, ui+i > rju, then a large curvature phase has ended and a 
new small curvature phase begins. Consequently, a large curvature phase can contain an 
arbitrary number of eras; In the following example, where the choice r/„ = 3.5 has been 
made, the large curvature phase contains two and a half eras. 



small curvature phase large curvature phase 



6.29 5.29 4.29 -> 3.29 -> 2.29 1.29 




Let us now briefly discuss how curvature phases are represented in the state space O. 
Consider a sequence of transitions St = {Ti)i^f^ on O and the associated Kasner sequence 
{ui)ieN- Suppose that I E [/imWt] is a small curvature phase (associated with a value 
r]u ^ 1) of the Kasner sequence (■uz)igN- Since the inverse of the Kasner parameter, i.e., 
u^^, is a measure for the (angular) distance (on K^) from T3 (and equivalently from 
Ti, T2), the condition > rju characterizing the small curvature phase is equivalent to 
(n™)~^ < T]^^, i.e., to the statement that the (angular) distance of the transitions from T3 
(Ti, T2) be less than a given distance r]~^. Suppose that the curvature transition number 
lin (which is the first curvature transition of the phase) is a Tat^ transition that takes place 
in the ^-neighborhood of the Taub point T3. In this neighborhood of T3 only three types 
of transitions are possible: TJvi , T/Vj , and Tr^ , see Figure [3l it is thus immediate from the 
figures that the initial 7/Vi transition generates a sequence where T/Vi transitions and 
or transitions occur alternately while the distance from T3 slowly increases (as u 
decreases by n 1-^ n — 1). For some time, therefore, the sequence is 'captured' in the r/~^- 
neighborhood of the point T3; since the curvature is small in this neighborhood, the phase 
is a 'phase of small curvature'. From (j45p we have that the length of a captured phase 
is given by L ~ when ti™ ^» in other words, the number of curvature transitions 
needed to leave the given ^-neighborhood of T3 is inversely proportional to the initial 
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Figure 5: A small curvature phase. Note that this particular small curvature phase starts 
with a 7}j3 transition followed by a Tr^ transition and a series of alternating TJvi and 
Tfjg transitions; the latter oscillatory behavior is of particular relevance for the billiard 
attractor discussed in Section [TTl 

distance (u™)^^ from T3. For an example of a small curvature phase in the state space 
picture, see Figure [5l 

In our description of generic asymptotic dynamics it will turn out to be important to 
characterize Kasner sequences from a stochastic point of view. In order to do so, in the 
subsequent section, we will consider the space of all possible Kasner sequences. 



8 Stochastic analysis of Kasner sequences 

The recursion formula (j43p that defines Kasner sequences typically generates sequences of 
pairwise different values ui. However, there also exist Kasner sequences (iiz)/GN that are 
periodic; for instance, the sequence 

5 + ^/13 6 7-\/l3 5 + \/l3 

Uq = —= > Ui = —= > U2 = —;= > U3 = -Uo = —;= > • • • 

^/13-1 ^/T3-l ^/V^>-\ \/l3-l 

(47) 

is a Kasner sequence with period 3, see [37| or [12l p. 236]. Note that this Kasner sequence 
is associated with a heteroclinic cycle (or rather, a set of entangled heteroclinic cycles) in 
the state space description of sequences. 

Consider an arbitrary Kasner sequence (ui)i^^. As discussed in Section [3 this sequence 
possesses a natural partition into eras j = 1, 2, . . ., where each era j is associated with an 
initial value u = > 1 of the Kasner parameter, which we denote by n™. The Kasner 
map ([13|) thus generates a sequence (M™)jgN obeying the recursion formula 

where frac(-) is a function that gives the fractional part of its argument. It is known 
that the transformation (j48p is associated with exponential instability: the distance be- 



-31- 



8 Stochastic analysis of Kasner sequences 



tween two points (whose initial distance is small) grows exponentially with the number of 
iterations j, see [5]. 

For the following we regard the sequence as a distribution of points (given by 

the values of u^°) in the interval (1,cxd) C M. By adopting this viewpoint it becomes 
natural to ask whether this distribution possesses a probabilistic description, i.e., whether 
it can be modeled by a continuous probability function. In general, this will not be the 
case, since there exist Kasner sequences that generate finite distributions {vi, . . . ,Vn}, 
i.e., li™ E {vi, . . . ,Vn} for all j £ N. In the example (H7|) we even obtain = v = 
(5 + ^/l3)/(^/l3 - 1) for ah j G N. 

However, the picture changes when we consider not the distribution of one single sequence 
(n™)jgN, but the entirety of all distributions generated by the collection of all possible 
Kasner sequences. This collective distribution is indeed modeled by a continuous proba- 
bility function, which arises as the stationary limit from an arbitrary probability density 
via the map i-^ [frac(ii™)]~"'^, cf. (j48p : see [5]. In other words, u^^ can be regarded as a 
random variable with a specific probability density given by 

w(u"') = . , ^ . , . (49) 

^ ^ log 2 u"i(l + ui'i) ^ ^ 

According to the probabilistic description with probability density (j49p . when we choose 
an arbitrary from an arbitrary era of an arbitrary Kasner sequence, the probability of 
li™ taking a value in (^1,^2) is given by 

1 1 

log2 7„^ u{l + u) 

In order to obtain a probability distribution on a bounded interval the inverse of the 
Kasner parameter is used: we denote the inverse of by x, i.e., x = (u™)~^; by 
definition, x G (0, 1). The probability density of the random variable x reads 

see [5]. 

It is straightforward to adapt the probabilistic description to the concept of small/large 
curvature phases. The Kasner sequence {ui)i^f>i consists of infinitely many small curvature 
phases j = 1, 2, . . ., where each small curvature phase is associated with a value u'" > t]u 
of the initial Kasner parameter (which naturally coincides with the initial parameter of 
the era that begins simultaneously with the small curvature phase j, see (|46p ). Let Xj be 
the inverse of it™, i.e., Xj- = {u™)~^ < r]~^; the distribution of the sequence (xj)jgN in the 
interval (0, 77"^) is described by a random variable x with probability density 

1 1 /"'"^ 

log(l + ??„') 1 + x Jq 



evidently, for r/„ = 1 the distribution t()(x) reduces to ([50]) . Note in particular that 
w{0) / 0. 
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Consider a Kasner sequence (u;);^^ whose distribution of the parameters (n™)jgN has a 
probabiUstic description (in terms of the probabiUty density (I5ip ). For such a Kasner 
sequence small curvature phases are expected to dominate over large curvature phases in 
the following sense: the probability that an epoch of the sequence chosen randomly among 
/ epochs belongs to a small curvature phase tends to one as / — > oo. The proof is based 
on stochastic arguments: 

We first compute the expectation value (L) of the length of a small curvature phase. Since 
L u™ — rju = >c~^ — rju, see ()l5|) . we obtain 

(L) = / w{>c)L{K)d>c t()(0) / (x^^"*^ — f]u)d>c = oo , 
Jo Jo 

i.e., (L) is infinite. Note that this feature is independent of the probability density ■w{x) 
as long as w{0) > 0. 

Second, we compute the expectation value for the length of a large curvature phase. A 
large curvature phase continues and ends the era that started simultaneously with the 
preceding small curvature phase, see (I46p . The probability for the large curvature phase 
to contain another (complete) era is given by 

p = l- T" u;(x)dx = w{K)dx = 1 - + ~ 1 _ ^ < 1 ; 
Jo Jvu' log2 log2 

the approximation assumes r/^ ^ 1. Accordingly, the probability that a large curvature 
phase contains at least another n (complete) eras is p"; the probability that a large curva- 
ture phase contains exactly n (complete) eras is (1 —p)p^. (Note here that the assumption 
r/t( ^ 1 ensures that there will be sufficiently many eras in one large curvature phase to 
permit a probabilistic description.) The expectation value for the length of an era that is 
contained in a large curvature phase is 

I / w{>()d>c\ / mt{>c~^)w{>c)d>( > / (x~^ — l) t(;(j<)(ix w (log 2)^"*^ log ry^ , 
\Jvu^ ' J'nZ'^ Jvu^ 

where again » 1 is assumed. Accordingly, a sequence of n eras within a large curvature 
phase is expected to be of length n(log 2)~^ logr/^. Combining the results we find that the 
expectation value for the length of a large curvature phase is given by 

oo 

rju + (log 2)~^ (log 7]u){l -p)^ np^ ^ + Vu log Vu ~ Vu log ??„ < oo . 

n=0 

Since the mean length of a large curvature phase is expected to be finite, while the ex- 
pectation value (L) of the length of a small curvature phase is infinite, we conclude that, 
asymptotically as / — > oo, small curvature phases dominate over large curvature phases. 

Recall that, in the state space description, small curvature phases are phases where the 
sequence of transitions Sr, which generates the Kasner series {ui)i^^ we consider, is cap- 
tured in a small neighborhood of the Taub points, see Figure El Accordingly, dominance of 
small curvature phases over large curvature phases means that the sequence of transitions 
Sj- spends more and more 'time' (as measured by the number of epochs) in a neighbor- 
hood of the Taub points as compared to the time spent close to ordinary Kasner states; 
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above all, the ratio diverges in the limit. Considerations of this type quite naturally bring 
the Taub points into focus when we aim at investigating generic asymptotic dynamics of 
solutions. 



9 Growth 



In this section we analyze the analytic and stochastic properties of a function that is 
naturally associated with a sequence of transitions: the growth function. 



Auxiliary differential equations and the growth function 

In our analysis we will encounter certain types of differential equations; the following 
auxiliary equation represents the common denominator: 

drA = 2{l-T,'^)A. (52) 

When a solution of (I35p on the oscillatory subset O is given, then the equation can be 
regarded as a linear differential equation with a prescribed time-dependent coefficient 
2(1 — S^). Let us now integrate (I52p along the different Bianchi type I and II transitions. 

Along frame transitions we obtain A = const, since (1 — S^) = 0. Single and mixed 
curvature transitions yield A{tq) i-^ A{t) = A{tq)[({t)/({tq)], which follows from ([39b]); 
note that C(''")/C(''"o) > 1 for all r > tq. Letting tq — > — oo and t ^ oo shows that the 
initial value of A (i.e., limT-^_oo ^(t)) and the final value of A (i.e., limT-^oo ^(''")) are 
related by the growth factor g; recall that g is given by 

S = -\ \ 2 ' 

1 — + 

where u is the (initial) Kasner parameter of the curvature transition; see (|4ip . In condensed 
form, we write 

A^ gA (53) 

for the effect of a curvature transition on the quantity A. Recall that, since u G (1, oo) for 
single and mixed curvature transitions, we have 1 < (7 < 3. 

Solving Eq. (I52p along a sequence of transitions Sf on O amounts to iterating (I53p . which 
we will elaborate on in the following. (Recall that, in our convention, St contains only 
Bianchi type I and II transitions, unless otherwise stated.) For the curvature transition 
of a sequence Si the growth factor is given by 

91 ■= r. ■ — 2 = \ — 2~ ^« ' where n := (mmjui -1,1}) ; (54) 

l-ui + uf 1 + Ul+i + uf^^ 

the quantity r; satisfies > 1, where r/ = 1 whenever ui > 2 — this is the case when 
= ui — 1, see (jlO|) . We define the growth function as the product of the growth 
factors: we let G : N ^ M be defined by G(0) = Gq = 1 and 
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Since gj > l,rj > 1 for all j > 1 (and gj 7^ 1 as j — > 00), Gi increases rapidly with / and 
— > 00 with I 00. In our convention, Gi denotes the value of the growth function at 
the beginning of the /^^ curvature transition; see Figure [H 
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Figure 6: Transition counting conventions. 
Iterating (j53p along / curvature transitions of the sequence St thus yields 



z-i 



Ao^G{l)Ao = GiA^=i\{gAA^- 



(55) 



j=0 



hence A is a function whose rate of growth is described by the growth function G. Finally, 
in order to write out Eq. (j55p in proper notation, let Ai denote the value of A at the 
beginning of the i^^ transition % of the sequence Sr = {To,Ti,T2, . . .); then we can write 



Gi(i) Aq ; 



(56) 



recall that = li is the curvature transition index associated with Sq-, i.e., the number 
of curvature transitions among {7o,Ti, . . . ,7^}, see Section [71 The use of this notation 
appears superfluous in the present context, since the function A does not change under 
frame transitions and it would thus suffice to focus on the curvature transitions / = 
0, 1, 2, . . . of St] however, next we will discuss a more general differential equation whose 
solutions are affected by frame transitions. 

Consider the following auxiliary differential equation, which is frame dependent: 



drB = 2(1 



Y?)B + 2Rl 



as with 



t3 - 2Ri ; (57) 

our aim is to integrate this equation along a sequence of transitions Sq-- 

For simplicity, let Sq- = (7o,7i,72, . . .) be a sequence of transitions that does not contain 
T/V2 and TN2R2 transitions; hence N2 = along Sq-, so that 

dri^2 - 2) = 2(1 - S2)(S2 - 2) + 2i?i - 2Rl 

holds for each transition %, which follows from (|35bp and p5jp . Consequently, (S2 — 2) is 
a particular solution of ([5] 



The general solution for B along Sq- can thus be constructed explicitly as the superposition 
of the particular solution (S2 — 2) and the general solution of the homogeneous equation. 
Since the latter is given by ()56l) we obtain 



Bi 



(S2(0-2)+G, 



Bn 



(S2(0) 



(58) 



where S2(i) denotes the value of S2 after i transitions (which, by the convention illustrated 
in Figure m is the value at the beginning of the transition with index i). In the generic 
case, i.e., Bq ^ (^2(0) — 2), we therefore find 



Bi ~ const 



00 



(59) 
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with const 7^ 0; in particular, — > cxd as i — > cxd. The result (j58p is slightly modified 
when T/V2 and TN2R2 transitions are taken into account, see Appendix [El however, the 
generic behavior of B as described by (j59p remains unchanged. 

In Section [7| we have argued that generically sequences of transitions do not contain double 
curvature transitions; in accordance with these considerations we have hence focused our 
analysis on sequences ^7- of this type. In the remainder of this subsection, however, we 
discuss non-generic sequences, in particular in view of the growth associated with double 
curvature transitions. 

From Section [6] and Appendix ID] we recall that the equation dr{2 — E3) = 2(1 — S^)(2 — 
S3) holds along double curvature transitions (i.e., Bianchi type VIq orbits) and analo- 
gously along Bianchi type VIIq orbits. Consequently, we find that A{t) = A{tq){2 — 
E3)(r)/(2 — S3)(to); the behavior of B is identical, since setting {Ri, R2, R3) = re- 
duces (j57p to (j52p . Since the a-limit of all orbits is a fixed point with S3 = 2, the 
expression (2 — S3)(t)/(2 — S3)(ro) diverges as tq —00. The growth factor g associated 
with a Bianchi type VIq /VIIq orbit can thus be considered as being infinite. It is imme- 
diate from these considerations that the growth function is ill-defined along sequences of 
transitions that contain double curvature transitions; it is natural, however, to consider a 
'piecewise' growth function: following the remarks of Section [7] we choose to view double 
curvature transitions as repeated interruptions of a regular Kasner sequence. Thereby 
the concepts introduced in Section [7] (i.e., eras, small/large curvature phases) are still 
relevant; the interruptions, however, correspond to two successive eras not being 'causally 
connected' — the initial Kasner parameter of an era that arises from the final state 
of a TjviATa transition is not related to the Kasner parameter of the preceding era. This 
suggests that we treat sequences of transitions that contain double curvature transition as 
a concatenation of 'pieces' of standard Kasner sequences glued together by Tn^n2 transi- 
tions. Along these fragments the growth function is defined in the standard way; however, 
joining the pieces is not possible, since this is connected with a divergent growth along 
the joining T/ViAf2 transition. 

In anticipation of our considerations in Section [10] we conclude these remarks by consider- 
ing orbits that are not exactly of Bianchi type VIq /VIIq, but arbitrarily close (in a higher 
dimensional state space where the a-/ti;-limits of the VIq/VIIo orbits are saddles). It is 
immediate that the growth of A and B along such orbits can be arbitrarily large depending 
on the degree of the type VIq/VIIo approximation. 

Properties of the growth function 

Since the growth function is the key element in the description of the behavior of A and 
it is appropriate to study this function in more detail. Let us begin by investigating 
the properties of the growth function G{1) along a sequence of transitions Sr with regard 
to the concepts of small and large curvature phases, and eras. 

Consider first a small curvature phase [/in,^out]) where we assume rju > 2. For all I the 
growth factor gi is given by 



91 = 



1 + ui + uf 



l + ui+uf 



(60) 



1 — ti; + n? 



1 + ui+i + ul 



■l+l 
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see The growth function satisfies the recursive relation G^+i = Gigi, so that we 

obtain 



G(/out + 1) ^ TT ^ 1 + n'" + {u'^f 

^ — ''in 



1 + f?« + 



(61) 



where the approximation requires r/„ S> 1. Eq. (j6ip describes the growth of G{1) during 
one small curvature phase. 

In a large curvature phase [k, since ui < t]u for all I £ [k, k], we have 



l-ui+uf l-r]u + r]i 



Cp pa 1 + 2r/-i if ?7« > 1 • (62) 



Therefore, in a large curvature phase, the growth function G{1) behaves like 

I 

G{1 + 1) = G{k) n 5fe > G{k) (Cp)'+i-'' , 

k=k 



(63) 



i.e., we observe exponential growth of G{1) with a given base Cp. 

Let us now discuss the growth of G{1) in an era [/in, /out]- We define the era growth factor 
as 



g(/out + 1) 
G(/in) 



n 91 j 5«out = - 



l + n"^ + (u'°)^ l + u 



)2 1 - tiout + (MOut)2 ' 



(64a) 



so that 



> 



1 + + (u''^)2 



(64b) 



Consider the partition of a sequence into a = 0, 1, 2, . . . eras. For the era with number a 
the era growth factor can be written as 



l + < + (n- 



m^2 



«+l)' 



1 - nr + (nr ) 



(65) 



"a+D 



The era growth function Q is defined as the product of era growth factors. From Qa+i 
QaQa we obtain 



a-l 



/a-1 



Q{a) = Qa = QqW^c = Qq [W{r 



,in\2 



c=0 



\c=0 



l + n'o" + «)2 (ni-)2 



(n-)2 



l + n- + (tx. 



in^2 



where we set = 1 by convention. The expression in brackets is a number in the interval 
(1/3, 3); when and are large the number is close to one. 

Let us now compute the expectation values of the growth factors. The era growth factor g, 
cf. (|64p . does not possess a finite expectation value. This is because the equations contain 



u 



\n\2. 



x^^) and (n"^) 



in^2\ 







w(x)x dx = oo 



(5) = 00 
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where we have used the probabihty density tt;(x) from Section [8j The same result appHes 
for the growth factor describing the growth of G in a smah curvature phase, cf. (|6ip . 
However, the expectation value of the logarithm of the era growth factor can be estimated 
straightforwardly: 



Jo 

By definition, the logarithm of the era growth function Q is the sum of the logarithms of 
the era growth factors, therefore we expect log Q to grow according to 



For further results on growth functions, see Appendix IFl 

10 Asymptotic shadowing 

Shadowing — asymptotic sequences of orbits 

The subsequent considerations constitute the core of our analysis: we provide evidence 
that the asymptotic dynamics of solutions of Einstein's equations which exhibit a generic 
spacelike singularity is represented by sequences of Bianchi type I and type II transitions. 

The basic assumption we make is that a solution of Einstein's equations X{t,x^) that 
possesses a generic spacelike singularity is asymptotically silent. We further assume that 
the dynamics of X{t,x^) along a timeline associated with a generic spatial point of 
such a singularity is asymptotically local. This entails that the asymptotic behavior of 
the solution is described, with an increasing degree of accuracy, by the dynamics on the 
SH part of the silent boundary. The SH silent boundary consists of an infinite number 
of copies of the finite dimensional SH state space that is spanned by the state vector 
S = {'Eap, Aa, Na^), where the spatial coordinates act as the infinite index set. Our 
assumptions thus lead to the conjecture that the asymptotic behavior of the solution 
X{t) at a generic spatial point is governed by the dynamics on the SH state space. 

In Section [5] we have analyzed one structure on the SH silent boundary that is expected 
to be of fundamental importance: the Kasner circle of equilibrium points K°. The local 
dynamical systems analysis of suggests that E^^ — > as r — > oo, which indicates 
intrinsic consistency of our considerations. Furthermore, we are led to the conjecture 
that there exists a 'stable part' instable of the state vector S, whose variables vanish in 
the asymptotic limit, i.e., -^stable = {N12, N23, N31, A^) ^ as r ^ 00. This conjecture 
is intimately connected with the assumption that the asymptotic evolution of a generic 
solution is largely dominated by Kasner states, in the sense that the solution spends an 
increasing amount of r-time in a small neighborhood of the (non-flat part of) the Kasner 
circle. 

Since we expect Sgtabie 0, it is the asymptotic behavior of the remaining 'oscillatory 
variables' Sqsc = (J^a, Ra, ^i, N2) that represents the essential asymptotic dynamics of 




(66) 




a (logg) > 1.77a . 
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the solution X{t). Under the above assumptions the Codazzi constraints (|21cp imply 
that 

NiR2^0, iV2i?2^0, {Ni-N2)R3^0 (67) 

as r — > oo, which is to be compared with (j35kl[) : in the same limit we obtain that Sosc sat- 
isfies the oscillatory dynamical system (j35ap - ()35hp and the constraint pSjP . This suggests 
that the asymptotic dynamics of Sosc = C^a, Ra, -^i; -^2) is governed by the dynamics on 
the oscillatory subset O, which we have shown to be represented by Bianchi type I, II, 
VIq transitions — frame and curvature transitions — and Bianchi type VIIq orbits. In brief: 
asymptotically, the state vector Sosc of a solution X shadows the flow on O. 

Along the orbit S'osc(t") the SH Codazzi constraints ([350) are satisfied only in the limit 
r — > 00; hence the disconnectedness of the Bianchi components of the oscillatory sub- 
set O is irrelevant for the solution S'osc(''")- As an immediate consequence of asymptotic 
shadowing it thus follows that in the asymptotic regime t — > 00, the orbit *S^osc('^) re- 
sembles a sequence of Bianchi transitions and orbits: S'osc(''") can be partitioned into a 
sequence of segments, where each segment is associated with a heteroclinic orbit on O (of 
type I, II, VIq, or VIIq). Two subsequent segments of SosAt) are joined smoothly in a 
neighborhood of an equilibrium point of the Kasner circleo As r — > 00 the segments of 
'S'osc(''") are approximated to an increasing degree of accuracy by heteroclinic orbits on O; 
simultaneously, the joining of segments occurs increasingly closer to K^; hence S^osc(''") 
is described by a sequence of Bianchi transitions and orbits with an increasing degree of 
accuracy. We may call this representation of the orbit S'osc(''") in the asymptotic regime 
an asymptotic sequence of O -orbits ASq (or, for brevity, simply asymptotic sequence). 
Note that the concept of asymptotic sequences ASq gives a precise meaning to BKL's 
'piecewise approximations'. 

As it will turn out in the following, generic asymptotic sequences ASo fall into segments 
that correspond to Bianchi type I and II transitions; it is suggestive to call asymptotic 
sequences of this kind asymptotic sequences of transitions ASq-- With slight abuse of 
notation we write ASj- = (7i)jgNi where % does not denote an exact transition in present 
context, but a segment of ASq, which we can view as an approximate transition. 

However, a priori we cannot exclude that Bianchi type VIq transitions (double curva- 
ture transitions) and Bianchi type VIIq orbits play an important role in the context of 
asymptotic sequences ASq- A segment of ASo that is approximated by a type VIo/VIIq 
orbit — an approximate type VIo/VIIq orbit in our nomenclature — is expected to connect 
a neighborhood of the point T3 with a neighborhood of one of the fixed points on sectors 
(312) or (321) on K^. Clearly, since the accuracy of the approximation, i.e., the degree of 
shadowing, improves as r — > 00, these neighborhoods are expected to shrink with increas- 
ing T. However, this does not automatically lead to an exclusion of type VIq /VIIq orbits 
from generic asymptotic sequences. Instead the exclusion of type VIq/VIIo orbits follows 
from a bootstrap-like argument which we present next. 

^■^The partition of Sosc(t) into segments involves a certain arbitrariness: which point in the neighborhood 
of K° is chosen to be defined as the end point of one transition and the initial point of the subsequent 
one is a matter of convention, however, asymptotically any reasonable convention leads to the natural 
segmentation into the heteroclinic orbits on O. 
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R2 — > and N2 — > as r — > CX3 

Consider an asymptotic sequence ASq- Equations (I35ap ~ (l35h|) imply that the product 
{R1R3) satisfies 

dr{RlR3) = 4(1 - S2)i?ii?3 + (Si - S3)i?li?3 + 2{RI - i??)i?2 

with increasing accuracy as r — > cxd. Using the auxihary quantity B, which is subject to 
the auxihary equation (j57p . we obtain an identical equation for {R2B), i.e., 

dr{R2B) = 4(1 - ^^)R2B + (Si - ^3)R2B + 2(i?i - rI)R2 . 

By choosing appropriate initial data for B, it follows that {R2B) coincides with {R1R3). 
Accordingly we obtain 

R2B = RiR3 or |i?2| ~ l^i^a] B"^ . (68) 

Generically, along an asymptotic sequence of transitions ASj- the quantity B satisfies 
Bi ~ const G^(j) in the asymptotic regime i ^ 00, where Gi is the growth function of 
Section [9l this is because ASj- is approximated by a sequence of transitions Sj- in the 
limit r 00. Along an asymptotic sequence ASq (which, in contrast to ASq-, might 
also contain type VIo/VIIq orbits), the behavior of B described by the growth function 
holds at least piecewise, where the interruptions are associated with the approximate type 
VIo/VIIq orbits of ASq', when the solution shadows a Bianchi type VIq/VIIq orbit, a large 
increase in B ensues, see Section [9l Since type VIq/VIIq orbits thus strengthen the growth 
of B in relation to the normal growth described by G;(j), we conclude that the relation 
Bi ~ const can be regarded as representing a lower bound for the growth of B. 

We can therefore write (168]) as \R2\i ^ const |i?ii?3|. Gj'^j-^ as z ^ 00, or in simplified 
notation, 

\R2\ ~ const \RiR3\ G'^ . (69) 

Since I-R1-R3I is bounded because of the Gauss constraint ( |35jP , Eq. (f69l) implies that 
\R2\ ^ as T — > 00, i.e., generically R2 is asymptotically suppressed. Analogously, 
|A^2| — > as r ^ 00, which we show in the following. 

Consider the product {NiR"^); from ()35a|) - (|35h|) we have 

dr{NiRl) = 2 [2(1 - S^) - - S2] NiRl - AR1R2R3N1 . 

bmce, generically, R2 ~ const {R1R3) G'^ by we obtain 

dr{NiRl) = 2 [2(1 - S^) - - S2] NiRl - ARjG-^{NiRl) . (70) 

In the asymptotic regime r — > 00 the last term in (I70p vanishes; in fact, by using the 
results of the previous sections we show in Appendix [F] that the term falls off rapidly 
enough so that 

dr {NiRl) = 2 [2(1 - S2) - S2 - S2] NiRl (71) 
holds asymptotically. Similarly, using the auxiliary quantity A leads to the equation 

dr{N2A^) = 2 [2(1 - S^) - - S2] iV2^^ 
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for {N2A'^). By choosing appropriate initial data for A, we conclude that 

\N2A^\ ~ const \NiRl\ ; 
equivalently, by using ([56]) . we can write |A''2|i ~ const lA^i-RHi or simply, 

lA^al ~ const \NiRl\ . (72) 



bmce \Rl\ and |iVi -iVal are bounded because of p5jp , Eq. ([72]) implies that, generically, 
|A''2| — > as r — > oo. 

Exclusion of R2 and N2 

Here we establish the following result: the statement (i?2)-^2) ^ (r — > 00) does not 
merely mean that the amplitudes in R2 and are decreasing, but it also leads to a 
complete or stochastic exclusion of the shadowing of orbits that involve R2 or A''2, i.e., 
Trj) '^2! '^2A^2' '^i/?2-R3' '^NiN2 transitions and Bianchi type VIIq orbits. 

In Section E] and Appendix [Dl we have proved that there are no Bianchi type VIq (T/VjATj) 
and VIIq orbits with arbitrarily small amplitudes: there exists e > such that max |A'^2| ^ e 
uniformly for all type VIo/VIIq orbits. Since, generically, |A''2| along asymptotic 
sequences ASq, this entails that Bianchi type VIq and VIIq orbits cannot be contained in 
generic sequences ASq in the asymptotic regime, as for sufficiently large r, |A'^2| satisfies 
|A^2| < £• Since the shadowing of Bianchi type VIq and VIIq orbits is thus completely 
excluded, asymptotic sequences ASq are sequences of transitions ASt, which consist of 
Bianchi type I/II transitions only. In brief: a generic ASq must be an ASq-- 

Let us now consider the remaining transitions that involve the excitation of one of the 
variables R2 or N2, i.e., Trj, Tjy^, Tr^n^, Tr^r^r^ transitions. Consider, for instance, 
the case of Tr,^ transitions: for a given transition the quantity |i?2| goes through 
a finite maximum r along the transition orbit, which can be computed via the Gauss 
constraint P5jp . However, since |i?2| ^ as r ^ 00 by (|69]) . |i22| can never reach the 
value r, and hence the given transition cannot take place when r is sufficiently large. 

It follows that more and more transitions are excluded as r — > 00, so that merely 
transitions in increasingly smaller neighborhoods of T2 and Q2 are possible; however, it 
is a priori not clear whether these transitions 'die out' completely after some time, or if 
there are infinitely many as r ^ 00. Consider, e.g., a neighborhood of the point T2: the 
repeated occurrence of alternating sequences of TJvi and Tr^ transitions (generating small 
curvature phases) is quite plausible a priori, particularly in the light of the statement that 
small curvature phases dominate in sequences of transitions St and thus in the asymptotic 
dynamics of ASq-', see the discussion in Section [8l However, in the following we establish 
a proposition that implies that Tr^ (and TJvj) transitions are excluded stochastically for 
generic sequences ASr'- 

Proposition (Stochastic exclusion of R2 and A'^2)' A generic asymptotic sequence ASq- 
contains a finite number ofT^^, Tr^h^h^, T/Vj, andTj^^^,^ transitions only. 

As for the previous results, the word 'generic' is important also in this context; there might 
exist sequences ASt that contain an infinite number of transitions that involve R2 or ^"2; 
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however, these sequences form a set of measure zero in the space of ah sequences. Recah 
further that asymptotic sequences of transitions ASt are the generic case of asymptotic 
sequences ASq', hence, in the proposition it is equivalent to write ASq instead of ASt- 

The vahdity of the proposition rehes on stochastic arguments, which we present in the 
remainder of this subsection. For simphcity we focus on the statement that there is a 
finite number of Tji^ transitions in a generic asymptotic sequence ASq-- The treatment of 
the other transitions is analogous. 

Consider an arbitrary frame transition. Along the orbit the quantity |i?2| goes 



through a maximum r that is determined by the Gauss constraint ( 35j ) ; it is straightfor- 
ward to express r in terms of the Kasner parameter u that characterizes the transition: 
in particular, if u ^ 1 (which corresponds to a Tjjj transition that takes place in a small 
neighborhood of T2), we obtain r ~ 3m~^; if (u — 1) <^ 1, then the corresponding 
transition is close to Q2 and r ~ (li — 1); cf. Figure [3l 

By Eq. ([69]) . along ASq- the quantity |ii2| decreases with time at a rate given by |-R2|i(j) oc 
G^^; hence, a necessary condition for an (approximate) transition to take place after 
the l^^ (approximate) curvature transition of ASq- is that u be sufficiently large (or (u — 1) 
be sufficiently small): 

An occurrence of a Tr, transition close to is only possible if I ' ^ ' _i 

IQ2J [ui-l<cGi 

(73) 

where c is a constant. Based on these considerations we find that establishing the propo- 
sition amounts to proving that for generic sequences the conditions in ([73]) are satisfied 
only for finitely many I as I ^ 00. For simplicity we focus on the first case of (j73p . i.e., 
we investigate the condition < c G~[ ; it is straightforward to adapt the results to the 
second case. 

Since asymptotic sequences ASr shadow sequences of transitions Sr, we can adopt the 
concepts introduced in Section [71 in particular, to obtain a more transparent formulation 
of the problem, we make use of the concept of eras: we view ASq- as a succession of 
eras [(/in)a, (^out)a] (a = 0, 1, 2, 3, . . .), and denote by G (0, 1) the value of at the 
beginning of era a, i.e., = {u™)^^ . Furthermore, we define 5a to be the value of the 
function on the r.h.s. of condition (j73p at the beginning of era a, i.e., 

where Qa is the era growth function defined in Section [9l Without loss of generality we 
assume the initial value of the sequence to be sufficiently small, i.e., < e. Based on 
these concepts, the first case of (|73p translates to the following necessary condition for the 
occurrence of transitions: 

An occurrence of a Trj transition close to T2 in era a is only possible if >Ca < 5a ■ (74) 

To establish the proposition we show that for generic sequences the condition < 5a is 
satisfied for a finite number of eras only. 

The arguments we present in the following are stochastic in nature. Our basic assump- 
tion is that we can model the sequence {>Ca)a&i by a random variable x on (0, 1) with a 
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well-defined probability density. The justification for this assumption and its origin are 
discussed in the subsection 'Randomized Kasner sequences' below. For the subsequent 
computations we employ the density w{>c) that describes the stochastics of Kasner se- 
quences; however, the basic qualitative results are largely independent of the concrete 
form of the probability density, see the remarks at the end of this subsection. 

Since x is a random variable, the sequence (xa)aeN represents an infinite series of trials 
with X. The sequence {5a)aeN, on the other hand, represents a sequence (0, 6a) of intervals 
whose length decreases as a ^ oo. The condition (j74p thus reads: the number of eras that 
can contain transitions coincides with the number of 'hits' x^ G (0, 6a) as a ^ oo. 

In Appendix [G] we demonstrate that the behavior of the sequence ((5a)aeN is decisive for 
the 'hit rate': 6a ~ leads to a diverging expectation value for the number of hits as 
a ^ oo. In our case this would lead to the conclusion that an infinite number of eras can 
contain transitions (although the recurrence intervals between transitions would 
become increasingly large). A sequence 6a ~ a~^, on the other hand, generates a finite 
expectation value; in that case only a finite number of eras is expected to contain 
transitions. 

The present situation is more complicated since (5a)aeN is not given explicitly, but is 
instead determined by a stochastic process. We have the relation 

6a+l = CQ^I^ = cQa^Qa^ = Qa^ 6a , 

see Section [9l and by using ([M]) and ([65]) , 

A 3 3xg 

Equivalently we can write 

log 6a+l - log 6a < log 3 + log X^ - log(l + X^ + X^) . 

Employing the probability density w{>c) for x, see ()50p . we find an estimate for the ex- 
pectation value of log 6a+i/ 6a, 

(log(5a+i -log^a) < log3+ / ?i;(x) (logx^ - log(l + X+ x^)) dxRi -1.77 < -1 , 

Jo 

see also Accordingly, for a S> 1 we expect 

(log 6a) < -a + log 6o < -a + log e . 

Consequently, the condition ([Till , i.e., Xq < 6a (or, equivalently, logx^ < log 6a) translates 
to the following statement: 

An occurrence of a transition close to T2 in era a is only possible if x^ < ee~° . 

Based on these considerations we can compute the probabilities for 'A^ hits' (cf. Ap- 
pendix[G]). For a 3> 1 the probability that x G (0, ee~") is given by 

P(x e (0, ee-")) = / t(;(x)dx =: 6a ~ ew(0)e-" . 
Jo 
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For a series of A trials we find that the probabihty P'^(O) ('no hit', i.e., the probabihty for 
the scenario x [0, ee"""] Va = 1, . . . , A) is given by 



A 

logP^(O) = log(l -Sa)^J2i- ^0)6-'^ + 0(^'e-2-)) ^ -ew{0) 



1 



.-A 



a=l 



a=l 



Similarly, the probability for 'one hit' is 

A 



a=l ^ "a 



P^(0) 



£W 



a=l 



P^(0) 



-Al 



ew{0)- 



and the probability for 'two hits' in a series of A trials, P'^(2), is given by 



P^(2) 
P^ 



a<b 



5a 



l-5al-5i, 



5a 



I -5a 



E 

a 



5a 



I -5a 



< 



1 r 



ew{0)- 



Letting A — > oo we obtain 



P°°(0) 



1 



P°°(2) < 



(75) 



where p = et(;(0)(l — e" 
hits is finite: 



p, P°°(l)<p, 
We conclude that the expectation value of the number of 



oo 

(#hits) = E^p°°w < E^T =^''^ < °° 

i=l 



Since the 'number of hits' corresponds to the number of eras that contain transitions, 
we have shown that the number of eras that contain transitions is expected to be 
finite and thus we have established the proposition. 

It is important to stress that the arguments of the proof rely on the assumption that the 
sequence ASr is generic. By stochastic methods it is impossible to exclude that there exist 
(non-generic) sequences that contain an infinite number of (or T/Vj, T];<i^^^, '^RiR2R3) 
transitions. 

We conclude this subsection with the observation that the above results are to a large 
extent independent of the concrete form of the probability density w{>c); e.g., it is sufficient 
to assume that w^O) is finite to obtain the result (j75p . This is essentially due to the rapid 
decrease of 5a, which is an exponential decay that stems from the (era) growth function. 

Let us finally summarize the results of this subsection: Bianchi type VIq and VIIq orbits 
are excluded from generic asymptotic sequences ASq, therefore a generic asymptotic se- 
quence ASo is an asymptotic sequence of transitions ASr- Generically sequences ASr 
do not contain transitions that involve an excitation of the variables i?2 and i-e., the 
transitions T^^r^, Trj, TJv2! and Tr^r^r.^ do not appear. The arguments leading to the 
stochastic exclusion of R2 and N2 are based on the stochastic properties of asymptotic 
sequences; next we discuss the stochasticity of asymptotic sequences in more detail. 
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Randomized Kasner sequences 

In Section [7] we have seen that generic sequences of transitions Sq- are associated with 
Kasner sequences (m;);^^! which are determined by the recursion (03]), where the initial 
value uq Q. In Section [8] we have presented a stochastic analysis of Kasner sequences: 
the collective distribution of all sequences (or, equivalently, (xj)^^^ = 

possesses a probabilistic description in terms of a probability density w. In the context 
of asymptotic sequences of transitions ASr-, however, we do not deal with exact Kasner 
sequences; in the following we investigate the resulting consequences for the probabilistic 
aspects. 

An asymptotic sequence of transitions ASt generates an approximate Kasner map, i.e., 

~ ~ ^ if > 2 

[(u; - 1) ^ if l<ui<2. 

Since asymptotic sequences ASr shadow sequences Sr to an increasing degree of accuracy, 
the errors in (j76p converge to zero as r — > oo. We presume that the nature of these errors is 
random, which is due to the fact that the approach ASr — > Sr is effectively unpredictable. 
This random element in the approach ASr — > Sr is the key to a probabilistic description 
of an approximate Kasner sequence (j76p : while an exact Kasner sequence might not be 
described by the collective probability density w[x) (or alternatively w{>c)), the random 
errors in (j76p smear out the probability distributions of exact sequences over the entire 
collective; we thus expect the sequence (>fj)jgN associated with an individual generic 
randomized asymptotic Kasner sequence to be described by the probability density w{>c) 
of Section [8l In brief, we can regard (i76|) as a randomized asymptotic Kasner sequence 
which admits a probabilistic description in terms of the density w{>c). 

The errors in (j76p are not uniform in u (i.e., independent of the value of u): the Kasner 
map u 1-^ u — 1 (or u {u — is relatively robust when u takes a value that is far 

from the extremes (which are u = 1 and u = oo), i.e., when n G (1 + e,e~^) for some e; 
in contrast, the errors in the Kasner map might be large when or (n — 1) is small. 
To see this we consider the (approximate) curvature transitions that generate (j76p . When 
u G (1 + e,e~^) the associated curvature transition connects two points on the Kasner 
circle K^ where the (angular) gradient of u (i.e., the gradient of u along the Kasner circle) 



is comparatively small, see Figure 1(b) A perturbation of the curvature transition, i.e., a 
deviation from the exact orbit, thus leads to a small error in the Kasner map. In contrast, 
when u > e~^, the associated curvature transition takes place in a small neighborhood of 



one of the Taub points, where the angular gradient of u is large, see Figure 1(b); hence, 
any deviation from the exact curvature transition orbit leads to a large error in the u- 
map. This effect is amplified by the local characteristics in the approach ASr Sr- 
we expect the accuracy of the shadowing to be smaller in the neighborhood of the Taub 
points (and the points Qq) than elsewhere. The dynamics in the neighborhood of the 
Taub points is rather complicated which is due to the fact that these fixed points are 
not transversally hyperbolic: the approach to the attractor is expected to be rather slow, 
and for fixed points close to the Taub points the actual dynamics is expected to display 
relatively large deviations from the linearized dynamics (more accurately described by 
'eras of small oscillations'). It is evident that this leads to an additional blow-up of errors 
in the n-map for large u. 
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Based on these consideration we expect that the Kasner map is not well-defined unless 
G (1 + e,e~^); only for this 'safe range' of u we can assume that ()76p holds. However, 
in the asymptotic regime the safe range increases, i.e., e — > as r — > oo, which is simply 
because ASr — > Sr- Suppose that the safe range grows rapidly enough. Then we can 
draw the following conclusion: the probability that the sequence (Mi)zeR ever leaves the 
safe range is zero. To see this, we use the same arguments that led to the stochastic 
exclusion of R2 and A^2: the 'number of hits' Xj- E (0, e^) (corresponding to the Kasner 
sequence leaving the safe range) has a finite expectation value. 

It is plausible that the increase of the safe range is in fact an exponential growth, since 
it is directly related to the accuracy of the shadowing. The decay rates of R2 and N2 
suggest that the overall accuracy of the shadowing behaves in the same way, which in turn 
leads to the required rapid growth of the safe range. Presuming the correctness of these 
considerations we find that generic asymptotic Kasner sequences are indeed described 
by (|76p . where the errors converge to zero (uniformly) as r — > 00. 

11 The billiard attractor 

In the previous section we have argued that for generic solutions the variables R2 and 
are not excited in the asymptotic regime, i.e., we have established the exclusion of 
transitions that involve the variables R2 and N2 from generic asymptotic sequences of 
transitions. Accordingly, the essential asymptotic dynamics of generic solutions X{t) 
is described not by the full oscillatory dynamical system ()35ap - (|35hp . but instead by a 
reduced system that is obtained by setting i?2 = and = 0. 

The billiard subset 

We define the billiard subset Oq as the invariant subset {R2 = 0) A (A'^2 = 0) of the 
oscillatory subset O, which leaves Sbiiuard = (5^1, 5]2, S3, -Ri, -R3, A''i) as the state vector. 

The system of differential equations on Oq is given by 



dr^2 

drRl 
drRs 

drNi 



2(1 



s2)Si - 2RI + iNf 



(77a) 
(77b) 





2(1 



S2)fli + (S2 



S3)i?l 



(77d) 





(77f) 



where 



i(s2 + s2 + si + 2i?2 + 2i?i). 



The Gauss constraint ( 35j ) reads 



1 _ y2 _ J_ 
^ ^ 12 



N? = 




and the Codazzi constraint ()35kp takes the form 



i?3 A^i = . 



(77h) 
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Due to the Codazzi constraint the bilhard subset consists of two invariant components: 
the component A''! = is the Bianchi type I subset Br^r^, which consists of K° U Br-^ U 
Br^ U Br^r^; the component A''i 7^ is the Bianchi type II subset, which is given by 
Br^ U Bni U Brj^Ni ■ The flow of the dynamical system on these subsets has been analyzed 
in Section [6l A depiction of the flow on the boundary subsets Br^ , Br^ (single frame 
transitions Tr^, Tr.^) and B^i (single curvature transitions T/vJ is given in Figure [3l for 
the flow on Br^r^ (double frame transitions Tr.^r^ ) and Br^n^ (mixed curvature transitions 
TriNi) see Figured! The analysis of the stability of the fixed points on K^, as regards 
Ob , is summarized in Figure [71 




Figure 7: Active triggers on associated with the billiard subset; note that R3N1 = 0. 

A sequence of transitions on the billiard subset, which we denote by Bq-, is an infinite 
concatenation of transition orbits on Og. Each transition % of Bq- is either a frame 
transition Trj, T^g, Tr-^r^, or a curvature transition TJvi, Tr-^Nx- It is understood that the 
final state of each transition Ti (which is a fixed point on K^) coincides with the initial 
state of the transition Ti+i. Clearly, billiard sequences Br can be viewed as a special case 
of the more general sequences of transitions Sr on O. 

The asymptotic suppression of multiple transitions 

The considerations of Section [10] can be summarized as follows: the dynamical evolution 
of generic solutions X is asymptotically silent and local, and the asymptotic behavior 
along generic timelines is governed by the dynamics on the SH silent boundary. At 
a spatial point x*, the asymptotic evolution of X{t) is thus represented by a solution 
S = {T,a/3, Aa, Nai3) of the finite dimensional SH dynamical system. The stable variables 
'S'stabie = (A^i2) A^23) A^3i) ^o) Vanish in the limit r ^ 0. Among the remaining oscillatory 
variables, which are Sosc = (Sq., -Rq,, A'"i, A''2), two more variables converge to zero: R2 
and A^2- This leaves the billiard variables SuiHnvd = C^a, Ri, R31 ^1) as the variables that 
capture the essential asymptotic dynamics. 
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The bilhard variables are constrained by the asymptotic condition 



which is the asymptotic Codazzi constraint that ensues from (j67p . In the same limit the 
variables Sbiiuard satisfy the billiard system (|77ap - (|77fj) and the Gauss constraint ( |77gP 



holds. The asymptotic evolution of Sbiiuard is thus governed by the dynamics on the 
billiard subset Og; in other words: the state vectors Sbiiuard of solutions X{t) shadow 
sequences of transitions Bq- on the billiard subset. Accordingly, in the asymptotic regime 
the orbit -S^biiiiard(''") can be partitioned into a sequence of segments, where each segment is 
associated with a transition on Og. We call this representation of the orbit S^biiiiard(''") an 
asymptotic billiard sequence ABr- Since a billiard sequence Br on Os is a concatenation 
of transitions of the types Tjvn '^iJu '^R-i, Tr^r^, 'TTViiJu an asymptotic billiard sequence 
ABt is a concatenation of approximate transitions, i.e., ABq- = {Ti)i^fq, where each % 
denotes an approximate transition in the present context. 

In asymptotic billiard sequences ABq-, however, multiple transitions are suppressed asymp- 
totically, i.e., transitions of the types Tp^^R^, Tn^Ri do not occur in the asymptotic regime. 
In the following we discuss this statement; however, we focus on the simpler case of Tpt-^Fii 
transitions. 

From the differential equations ()77dp and (j77ep we obtain 

dr log ^ = 3S2 and thus log ^ oc 3 / S2 dr , 

where we assume for simplicity and without loss of generality that i?i,i?3 > 0. In an 
asymptotic billiard sequence ABr the variables Ri and i?3 (as well as A^i) are characterized 
by oscillations between zero and finite maxima (which are bounded above via the Gauss 
constraint); in particular, Ri -/^ and i?3 7^ as r ^ cxd (as well as A''i 7^ 0). For Tr^ 
transitions, while i?i is of order unity, -R3 is small, so that Ri/R^ is large; conversely, for 
transitions, R1/R3 is small. 

Consequently, log(i?i/i?3) and thus J T,2dT is a function that oscillates in the range 
(—00,00). The central property of these oscillations is that the amplitude increases 
(rapidly) with increasing r, which is due to the increasing accuracy of the shadowing. In 
the asymptotic regime we can use a simple model for these oscillations which is based on 
the observation that the orbit ABr spends an increasing amount of r-time in a neighbor- 
hood of the Kasner fixed points, while each transition always takes a fixed At (depending 
on the transition). We obtain J Ti2dT = ^(S2)j(AT)j, where i £ N consecutively numbers 
the Kasner fixed points that the orbit passes, and where, accordingly, (Ar)j is the time 
the orbit spends in a neighborhood of the Kasner point i. The quantities (^2)1 range 
in (—2,2); in conformance with our assumptions we expect the distribution of the {"^2)1 
in (—2,2) to have a probabilistic description. Therefore, the sum ^(S2)j(AT)j can be 
thought of as arising from a stochastic process resembling a random walk ^Sj; in the 
present case, however, the step size Sj is not fixed, but itself random; most importantly, 
the step size increases with the number of steps taken. 

The growth of the step size is of crucial importance. Consider an interval (—a, a) (where a 
is small compared to the first step size). If Sj increases slowly with z, the probability that 
Sr=i •Si £ o) for some n (recurrence probability) is one. If the growth is sufficiently 
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fast, however, this probabiUty is zero; this is case, for instance, if the growth rate is 
geometric. The behavior of the quantity log(i?i/i?3) is determined accordingly: for every 
a there exists such that the average step size is large compared to a for all t > t^, 
the probability that log{Ri/R3) lies in the interval (—a, a) for some r > is expected 
to vanish; this is because the step size increases rapidly, which is in turn due to the 
expectation that the orbit ABj- shadows the attractor with a rapidly increasing accuracy 
(so that (AT)j grows fast). Consequently, in the asymptotic regime, the probability of the 
quantity log(i?i/i23) being in any given interval around zero vanishes. Since, therefore, 
either Ri is small and R3 large or Ri large or R3 small, our arguments suggest that 
multiple transitions are suppressed as r — > 00. 



The billiard attractor 



The asymptotic suppression of multiple transitions in the asymptotic regime entails that 
asymptotic billiard sequences ABr do not approach general billiard sequences Br on the 
billiard subset but sequences on the billiard attractor O/sj^. 

We define the billiard attractor subset as 

ObA = KO U Bn, U Br, U Br, , 
i.e., Oi3_A is the boundary of the billiard subset Ob- 

An attractor sequence of transitions At (or, for brevity, attractor sequence) is a sequence 
of transitions on the billiard attractor Osa- By definition. At is an infinite concatenation 
of single transitions of the type T/Vi, and Tr,, i.e., one type of single curvature 

transition and two types of single frame transitions; to get an intuitive picture of attractor 



sequences we simply refer to Figure 3(a) -3(c) 



The collection of our results leads to the formulation of the dynamical systems billiard 
conjecture: 

Conjecture. The asymptotic dynamical evolution of a generic timeline of a solution of 
Einstein's vacuum equations (expressed in an Iw as aw a frame) that exhibits a generic space- 
like singularity is characterized as follows: 



(i) It is asymptotically silent and local. 

(ii) In the asymptotic limit the essential dynamics is represented by an attractor sequence 
At on the billiard attractor O/sj^. 



The dynamical systems billiard conjecture can be viewed as the 'dual formulation' of the 
cosmological billiard conjecture by Damour, Hennaux, and Nicolai jl3j . In Section [T2] we 
will establish in detail the correspondence between the two approaches. In the remainder 
of this section we will give a brief summary of our analysis that has led to the dynamical 
systems billiard conjecture. 

Our previous analysis constitutes a derivation of the dynamical systems billiard conjecture; 
however, although a considerable part of our treatment meets the criterion of mathemat- 
ical rigor, some of our arguments are heuristic rather than rigorous (despite their being 
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mathematically convincing); tlicrcfore, the presented derivation of the billiard conjecture 
does not represent a rigorous proof but merely a first step toward a rigorous treatment. 
Nonetheless, apart from the fact that our analysis provides strong support for the billiard 
conjecture, we expect that many aspects of our considerations are fundamental for a deeper 
understanding of the asymptotic dynamics associated with generic spacelike singularities; 
in particular we believe that we have identified several issues that have to be taken into 
account in the pursuit of rigor. 

Let us therefore recapitulate the main steps in our derivation. Hereby we are guided by the 
uncovered hierarchical structure of (boundary) subsets: the essential asymptotic dynamics 
can be restricted successively to subsets of subsets, to boundaries of boundaries. 

A Assumption: Asymptotic silence and asymptotically local dynamics. Asymptotically 
local dynamics is defined by the requirement that — > 0, da{S,rp,lJjj) 0, 
{ra,Ua) 0, which leads to a spatial decoupling of the field equations. As a 
consequence, a generic solution X{t) of the Einstein vacuum equations at a generic 
spatial point approaches the SH part of the silent boundary = 0, which is 
spanned by the state vector S = (^^13, A^, N^is)- 

=^ B Assumption: Asymptotic dominance of the Kasncr states in the asymptotic evolution 
of solutions. Linear analysis of the Kasner circle K'-^ leads to the conjecture that the 
variables S'stabie vanish in the asymptotic limit r — 00, where Sgtabie = (-^a/3;^a) 
{a ^ (3). Hence the essential asymptotic dynamics is represented by the dynamics 
on a subset of the SH part of the silent boundary, the oscillatory subset O and 
its associated state vector Sosc = (^0,^0,^1,^2)- Orbits on O are heteroclinic 
orbits (in general: transitions) that can be joined to form infinite sequences. These 
sequences of orbits/transitions are accessible to a mathematically rigorous treatment, 
where stochastic aspects become important. 

C Assumption: Asymptotic shadowing. A generic solution X{t) shadows a sequence 
of orbits / transitions on O with an increasing degree of accuracy. Our analysis shows 
that we have R2 ^ and A^2 — ^ along generic orbits X (r) ; this implies that certain 
types of transitions that involve the variables R2, N2 are ruled out a priori (complete 
exclusion), while it can shown that the probability for the occurrence of any of 
the other transitions involving these variables is zero (stochastic exclusion). The 
asymptotic exclusion of R2 and A^2 entails that the essential asymptotic dynamics 
of X is described by the dynamics on the billiard subset Ob with its state vector 

'S'billiard = (^i , S2 , II3 , i?l , -R3 , A?^i ) . 

D Asymptotic shadowing of billiard sequences. A generic solution X{t) shadows a 
sequence of transitions on Ob (billiard sequence) with an increasing degree of ac- 
curacy. Our arguments indicate that the variables {Ri,R^,Ni) of ^billiard cannot 
be excited simultaneously. This leads to the conclusion that the generic asymptotic 
dynamics is governed by the flow on the boundary of the billiard subset, which is 
the billiard attractor Oba = K° U Bn^ U Br^ U Br^. The dynamics on this attractor 
set is represented by billiard attractor sequences (consisting of one type of curvature 
transitions and two types of frame transitions: T/v^, and Tr^, Tr^), which thus form 
the attractor for generic solutions X. In other words: the essential dynamics of a 
generic solution X is characterized by its approach to an attractor sequence At- 
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In this paper we have given heuristic arguments for A =^ B, but our main focus has been 
the other steps. Assuming that generic solutions X are asymptotic sequences of O-orbits 
(i.e., assuming shadowing of the flow on O) has inevitably led to the implication B ^ D. 
Since stochastic considerations are essential for parts of our considerations, we obtain only 
statements for generic solutions. The possibility exists that there are solutions violating 
one or the other properties; however, such solutions are expected to form a set of measure 
zero in the space of all solutions. 



12 'Duality' of Hamiltonian and dynamical systems billiards 

In Section [2] we have outlined the Hamiltonian approach to cosmological billiards, where 
we have followed Damour, Hennaux, and Nicolai [131 I14j . This approach is based on 
an analysis of the Hamiltonian Ti. in Iwasawa frame variables, which are: the diagonal 
degrees of freedom b", the off-diagonal variables M^i and the conjugate momenta. An 
asymptotic Hamiltonian TCoo is constructed from Ti. by taking the limit p = \/—bab'^ — > c>o 
and dropping all terms except for three terms that are identified as being 'dominant'. The 
asymptotic Hamiltonian Tioo is assumed to describe the generic asymptotic dynamics of 
generic solutions that exhibit a spacelike singularity. 

A direct consequence from the form of the asymptotic Hamiltonian Tioo is that the off- 
diagonal degrees of freedom are asymptotic constants of the motion, as is a momentum 
variable associated with a projection of the diagonal degrees of freedom — a phenomenon 
referred to as 'asymptotic freezing'; for a derivation of these results using the present dy- 
namical systems formalism and the results from the previous sections, see Appendix iHl 
Therefore, the non-trivial asymptotic dynamics is encoded in the diagonal degrees of free- 
dom, or, more precisely, their projections 7°^: the asymptotic dynamics at each spatial 
point is described as a geodesic motion in a portion of hyperbolic space that is bounded 
by sharp walls, see Section [2] and [13^ I14j. This picture prompts the terminology 'billiard 
motion' and 'cosmological billiards'; see Figure 8(i 



It is of interest to consider a Hamiltonian that is in a sense a link between Ti. and Tloo' 
the 'dominant' Hamiltonian TCdom- It is obtained by dropping the 'subdominant' terms of 
Ti. before the limit /) — > 00 is taken. Since the asymptotic Hamiltonian constructed from 
Tidom naturally coincides with Tioo, the asymptotic dynamics described by Ti and Wdom is 
the same. 

From Section [2] and [13\ [T^ we obtain the dominant Hamiltonian: 



"^dom = N 



where Ri := —V3, R3 := —{Vi + n3V2), and where Ni is a function of the variables 
given by Equation (j93ap . The non-trivial dynamics is represented by the diagonal degrees 
of freedom b"' (or alternatively by p and the projected variables 7"): the dynamics is 
described as a motion in a portion of Lorentzian space that is bounded by 'exponential 
walls'. As the solution approaches the singularity these walls become increasingly sharp 
with an infinitely high potential as limit; we refer to the remarks at the end of this section. 
The asymptotic dynamics is thus a free null geodesic motion in Lorentzian space inter- 
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rupted by bounces at the walls, which yields an alternative representation of cosmological 
billiards. 

The asymptotic dynamics as described by Hamiltonian cosmological billiards can be in- 
terpreted in terms of generalized Kasner solutions: Fermi propagated generalized Kasner 
solutions appear as straight lines (geodesies) in flat Lorentzian space in the 6° description, 
or, equivalently, as geodesies in hyperbolic space in the projected 7" description. A free 
motion between two bounces at walls thus translates to a (Fermi) 'Kasner epoch', a phase 
where the solution evolves as a generalized Kasner solution, see (j97p . The walls and the 
associated bounces, on the other hand, are of two kinds: (i) frame/centrifugal/symmetry 
walls and associated bounces; (ii) curvature walls and associated bounces. Bounces of the 
former type merely result in axes permutations of a Kasner solution; bounces of the latter 
type are Bianchi type II bounces: the Kasner state changes in accordance with a change 
generated by a Bianchi type II solution. 

The Hamiltonian approach emphasizes the dynamics of the configuration space variables 
(6" and the projected variables 7", respectively); therefore one may say that the Hamil- 
tonian picture yields a 'configuration space' representation of the asymptotic dynamics. 
In the following we will see that the dynamical systems approach can be viewed as a dual 
representation, i.e., a 'momentum space' representation. 

Consider the dominant Hamiltonian from the perspective of the dynamical systems ap- 
proach. The terms i?iexp(6^ — 6^), i?3exp(6^ — 6^), and iVi exp(— 26^) can be associated 
with i?3, and Ni by using the relations ()91bp and (|91cp . Accordingly, Wdom contains 
^{Rj + Rl) + lNf modulo a common factor. Recall in this context that the Hamiltonian 
is proportional to the Gauss constraint, cf. (|77gp. This indicates that there is an inti- 



mate connection between: (i) centrifugal bounces in the Hamiltonian picture and single 
frame transitions Tr^ and Tfjg in the dynamical systems picture; (ii) curvature bounces 
and single curvature transitions 7/Vi. 

When we calculate the Hamiltonian equations of Wdoih = associated with the non- 
trivial degrees of freedom b", adopt the appropriate time gauge, and then perform suitable 
variable transformations according to the formulas of Appendix [Aj we obtain: 

(i) The billiard system, i.e., the differential equations (j77ap - (j77f|) and the Gauss con- 
straint (|77gP; however, without the Codazzi constraint R3N1 = 0. 



(ii) Evolution equations for Ea = {Ei^ , , E^^) and the Hubble variable H, with all 
spatial frame derivatives and all variables X (and Ua,ra) set to zero except E^, H, 

and SbiUiard- 



The unconstrained billiard system (i) forms an independent coupled system of differential 
equations. The equations (ii) for H and Ea, on the other hand, decouple from the uncon- 
strained billiard system. These equations can be regarded as linear equations with time 
dependent coefficients, which are provided by the solutions of the billiard system. (Com- 
pare with the treatment of the auxiliary quantities in Section [9j) The evolution equations 
for Ea constitute the lowest order perturbation into the physical state space of the full 
equations for Ea in the vicinity of the silent boundary; together with the Hubble variable 
H they yield 6°^ according to Equation (1791 1) in Appendix [Aj The off-diagonal degrees of 
freedom can be treated similarly: the solutions of the unconstrained billiard system pro- 
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vide the time dependent coefficients for the lowest order perturbation equations. In this 
manner we obtain asymptotic freezing, see Appendix [Hi This result further strengthens 
the connection between the Hamiltonian and the dynamical systems approach. 

The asymptotic limit of the dynamics of the unconstrained billiard system is represented 
by attractor sequences Ar, see Figure [8(b)0 The duality of the Hamiltonian and the 



dynamical systems approach thus becomes apparent: one phase of free motion in the 
Hamiltonian billiard picture corresponds to a Kasner fixed point on the Kasner circle K^. 
The underlying reason for this is that the dynamics has been projected out of the self- 
similar Fermi propagated Kasner solutions by means of a conformal normalization that 
has yielded scale-invariant variables. The bounces at walls in the Hamiltonian billiards 
correspond to motion in the dynamical systems picture represented by the transition orbits 
(heteroclinic orbits) between Kasner points. The curvature transitions T/Vi ai"e non-scale- 
invariant solutions and correspond to curvature bounces in the Hamiltonian billiards; 
the frame transitions Tr^, Tr,^ rotate axes w.r.t. a Fermi frame; see Appendix [Cj The 
transitions appear as straight lines (in Euclidian space) when projected onto S^-space. 
Note that the transitions and the Hamiltonian wall bounces yield exactly the same rules for 
changing Fermi propagated Kasner states. (However, the Hamiltonian billiard only gives 
a consistent picture when we consider generic attractor sequences, which are characterized 
by irrational values of the Kasner parameter u. Rational values of the Kasner parameter 
u leads to motion straight into the sharp corners of the asymptotic billiard where the 
dynamics is undefined.) 

Changing from the Hamiltonian billiards to dynamical systems billiards, the notions 
'bounces at walls' and 'motion along straight lines' switch places. In the Hamiltonian 
picture, free motion is followed by a bounce at a wall, which again gives rise to a phase of 
free motion. In the dynamical systems picture, we observe 'bounces' at the fixed points 
on the Kasner circle K^; in fact, since the Kasner circle constitutes the boundary of S^- 
space, it can be viewed as a 'wall' in the dynamical systems picture. Between the bounces 
we have motion along straight lines in S^-space (transitions). In brief: the Hamiltonian 
'motion-bounce-motion-bounce' is translated to 'bounce-motion-bounce-motion' in the 
dynamical systems billiards. 

Since the variables are intimately connected with the variables tTq , it is natural to refer 
to the projected dynamical systems picture as a 'momentum space' representation of the 
asymptotic dynamics, which complements the Hamiltonian configuration space represen- 
tation. To compare the two pictures, see Figure [HI 

The correspondence between the dynamical systems and Hamiltonian pictures offers a 
possibility for mutual support. Let us give a few examples: (i) As can be inferred from 
Figures [5] and [8l the typical 'captured' T/v^-T^g oscillatory behavior (a small curvature 
phase) corresponds to the many oscillations in the upper right corner in the Hamiltonian 



billiard, see Figure 8(a) (ii) Figure [7] indicates that Tr^ transitions are less frequent 
than transitions: Ri and are unstable in three sectors of K^; however, i?i has 
to compete with A'^i in sector (132), which 'depletes' Tr^ transition at the cost of T/v^ 
transitions; T^g is unaffected by this. The correctness of these arguments is corroborated 
by the features of Lobashevski geometry close to the circular boundary and by the fact 
that the wall associated with Tfjg bounces is longer than the one associated with Tr^ 



^^A more detailed discussion of this statement will be given elsewhere. 
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(a) The cosmological billiard (b) The billiard attractor 



Figure 8: Figure (a) shows part of an orbit in terms of free Kasner (Fermi frame) motion 
and frame and curvature bounces. This represents a 'configuration space' projection of 
the asymptotic dynamics. The disc here represents hyperboHc space. Figure (b) shows 
part of an orbit in terms of single frame and curvature transitions, i.e., it shows a part of 
an attractor sequence. Note that the solution does not quite return to any of the Kasner 
points it has 'visited' before. This description represents a 'momentum space' projection 
of the asymptotic dynamics. The circle here is the Kasner circle K^. The dashed lines 
correspond to the two possible single transitions that are possible at this stage; which 
one is realized depends on initial data. This corresponds to that free motion in a given 
direction in Figure (a) may either lead to that one hits the wall associated with Ri (the 
short wall) or Ni (the curved wall). 

bounces, (iii) Consider the Hamiltonian T^dom and perform a projection of the dynamics 
onto spatial constant volume slices in 6"-space. The walls are 'exponential walls' whose 
contours are of a triangular shape where the corners are smoothed out. Toward the 
singularity (i.e., with shrinking volume) the triangle increases in size, where simultaneously 
the triangle's sides become increasingly straight and the corners increasingly sharp (relative 
to the size of the triangle), see for similar SH examples. Particularly the corner shape 
is sensitive to the additive effects from the different walls potentials. The increasing size 
of the potential contour triangle and the comparatively sharper corners suggests that it 
becomes increasingly less likely that the 'billiard ball' hits the corners. This geometrical 
feature suggests that, statistically, multiple transitions are asymptotically suppressed, in 
agreement with our quantitative analysis. 

It is of interest to also compare the dynamical systems picture with previous (Hamiltonian) 
treatments of so-called asymptotic velocity term dominated dynamics; this is done in 
Appendix [H 
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13 Models with symmetries in an Iwasawa frame 

In this section we give a brief overview on models with symmetries. Such models are of 
obvious significance for the study of cosmological singularities; since the equations simplify 
we can hope to obtain more substantial evidence and perhaps even prove theorems on the 
asymptotic dynamics of solutions. This in turn can be expected to shed new light on the 
heuristic statements concerning the generic case without symmetries. 

Gi models 

By definition, cosmological Gi models possess one spacelike Killing vector; let us choose 
this Killing vector to be d^i . The Iwasawa decomposition is compatible with Gi symmetry 
and we are therefore able to choose an Iwasawa frame that is symmetry adapted; we set 
da = Ead^t, where Ea is a lower triangular matrix, see Section [3l accordingly, 

di = Ei^d^i , 82 = E2^d^i + E2^d^2 , 83 = Es^d^i + ^3^5^2 + E3^d^3 , 

where Ea^ as well as H and M are independent of x^. Consequently, when / is a function 
that is independent of x^, it follows that 

dif = 0, d2f = E2^d,2f , ^3/ = {E3^d,2 + E3^d,^)f . 

This in turn leads to a decoupling of the Ea^ equations from the system (fT9]) - (f23]) . since 
Ea^d^i disappears from the r.h. sides of the evolution equations and the constraints. Note 
also that ()20al) and (I21al) do not mix the components E^^ with Ea'^ and -Eq,^. In addition, 
using the expressions in Appendix [Aj the symmetries are seen to further reduce the state 
space (I24b|26p : since 

N22 = N23 = Ai = Ui = ri = 0, 
the state space associated with the reduced coupled system of equations is given by 

X = {E2^, £"3^, E'a^, Sq, Nia, A2, A3) ; 

in addition U2, U3,r2,r3 are non-zero for most temporal gauge choices. 

Apart from the above generic Gi case there exists special Gi models. The hypersurface 
orthogonal Killing vector case, also referred to as the polarized case in a U{1) context, 
see |38) and references therein, is given by setting E2^ and £3^ to zero, so that the metric 
takes a diagonal form. Using the expressions in Appendix [Aj it follows that A'^n = R2 = 
R3 = 0. In this case, the sectors (231), (213), and (123) on the Kasner circle K° become 
stable, which in turn implies that solutions asymptotically approach a Kasner state; for 
proofs, see [38| [39]. 

G2 models 

Models with two commuting spacelike Killing vectors are known as G2 models, see, 
e.g., [T^. An Iwasawa frame that is adapted to the symmetries is obtained by letting 
di and 82 be tangent to the group orbits. Accordingly, E'^* and H, J\f are functions that 
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are independent of and and thus depend on and only. For any / 
we get 

dif = 0, d2f = 0, d3f = Es^d,3f; 
furthermore, as follows from Appendix E], 

N22 = Ni3 = N23 = N33 = Ai= A2 = Ui = U2 = ri=r2 = 0. 

In addition one can exploit the remaining freedom in choosing the frame at a spatial point 
to obtain R2 = 0, and then the Codazzi constraint and the evolution equation for R2 force 
R2 to be identically zero everywhere [40j. Using a symmetry adapted Iwasawa frame thus 
leads to a decoupling of all Ea^ equations except for the equation for £'3^. Accordingly, 
the reduced state space for the problem is represented by the state vector 

^ = {^3"^, Sq, -R3, A^ii, A^i2, A3) ; 

in addition f/s, ra are non-zero for most temporal gauge choices. If one in addition chooses 
the timelike separable area gauge |41j . so that the area density of the G2 symmetry orbits, 
i.e., A = (ei^e2^)~^, is a function of x^ only, then ^3 = r3, as follows from (I21ap . 

There exists several special classes of G2 models: 

(i) Hyper surf ace- orthogonal models are characterized by the restriction A'^n = R3 = 0. 

(ii) Orthogonally transitive models are obtained by setting Ri = 0. 

(iii) The diagonal case is given by A'^n = Ri = R3 = 0; this case is also referred to as 
the polarized case in literature dealing with a topology |42l I43j . 

(iv) Plane symmetric models are determined by A^n = = -Ri = -R3 = 0, Sn = ^22- 

In [22] slightly different variables were used to discuss the general case with T^-topology; 
these can be obtained as follows: firstly, note that Ua = Ua/H—ra , A^ = aa/ H+ra^ where 
the acceleration and spatial connection coefficients associated with g, respec- 
tively. Secondly, to identify the present frame with the one used in [22] we permute both 
frame and coordinate indices from the present ones to the ones used in [22] according to 
1^3, 2— >1,3— >2. Then, in addition, we make the few scalings and notational changes 
that are needed in order to obtain the variables (£^1^, r, [/, S+, S_, Ex , ^2) A^-, A^x ) used 
in [22] in terms of the present ones; £3'^ is replaced with £1^, which was used in |22j . while 
the other variables are related as follows: 

r = r3 = A3 , U = r3 + U3, (78a) 

^+ = ^(^11 + ^22) , 5]_ = ^^(Sii — S22) , = ~:^-^3 1 (78b) 

S2 = -^i?i, ^^- = i73^n' ^x = ;)5Ar,2. (78c) 

For our present purposes it is of interest to note that numerical experiments performed 
in the general G2 case in [22| gave evidence for the claim that multiple transitions be- 
come increasingly rare and that the billiard attractor gives a correct generic asymptotic 
description. However, the results in the present paper suggest that it would be useful to 
numerically check our more detailed statements about asymptotic suppression and decay 
rates (as well as freezing, as discussed in Appendix |H]) ; the same also holds for the next 
class of models. 



= /(x°,x^) 
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Generic Bianchi type VI_i/9 models 

The evolution and constraint equations for SH cosmologies are given by setting all terms 
involving spatial derivatives to z6ro, which, includes setting J/q, — — Tq,, see also Sections [3] 
andm This implies that the conformal Hubble normalization becomes equivalent with the 
usual Hubble normalization used in, e.g., [I2]. While the Iwasawa gauge is compatible 
with the symmetry adapted frames of Bianchi types I-VH, since they admit a subgroup 
of two spacelike commuting Killing vectors, it is incompatible with the symmetry adapted 
frames of Bianchi types VHI and IX. The most general model with a symmetry adapted 
Iwasawa frame is the general Bianchi type VI_i/9 model, also sometimes referred to as 
the exceptional type Ylh model. 

The asymptotic dynamics of generic Bianchi cosmologies of type VI_i/9 has been studied 
in [36]. The compatibility of these models with the Iwasawa gauge is responsible for their 
particular significance in the context of generic inhomogeneous cosmologies: the general 
type VI_i/9 case is the only SH vacuum case where models possess the same attractor 
toward the singularity as generic inhomogeneous cosmologies in an Iwasawa frame; in 
particular they exhibit oscillatory asymptotic behavior. The type VI_i/9 models admit 
a G2 subgroup; adapting the spatial Iwasawa frame to the three-dimensional symmetry 
surfaces so that the first two vectors are adapted to the G2 subgroup implies that the 
same quantities as in the G2 case become zero, however, the extra symmetry forces some 
additional quantities to be zero as well. The Codazzi constraint with a = 2 yields Ri{Ni2 — 
3A3) = 0, but, in contrast to other type VI/i models, h = —1/9 implies, with the present 
conventions, that N12 = SA^; hence Ri is not forced to vanish since the Codazzi constraint 
is automatically satisfied, see, e.g., [12]; the general type VI_i/9 case is thus characterized 
by nonzero off-diagonal shear components Ri and R3. It follows that we can choose 

S = (Sii, S22, Ri, R3, Nil, As) , 

as the reduced state space variables. 

Using the variables (fTHll with the opposite sign^ for S2 and Sx, and Va = Ua = 0, 
yields the system given in [36], where a numerical investigation indicated that multiple 
transitions become increasingly rare and that the billiard attractor gives a correct generic 
asymptotic description, in agreement with the present general analysis. 

14 Gauge considerations 

It is a matter of course that there exists many useful spatial frame choices; while our 
present considerations are based on the choice of an Iwasawa frame, other prominent 
gauges are: the frame used in UEWE [2Dj; the SO{3) choice used by Benini and Montani 
as the starting point for a Misner/Chitre billiard analysis [44]; the Fermi choice {Ra = 0). 
Each of these gauges exhibit some computational advantages over the others in special 
contexts, and each has its physical merits: e.g., Fermi observers associated with Fermi 
frames do not get dizzy. But it is only in special situations, when one considers an 
important class of spacetimes that share a special, often global, feature, that one can 

^■'This is because the opposite sign was used in [121 136) for the definition of Ra- 
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argue that a particular frame choice is preferred. It is probably safe to conclude that no 
single frame choice will be optimal as regards all features one may be interested in; e.g., 
Iwasawa frames are incompatible with the symmetry adapted SH frames of Bianchi type 
VIII and IX, which are historically the prime examples of 'Mixmaster' dynamics. 

The asymptotic causal structure associated with asymptotic silence indicates that local 
and quasi-local aspects, and not global ones, are most important when trying to under- 
stand the physics of generic spacelike singularities in GR. Since this suggests that there 
is no preferred choice of frame, it is natural that results should be translated to a frame 
independent description in order to separate gauge features from physics. In particular, 
the attractor should be characterized within the framework of a spatial frame invariant 
description of the union of the silent Kasner and Bianchi type II subsets; it follows from 
this that we are interested in quantities that do not involve Rf^ • 

Let us therefore characterize this subset in a frame invariant way. It is described by the 
equations 

G'^ = U^ = r^ = A^ = iV„^E,^^E^^(iV°T - E"^) = Dn = An = 0, 

where we have introduced the following quantities: 

G'^ := E^' Ep^ , Dn := det(iV„;3) , := i [N^pN'"" - [N^'^f] . 

The quantity G^^ is the conformal Hubble-normalized contravariant three-metric0, see 
plj : if N'^^ is of rank two, then = 0; if A^"^ is of rank one (as it is for Bianchi type 
II), then = Atv = 0. The quantity N^pT.^^T.'^'' {N'^^ - S"^) is the square of the 
expression Ca^^ Njjs^-y^ , which is associated with the asymptotic Codazzi constraints. 

To explore the frame invariant description of the attractor it is of interest to consider the 
following Weyl spatial and spacetime scalar invariants 

c2 1 c caB -1/2 1 -7/ ijctl3 

= 6 '^apo , n = Q rLafSiL , 

f^ahcd /~i */-iabcd 

where {£af3,'Ha/3} = {Ea/3 , Ha/3} / H'^ are the electric and the magnetic part of the Weyl 
curvature tensor. (Here, Eap, Hai3 refer to the orthonormal frame components of the Weyl 
tensor of g while Sa/s^'Ha/s are the orthonormal frame components of the Weyl tensor of 
G, or, equivalently, the conformal orthonormal frame components of the Weyl tensor of 
g)- 

The attractor consists of a sequence of Kasner states connected by Bianchi type II solutions 
(i.e., attractor sequences). To formulate this in terms of curvature frame invariants we 
proceed as follows: we first give the electric and magnetic parts of the Weyl tensor for the 
Bn, subset: 
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iA^i (S3 + 2Si 



^^Incidentally, this suggests that it is this object that is natural to use if one wants to pursue a metric 
approach to the Einstein field equations and generic spacelike singularities. 
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The equations (|39p then allow us to express the above equations in terms of the frame 
invariant and subsequently we can compute the curvature scalars in terms of C, i-e., we 
obtain a frame invariant description of the scalars. However, the resulting expressions are 
quite complicated and we will refrain from giving them explicitly. The Kasner states yield 
Wi = 27n^(l + u)'^/{l + li + li^)^ and W2 = 0; Wi is monotonically decreasing in u with 
the range (0,4) with Wi = when u = 00, i.e., for the Taub state, while Wi reaches its 
maximum value Wi = 4 for the non-flat LRS Kasner state u = 1. 

The scalars can be expressed implicitly in r-time for a Bianchi type II transition since it 
is possible to integrate (f39l) : we obtain 

(c-c-)^(c+-c)^r'/' 



r(C) = In 



In Figure 9(a) we plot a typical attractor sequence; note that a change of era occurs. We 



then plot the S"^, Ti.^ and Wi, 'W2 scalars against each other in Figs. |9(b)| and 9(c) , and 



obtain a frame independent description of the same orbit in terms of curvature properties. 

Apart from spatial gauge considerations, it is of interest to ask oneself: to what degree do 
asymptotic silence and local dynamics for a generic spacelike singularity depend on the 
temporal gauge one is using? Note that asymptotic silence has been defined gauge invari- 
antly, and hence the issue is whether a generic spacelike singularity is asymptotically silent 
(which our results indicate) and to what extent asymptotic local dynamics is a temporally 
gauge robust feature for a generic spacelike singularity. The asymptotic causal structure 
associated with asymptotic silence suggests considerable temporal gauge robustness, and 
we expect a large class of time choices to be compatible with asymptotic local dynamics. 

We first note that if one considers a foliation associated with a time choice, and if it is 
possible to choose a parametrization of this foliation so that N is bounded and greater 
(or smaller) than zero asymptotically toward the singularity, then the billiard attractor is 
a local attractor in the full state space; the difference between using such a time variable 
and the local reparametrization that uses r only amounts to an unimportant conformal 
factor on the right hand sides of the billiard equations for each fixed x*. We then note that 
the billiard attractor is a local attractor if one uses any temporal choice compatible with 
the asymptotic local dynamics condition, which, in particular, implies that the timelike 
congruence is asymptotically conformally geodesic since — > (which includes the syn- 
chronous choice used by BKL). However, presumably it is possible to choose a time gauge 
such that e.g., 7^ 0, which is hence non-compatible with asymptotic local dynamics, 
even though the spacetime may exhibit a singularity for which there exists a large class of 
time choices for which the dynamics is asymptotically local. 



15 Concluding remarks 

In this paper we have considered the dynamics of solutions of the Einstein equations in 
the asymptotic limit toward a generic spacelike singularity: the asymptotic behavior of 
solutions is represented by the billiard attractor and the associated attractor sequences. 
We have derived the billiard attractor based on the dynamical systems formulation of 
the Einstein equations in the conformal Hubble-normalized Iwasawa frame approach; in 
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(a) A typical attractor sequence. 





(c) Wi - m 



Figure 9: In Figure 9(a) a typical attractor sequence is shown, which is associated a 



sequence of Kasner epoch characterized by the values of u given by {uo,ui,U2,U3,U4} = 
{3.18,2.18,1.18,(1.18-1)"^ = 5.56,4.56}. Figures |9(b)| and |9(c)| display the corresponding 
Weyl scalar properties (note that the scalars are constant during frame transitions); the 
bullets pinpoint the Kasner epochs. The dashed line in |9(c) indicates a new era and the 
arrows show the direction of increasing time toward the singularity. 
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addition we have established the 'duahty' of the presented framework with the Hamiltonian 
approach to cosmological bilhards of Damour, Hennaux, and Nicolai |13j . 

The cornerstones of our derivation of the bilhard attractor have been (i) the identification 
of a hierarchy of state spaces and invariant subsets; in the asymptotic hmit, dynamics is 
restricted to subsets of subsets, to boundaries of boundaries, descending from the fuh state 
space via the SH silent boundary and the oscillatory subset down to the billiard attractor; 
(ii) a thorough analysis of sequences of transitions; in particular, we have introduced the 
concepts of small and large curvature phases that appeared prominently in our consid- 
erations; (iii) stochastic analysis; we have investigated the probabilistic nature of Kasner 
sequences and we have discussed the effects leading to randomization; randomized Kasner 
sequences and randomized asymptotic attractor sequences are fundamental concepts in 
our analysis; (iv) the computation of decay rates; in connection with probabilistic aspects 
we have given the decay rates of several quantities; this decay has turned out to be the 
underlying reason for the restriction of the dynamics to subsets and boundaries of subsets. 
The presented derivation of the billiard conjecture is not mathematically rigorous, but 
depends on arguments that are merely heuristic — despite their being convincing. How- 
ever, we are positive that in any endeavor aimed at obtaining proofs, several concepts and 
methods introduced in this paper will play a prominent role. Let us elaborate on this. 

Generalizing the concept of an era has played an important role in our analysis. The 
partitioning of sequences into small and large curvature phases has shown that small 
curvature phases dominate over large curvature phases from a stochastic point of view; 
hence, in the asymptotic evolution of solutions the phases when the solution is close to the 
Taub solution are crucial — a fact that might not have been appreciated enough previously. 
Note that the division into small and large curvature phases has featured prominently in 
the more technical computations, see Appendix iFl for details. 

Asymptotic (attractor) sequences give a precise meaning to the notion of 'piece-wise ap- 
proximations' by BKL. The dynamics of a generic timeline of a solution with a generic 
spacelike singularity is expected to asymptotically shadow the attractor, and hence the as- 
sociated orbit can be partitioned into segments where each segment can be approximated 
with an increasing degree of accuracy by a heteroclinic orbit on the attractor. However, 
due to 'errors' associated with the approach toward the attractor, the solution does not 
exactly follow the heteroclinic orbit structure on the attractor; instead these increasingly 
small errors lead to an approximation that might be subsumed under the notion of a 
randomized sequence of heteroclinic orbits. 

Randomization constitutes an important ingredient in several contexts. The underlying 
structures that allowed one to obtain mathematical proofs about the attractor of the 
diagonal Bianchi type IX models in a Fermi propagated frame \18\ [T9] are specific for 
these models and are not available in other cases; since it is such cases that are relevant 
for the present general scenario, the diagonal Bianchi type IX models are misleading. Our 
present analysis suggests that one has to know the (asymptotic) history of a solution to 
unravel its asymptotic features; this causes a dilemma since this requires that one finds the 
solution, which seems unlikely. However, randomization makes it possible to stochastically 
examine the cumulative effects of small and large curvature phases and this allows one to 
estimate decay rates and give a description of what is going to happen generically. In our 
opinion statistical analysis will be an essential new ingredient in future proofs. 
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In this paper we have obtained decay rates associated with asymptotic suppression and 
freezing (see Appendix |H|) which will hold in the neighborhood of the billiard attractor 
as well. Since decay rates are likely to be ingredients in any proof, it is important that 
numerics not only check that the dynamics of a generic spatial point approaches the 
attractor, but that the attractor is approached in the way described by the decay rates; 
therefore, these results offer new input for numerical investigations. 

As we have seen, the asymptotic dynamics is inevitably restricted to subsets of subsets, to 
boundaries of boundaries, which underlines the importance of the hierarchical structure 
of invariant sets associated with asymptotically local SH dynamics. However, there also 
exists a hierarchy of invariant subsets that is based on the number of Killing vectors 
the spacetime admits. This symmetry-based hierarchy is of interest since it is associated 
with different levels of technical difficulty; taking highly symmetric, and thus technically 
simple, cases as a starting point, one can analytically and numerically explore, step- by- 
step, the cases of decreasing symmetry by making use of the gained experience; since 
several structures that are relevant to the non-symmetric case already appear in cases of 
special symmetries, this is likely to lead the way to the generic case. The present work 
shows how the Hubble-normalized Iwasawa frame approach naturally provides a framework 
for dealing with this type of hierarchy; the formulation comprises different symmetry levels; 
note, however, that the Iwasawa decomposition is natural only when one has commuting 
spacelike Killing vectors, and hence it is not appropriate for Bianchi types VIII or IX. The 
simplest model with an oscillatory singularity as described by the present formalism is 
the general Bianchi type VI_i/9 model. The next level of difficulty as regards oscillatory 
singularities is the general case with two commuting spacelike Killing vectors; then the 
general case with one spacelike Killing vector; and finally the general case with no Killing 
vectors. Incidentally, by imposing geometric restrictions on the Killing vectors such as 
hypersurface orthogonality, see, e.g., [ISdS], one can obtain relatively simple classes of 
models with non-oscillatory singularities. Exploiting this outlined hierarchy of symmetries 
might prove to be the key to a complete understanding of the nature of generic spacelike 
singularities. 

A dynamical system that is of particular interest in our analysis is the system ()77p : when 
we do not impose the Codazzi constraint ()77hp we refer to this system as the unconstrained 
billiard system, see also the discussion in Section [12] in connection with the 'dominant' 
Hamiltonian. We expect the unconstrained billiard system to possess the same generic 
asymptotic behavior as the full system X in the neighborhood of an asymptotically silent 
singularity; hence, in our approach, it is this system that is of central importance for the 
asymptotic dynamics, rather than, e.g., the Bianchi type IX system. We will provide a 
thorough analysis of this system and its properties elsewhere. 

Another interesting issue we have put little emphasis on in the present paper is the question 
of consistency. It is suggestive to take the billiard attractor and the associated billiard 
conjecture as a starting point for further developments by inserting attractor sequences 
At into the relevant equations; using similar reasoning as done here for the suppression 
of variables, we will be able to derive more substantial justification for the claim that 
Ej ^ (a major ingredient for asymptotic silence) and instable 

when r ^ oo (cf. 

steps A and B in Section llip . In addition we will obtain further information on the 
approach of solutions toward the asymptotic limit. We will pursue issues of this kind in 
future work. 
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In this paper we have been concerned with the vacuum GR case, however, the approach 
and methods we have developed will be applicable to the case when matter is included 
and also to other theories (see, e.g., [131 [Hj a-^d references therein). As regards possible 
sources this naturally suggests a study of the influence of matter on singularities. As 
a simple example we refer to |46] . which gives an indication that, e.g., Vlasov matter 
behaves differently in some respects than perfect fluids. Of particular interest is the 
question of structural stability of generic spacelike singularities, especially since generic 
singularities seem to uncover the essential properties of matter. There are indications that 
a classification of the influence of matter on generic singularity structure naturally rests 
on (i) energy conditions and (ii) whether the effective propagation speed is less than the 
speed of light or not. There are also indications that suggest a subclassification based on 
the behavior of matter in the case of speed of light propagation: e.g., massless scalar fields 
and electromagnetic fields influence the generic spacelike singularity in different ways; 
does this motivate a subclassiflcation based on spin? In the sublight case it seems to be 
natural to base a subclassiflcation on the question of what features are affected by matter 
and what features are not, e.g., the Hubble-normalized energy-momentum tensor may go 
to zero, which leads to an asymptotic description of the geometry by vacuum solutions, 
but this does not necessarily mean that, e.g., the Hubble-normalized rotation of a perfect 
fluid tends to zero. Another issue is how, and if, matter influences the connection between 
generic spacelike singularities and weak null singularities in asymptotically flat spacetimes, 
see j23]. Questions and issues like these lead to a variety of possible research projects that 
could make our understanding more substantial. 

In this paper we have assumed asymptotic silence and asymptotic local dynamics, and our 
work has provided results and evidence that support the consistency of this scenario for 
generic timelines. However, this does not mean that there could not exist interesting phe- 
nomena associated with special timelines, indeed, we believe that a set of measure zero of 
timelines exhibit spike formation and recurring 'spike transitions' [22] that are associated 
with non-local dynamics. Even so, it still seems that asymptotic silence prevails and that 
asymptotic local dynamics play an important role — remarkably, spike transitions seem to 
be governed by variations of the Kasner map, which is associated with asymptotically local 
dynamics, hinting at further more deeply hidden structures. Asymptotic spike formation 
is associated with the unstable variables A'"i , Ri, going through a zero at a spatial point 
(where a zero in A^i yields a 'true spike', seen in the curvature, while zeroes in Ri, R^ yield 
false spikes, i.e., gauge effects not seen in the curvature; presumably, zeroes in N2 and R2 
do not play a similar role since these variables are generically suppressed in the asymp- 
totic limit). There are hints that spikes may interfere — destructively and constructively. 
Although uncertain, constructive spike interference seems to dominate, which, if correct, 
leads to asymptotic 'spike cascading' [I^. There are many other unresolved spike issues 
as well: are there spikes that undergo infinitely many recurring spike transitions? Where 
do spikes form and in what way — how does spike interference work? Do spikes asymp- 
totically form a dense set? Can spikes leave observational imprints in, e.g., the cosmic 
microwave background? These question are clearly interesting in themselves, moreover, a 
clarification of some of these issues is of considerable interest in the context of eventual 
generic singularity proofs. 

Ultimately one might ask oneself the question why generic singularities should be studied 
at all in a classical GR context? Firstly, there exists a regime between the Planck era and 
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the GUT era where GR is expected to hold and where the approach toward the singularity 
is presumably described by the dynamics toward a generic singularity (recall that one of 
the points of inflation is to 'erase the effects of initial data' and that before this erasure a 
singularity is presumably generic according to this line of reasoning). Secondly, black hole 
formation is associated with initial data reflecting the complexities of the real universe; 
one would hence also in this case expect generic spacelike singularities to play a role before 
one enters the Planck regime. Thirdly, the formation of generic singularities is associated 
with considerable structure, even in the case of spike formation: can this structure be 
used to asymptotically quantize gravity where it needs to be quantized, namely in the 
ultra-strong gravitational field in the neighborhood of a generic spacelike singularity 'j^ 
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A Iwasawa variables and useful equations 

In this section we derive in some detail the connection between the Hamiltonian Iwasawa 
variables of Section [2] and the conformal Hubble-normalized variables of Section [31 Let us 
begin by recalling that for an Iwasawa frame 

e° = exp(— 5")AA"j , ea* = exp(6°)A/'*a (no summation over a) , (79a) 

where M^i and Af^a are the upper triangular matrices given by ([5]) and Q, respectively; 
the off-diagonal components of AA'^j {-Af^a) are denoted by Ua (n^). Therefore, due to ([16]) . 
we obtain for the conformal frame components: 

E^'i = Fexp(-6")AA"i , Ej = H'^ exp{b'')M\ (no summation) . (79b) 

In the following we use the notation = -Eq" (no sum over a) and analogously = ea". 
Note that the frame variables have been chosen to be positive, from which it follows 
that also Ea > 0, since we take the cosmological model to be expanding toward the future 
and hence H > 0. Equation (1790 can be inverted easily: 

6° = log(e,) = log{HE^) , ^1 = ^ , ^2 = ^ , ^3 = ^ . m) 

£/2 £/3 

In order to make contact between the Hamiltonian variables and the conformal Hubble- 
normalized variables we first consider the Hubble scalar H. Since gij = e'^^e^-Sap, we 
have 

^ = exp(-6i -1?- b') = exp (- 6") ; 

^^For the exploitation of some of these structures in the context of quantization of special models, see, 
e.g., [Mils]- 
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accordingly, we obtain for the derivative (w.r.t. coordinate time x^): 
Since the derivative of is proportional to the expansion 9 according to 
this implies that 

Using the relation (|9aP between 6" and tTq,, where we recall that A'" = N^Jg^ we obtain a 
representation of H in terms of the momenta, i.e., 

Equation ([52]) suggests that we define a quantity A, 

K:=2H^=\y] (83) 

which we will use frequently in the following. In terms of the frame components, A can 
be expressed as 

A = 2H^ = 2{H^EiE2E^r\ (84) 
since ^ = (eieaea)-^ = R-"^ {EiE2E^)~^ . 

Inverting the commutator equations ()17|] yields the conformal Hubble-normalized variables 
in terms of M and Ea and their derivatives: 

U^=d^\ogM , (? = i S", ao^„\ (85a) 

j^a _ 1 ^a^-y j^f3 ^ SqE^^ , Tia/S = —E'^i 5o-E'(a* ^l3)-y ) (85b) 

= i E^ daEp' - \ E^, dpEj , 7V"^ = e"^^' d^Es' . (85c) 

Inserting (179b[) into the equation ()85ap for g yields 

Since 5^.0 = HNdo, Equation (fST]) takes the simple form do 6" = —3, which in turn 
implies = 3y^ from ([80|) . Using the expression ([82]) to compute BqU we eventually 

obtain 

<7 = 2-^A-%V (87) 

Note that Equation ([55]) reproduces the dimensional equation 9o-?^ = —(1 + Q')-H^ for H. 
Inserting (|79bp into (|85bp yields 

i?i = -ie^'-^'aon3 , i?2 = ^e'''-''' (-ngaoni + ^0^2) , i?3 = -\e'''-''' d^ni , (88a) 
E23 = +\e''~'"don^ , Si3 = ie'''-''' {-n^B^n^ + ^0^2) , E12 = +\e''-''don^ (88b) 
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for Ra and {a ^ P), and 
^11 = IdoY.b'^-dob', S22 = ^ao5]6"-ao5^ S33 = ^doY.b'' -dob' (88c) 

a a a 

for the diagonal elements of Sq,/?. The main observation is that the off-diagonal components 
of T,af3 are given by the Fermi rotation parameters Ra according to 

(5^23i ^31, ^12) = {—Ri,R2, —R3) , 
which is a fundamental property of the Iwasawa gauge. 



In order to express to r.h. sides of p8|) in terms of the momenta tTq, and Va we write ([9]) 
in the form 

(90a) 

doni = I e2(^^-''') {Vi + m V2) , ^0^3 = | e2(^'-''') , (90b) 

where we use N~^d^o = {A/2)do. A straightforward computation then results in 

Tict = Sao = 2 — A^"*^ tTq , (no sum) (91a) 

Ri = -A"^ exp(62 - b') P3 , i?2 = exp(6^ - 6^) P2 , (91b) 

i?3 = -A-1 exp(6i - b'^){nsV2+Vi) , (91c) 
where A is given in terms of the Hamiltonian variables by (|83p . 
Inversion of (|9ip yields 

TTo = A [2 - Sa] , dSEI) 

Pi = -A i?3 E2/S1 + n3P2 , P2 = A i?2 ^3/^1 , dSS) 

V3 = -ARiE3/E2, m^) 

where A is given by (j84p . 

To obtain the conformal Hubble-normalized variables and Na/s in terms of the Hamil- 
tonian variables we proceed analogously. Based on (j85cp we find 

Nq. = exp (—26") Ao (no summation over a) , 

A„^ = A-i exp(-(6- + 6^))/^^ (a / /? / 7) , 

^„ - r„ = A-i exp (-(6- + 6^)) /^+ (a / /? / 7) , 

where 

Ai = 2{d2ns - Bsfii + (ni + l){nidins - ^2713)) , (93a) 
A2 = -2ain3 , (93b) 
A3 = , (93c) 
h± = di{b^±b'), 5i = 9i, (93d) 

f2± = d2{b'±b^)Tdini, d2 = nidi + d2, (93e) 

/si = dsib^ ± b^) =F 92n3 + ni5in3 - 9in2 , <93 = "-2^1 + n3(92 + ^3 ; (93f) 
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as usual we denote the diagonal elements of N^p by N^- In particular we see that 

N3 = 

in the Iwasawa gauge. 

Finally we note that the momentum conjugate to A = log p reads 

vr, = Alog[(ei)(2"S^)(e2)(2-^^)(e3)(^-^=')]. 
where = HEa = exp(6'^). 

We conclude this section by investigating the equations for the conformal spatial frame 
components and the consequences thereof in some detail; these equations will be useful 
in Appendix [Hi From (|17ap and (j20ap we see 

doE^' = Fj^Ep' , where = q6£ - ^£ - e£^W . 

In Iwasawa gauge ()89p we obtain 

doEi ={q- Si)Fi doE2 ={q- ^2)E2 doEs = {q - ^3)^3 , (94a) 

where Ei, E2, E3, and Si, S2, S3, are again the diagonal components of E^ and Sc^/j, 
and 

doE2^ = 2R3E1 + {q- T.2)E2^ doE3^ = 2R1E2 + {q- S3)F3^ (94b) 

^0^3^ = -2R2E1 + 2RiE2^ + {q- S3)F3^ , (94c) 

for the non-zero off-diagonal components of Ea^. Using that 

doH = -{l + q)H 

we find for the derivatives of = HEq,: 

aoei = -(1 + Si)ei does = -(1 + S2)e2 ^063 = -(1 + S3)e3 . (95) 

Furthermore, from 

V5 = H-^ {EiE2E3)-^ , 

we obtain 

from a direct computation based on (j94ap : therefore, 

doA = do{2H^) = (2 - q){2H^) = (2 - q)A . (96) 
Evidently, this equation is consistent with ([83]) and ([ST]) . 
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B Deriving the generalized Kasner line element of BKL 

The Kasner subset on the silent boundary Ea^ = is characterized hy Ua = Ta = (SH 
condition) and = 1 — = Aa = N^p- We speciaHze to the Kasner circle K^ by setting 
Ra = which corresponds to setting = (a / /?). Therefore, for every point on K^, 
= diag[Si, S2, S3], where the shear variables Sq, (q = 1 . . . 3) are temporally constant 
spatially dependent functions; we set diagpn, S22; S33] = diag[3pi — l,3p2 — 1,3^3 — 1]. 

In order to derive the leading asymptotic expressions for H and we first note that 
q = 2 which is due to (p2al) . Equations (HH) and (fTTa]) . (f20a|) yield 

H = He^\ EJ = EJ e-^^^-P'^)^ (no sum over a) , 

which leads to 

e'^j = e"j e"^^""^ (no sum over a) ; 
recall that hatted objects denote temporally constant spatially dependent functions. 
In the synchronous gauge of BKL, i.e., N = 1, M = H , equation ([29]) leads to 

t = i+\H-^e-^^ . 

In addition, BKL also choose the synchronous time coordinate to be a simultaneous bang 
time function. In the present derivation this amounts to setting t = as an initial value 
condition, which can be accomplished by making a coordinate transformation. Conse- 
quently we obtain t oc e"^"^. With suitable redefinitions, i.e., e^i oc li, e^j oc mj, e^j oc n^, 
and a = t^^, b = t^^, c = t^^, the 'generalized Kasner line element' of BKL ensues: 

ds"^ = —dt^ + [a^ li Ij + rrii mj + nj] dx^dx^ . (97) 

Thus, instead of ad hoc assuming the above line element, and inserting it into the field 
equations in order to analyze its consistency, we have now derived it as the lowest order 
perturbation of K° on the silent boundary. 

Finally, note that an equivalent derivation can be performed in a Fermi frame instead 
of an Iwasawa frame: we could start with = 0, which would lead to a 4-dimensional 
ellipsoid of Kasner equilibrium points characterized by temporally constant values S^,^, 
and then make a temporally constant rotation and diagonalize this matrix. 

C Kasner solutions in a rotating frame 

It is customary to represent a generalized Kasner solution in a Fermi frame with constant 
diagonal shear variables: diag[Si, S2) ^3] = const. However, when one considers a rep- 
resentation in a frame that rotates w.r.t. the Fermi frame, one obtains time-dependent 
non-diagonal shear variables ^^^(t) that are related to through time-dependent rota- 
tions 0|^(t), i.e., 

(02' Sa/3 Ot] = diag(Si, S2, S3) = diag(3pi - 1, 3p2 - 1, 3p3 - l) = const , 



-68- 



C Kasner solutions in a rotating frame 



where the quantities (Si,S2,5]3) can be regarded as the time-independent eigenvalues 
of Sa/3(T). Accordingly, Sq,^(t) generates time-independent frame invariants: the linear 
invariant reduces to the trace-free shear condition, i.e., tr (SQ,^(r)) = Si(r) -|- S2(t) -|- 
T,s{t) = Si -|- S2 + S3 = 0; the quadratic invariant yields the Gauss constraint, i.e., 

S°^S„^ = S? + S^ + S| + 2Rl + 2Rl + 2Rl = tl + tl + tl = Q . 

Finally, the conserved cubic invariant is given by 

det Sfj^(T) = S1S2S3 = const . 

This invariant can also be expressed in terms of a Weyl scalar. The magnetic Weyl tensor 
is identically zero for Kasner, but the electric part is non-trivial and thus one obtains one 
non-zero quadratic Weyl scalar, related to det(So/3) according to 

Wi = "^g^^ — = 2- det{T,afs) = const . 

Note that Cabcd.C'^^^'^ is the Kretschmann scalar associated with the metric g while the 
quantity CabcdC"'^'^'^ / is the Kretschmann scalar connected with G. 

Since det Sq,^ = S1S2S3 and (Si, S2, S3) E K^, det S^^ is a constant in the range [—2, 2]. 
Conversely, each value of det S^,^ from the interval (—2, 2) generates a unique ordered 
triple So < S/3 < S^ (a 7^ /? 7^ 7 7^ a), while det Sq,^ = ±2 leads to the LRS points 
S, = ±2, S^ = S^ = ^1 (a / /3 / 7 / a). 

The value of det Sq,^ thus characterizes a Kasner state frame invariantly and uniquely; 
however, it is usually more convenient to use the Kasner parameter u instead, which is 
defined implicitly through 

det(S„^) = 2 - — ■ — = const , u£[l,oo] , (98) 

(1 -|- ti -|- 

where det S^^ is a monotonically increasing function of u. The Kasner parameter u G 
[1,00] parametrizes the one-parameter set of Kasner states {pi,P2,P3) according to 

P°' = T~> — ~, — 2 ' = 1 2 ' Py = "^,^^2 ' (99) 

l + n-|-ti^ I + u + u^ l-|-n-|-n^ 

for sector (a, /?, 7) of the Kasner circle, i.e., Pa < < P"f (Sq < S/3 < S-y); here u £ (1, 00). 
It is easy to check that Pa+Pp+P'y = 1 and = 1- The value u = 1 represents 

the points Qi, Q2, Q3, i.e., the three equivalent representations of the non-flat LRS Kasner 
solution; n = 00 defines the Taub points Ti,T2,T3 and thus the Taub solutionis 

The considerations of this section simplify the analysis of the silent Kasner subset /C, 
see the subsection 'The silent Kasner subset' in Section [6l Since a solution of /C is a 
representation of a Kasner solution in an Iwasawa frame, properties such as n = const 
easily follow. 



^^BKL define the Kasner parameter u according to pi = — tt/(l + u + v?), p2 = {1 + u)/{l + u + u^), 
P3 = u{l + u) / [1 + u + u^) , cf. (|99p : here, the order of the Kasner exponents is fixed as pi < P2 < Pa by 
definition. We find it more natural to use the frame independent definition (|98p . and permute the ordering 
of Pa according to the sector one considers when deahng with frame dependent matters. 
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D Multiple transitions 

Here we give the proofs of some statements made in Section [6] concerning multiple tran- 
sitions; as an alternative to using frame invariants we will use elementary methods from 
the theory of dynamical systems. 

Multiple frame transitions are solutions of (j35p on the silent Kasner subset /C, for which 
at least two of the variables (Ri, R2, R3) do not vanish identically. The equations are 

dr'Si = -2Rl - 2RI dr^2 = -2Rl + 2i?| a^Sg = 2Rl + 2i?^ (100a) 

drRl = (5^2 ~ '^3)Rl + 2ii2^3 (?T-i?2 = (5^1 ~ ^3)R2 drR^ = (Si — S2)i?3 — 2i?ii?2 

(100b) 

together with = 1. The set -R2 = is an invariant subset of (jlOOp (while i?i = 
and i?3 = are not invariant), hence possible classes of multiple frame transitions are 
'^R-iRi (which satisfy i?i ^ 0, R2 = 0, R3 ^ 0) and Tr^r^r^ (where all three functions are 
different from zero). 

First, let R2 = 0, i.e., we consider Tr.^r^ transitions. In this case, Ri = d and i?3 = 
are invariant subspaces; we consider, without loss of generality, the state space defined by 
i?i > and R3 > 0; since T? = 1, the closure of this state space is compact. The functions 
Si and S3 are monotone on the state space, and hence the monotonicity principl implies 
that the o-limit and the w-limit of each orbit must lie on the boundaries Ri = Q (i.e., 
on Br^) or ii3 = (i.e., on Br^) of the state space. Using the known structure of the 
flow on the boundaries (which contain the single transition orbits) we conclude that the 
a-/a;-limit must be a source/sink on K^. The fixed points in four of the sectors of K*~^ 
are saddles; the fixed points in sector (321) are sources, since both Ri and R3 belong to 



the unstable subspaces, see Figure 1(a) , the points in sector (123) are sinks. Thus, Tr^r-^ 



orbits originate in sector (321) and end in sector (123), see Figure 4(a) (Recall that, in 
this paper, 'time' is directed toward the past singularity; the nomenclature 'a-limit set' 
and 'w-limit set' is used in accordance with the chosen time-direction.) 

Second, let R2 ^ 0, i.e., we consider Tr^r^r^ transitions; without loss of generality we 
assume R2 > 0. Since Si and S3 are monotone functions, the monotonicity principle 
yields that the a-/u;-limit of every orbit must lie on the boundary i?2 = 0. Applying the 
previous result on the limit sets of Tr^r^ transitions we see that the a-limit must be an 
equilibrium point in sector (321) and the w-limit a point in sector (123). 

For mixed frame/ curvature transitions T^^r^ the system (j35p reduces to 

dr^i = 2(1 - s2)Si + |iVf drNi = -2{T? + Si)A^i (101a) 

drRl = 2(1 - T?)Ri + (S2 - S3)i?i (101b) 
drT.2 = 2{l-T/)i:2-2Rl-\Nl dr^3 = 2{l - Ti^)Y.3 + 2RI + \nI (101c) 

together with (1 - S^) = (l/12)7Vf ; the remaining variables vanish identically. Without 
loss of generality we consider the case A^i > and i?i > 0. By (jlOip the function (Si + 4) 



^*The monotonicity principle essentially states that if there exists a function that is monotone along 
the orbits of a state space with compact closure, then the a/uj-Vami set of every orbit is contained on the 
boundary of the state space; see [121 [50] and references therein. 
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is monotone on the (relatively compact) state space, and hence the monotonicity principle 
guarantees that the a- and the w-limit of every orbit resides on the boundaries A'^i = or 
Ri = 0. Using the structure of the flow on the boundaries (as given by the single curvature 
transitions of Figure [3]) , we find that the a-limit of each orbit is a fixed point on sector 
(132), and that the w-limit is a point on sector (213) or (231) of K^, see Figure [4(b) [ 

From (jlOip we can derive the validity of the Kasner map (j40p explicitly. The equa- 
tions (llOlap form a two-dimensional decoupled system for Si and A^i by means of the 
Gauss constraint; the underlying reason for this is that Si and A^i are invariant under the 
frame rotations in the (e2, e3)-plane that are induced by Ri. Since the subsystem for Si 
and A^i is identical to the one for the single curvature transition case T/v^, the respective 
solutions are identical; in particular, as in (j39ap . 

Si = 4 + (l + tx2)C, 

where C is defined as in (I39bp and u = U-. This implies that (Si)_|_ is identical to its 
counterpart in the T^Vi case (provided that the initial values (Si)_ are the same). Note, 
however, that S2 and S3 differ from the functions given in (j39ap . since i?i 7^ changes 

the evolution of these quantities. In particular we find that (S2)_|_ ^ ^ = (S3)_|_ ^ and 

(S3)_(_^^^^ = (S2)^^^, which reflects a relative rotation of the axes. However, in terms of 
the Kasner parameter u the two final states are indistinguishable, which is a consequence 
of the frame invariance of u; hence n+ is given in terms of U- by the Kasner map (I40p . 

The equations on the silent Bianchi type VIo/VIIq subsets are given by 

a,S„ = 2(1 - S2)S„ + d^Na = -2(S2 + S,)iV, , 

a = 1...3, a = l,2; here, the quantities '^S^a are polynomials in A^i and A'^2 given by (I35ip . 
The Gauss constraint ( 35j ) reads 

Sf + S^ + Si + ^(Afi -iV2)2 = 6. (102) 

Since ^533 = -(l/3)(iVi - iV2)^ we obtain 

5,(2-S3) = ^(iVi-iV2)'(2-S3), (103) 

i.e., (2 — S3) is a strictly monotonically increasing function unless A'^i = A'^2- 

The Bianchi type VIq subset is given by {N\ > 0) A (A''2 < 0) (or, equivalently, by the 
reversed inequalities). Hence, the Gauss constraint (1102p leads to a compact state space 
given by Sf + S| + S| + (l/2)Ar2 < ^-^^^^ Nx = N wid. N2 = N - ^JVliX - Y?). The 
function (2 — S3) is strictly monotone, since (A^i — A'^2) > everywhere. It thus follows from 
the monotonicity principle that the Q;-/w-limits of orbits must reside on the boundaries of 
the state space, i.e., on YP U B^q^ U B^q^. Using the known structure of the flow on the 
Bianchi type II subsets, see Figure El and the resulting local properties of the fixed points 
on K^, we find that the a-limit of each orbit is the fixed point T3 on K^ and that the 
w-limit is a point on sector (312) or (321). 

The Bianchi type VHq subset is given by (A^'i > 0) A (A''2 > 0) (or, equivalently, by the 
reversed inequalities). In this case, the Gauss constraint does not enforce a compact state 



-71- 



E Behavior of an auxiliary quantity 



space. The derivative of the function (2 — S3) is not positive when A^i = A''2, see ()103p : 
however, 39^(2 - J:s)\n^=n2 = 4(Si - S2)^iVf (2 - S3), and hence (2 - S3) is a strictly 
monotonically increasing function except on the LRS subset A^i = and Si = S2, 
which yields that (Si,S2,S3) = (—1,-1,2) or (Si,S2,S3) = (1,1,-2). The latter is an 
invariant subset with drNa = —ANa, which represents the non-flat LRS Kasner solution in 
a Bianchi type VIIq symmetry foliation. The set (Si,S2,S3) = (—1, —1,2), on the other 
hand, is a one-dimensional set of fixed points, parametrized by Ni = N2 = const > 0, 
which we denote by C^; each fixed point represents the Minkowski spacetime in a Bianchi 
type VIIq symmetry foliation, see [12l pp. 130, 133]. It can be proved, see [511 Theorem 
5], that the a-limit of every orbit in the Bianchi type VIIq state space is one of the fixed 
points on C^. The w-limits are points on sectors (312) or (321) of K^. 

For each individual Bianchi type VIq or VIIq orbit, the quantity |A''i| (and, equivalently, 
1^2!) goes through a maximum value. It is of interest to note that there exists a uniform 
bound, i.e., there exists e > such that max^ |A'^i(t)| > e (and max,- |A''2(r)| > e) uniformly 
for all Bianchi type VIq /VIIq orbits. In order to see this, assume that the assertion is false; 
then there exists for all n S N a Bianchi type VIq (Bianchi type VIIq) orbit Tn, such that 
max,- 1 A''i(r)| < 1/n, i.e., |A^i(r)| < 1/n Vr. In particular, |A^i(rn)| < 1/n, where t„ is 
such that S3(t„) = —1. Since the state space is compact (or, in Bianchi type VIIq, the 
intersection of the state space with the set |A^i| < 1), without loss of generality the limit 
lim„^oo 5]a(r„) exists — otherwise we go over to a subsequence. If the limit is neither of 
the points Qi, Q2, then (Sq,, A^i, A^2)('^n) converges to a point on a TjVj single curvature 
transition orbit. Hence, for sufficiently large n, the corresponding orbit T„ shadows this 
T/V2 transition. However, the flow on the boundary of BniN2- ii3NiN2+) forces In to 
also shadow the TjVi transition preceding/succeeding it. Along these T/Vi transitions we 
have |A^i| ^ (1/?^) for large n, which is a contradiction. The argument is analogous when 
lim„_>oo SQ,(r„) is one of the points Qi, Q2; in this case we employ the saddle structure 
of these fixed points to show that for sufficiently large n, the orbit T„ shadows a T/Vi 
transition along which |A'^i| {^/n), which is again the desired contradiction. 



E Behavior of an auxiliary quantity 

In this appendix we consider Eq. ([57|) . i.e., 

8 

—B = 2(1 - Y?)B + 2Rl - 2Rl ; (104) 

we show that, generically, Bi ~ const as z — > cxd, cf. along a sequence of 

transitions St- (In Section [9] this statement was shown to be true for sequences St that 
do not contain TjVj and TN2R2 transitions, which is equivalent to assuming N2 = 0.) 

For the following it is convenient to use a continuous time variable along Sq-', we introduce 
'sequence time' A by A = (1/2) + (arctan r)/7r, i.e., r = tan [(2A — l)7r/2]; accordingly, the 
time interval AG (i, i + 1) (for i G N) corresponds to r ranging in (—00, 00) during the 
transition. By using sequence time we define the growth function G by 

cIt 

G{\): 9aG = 2(1-S2)— G, G(0) = 1 . 
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By construction, G{X) coincides with G^j) as defined in Section [9] for A = i; this is 
because Gi = nfc=o5'^' "^^ere the growth factors gk are defined as the ratio C+/C- (where 
^-1- = limT-^-too C) that arises from the equation d-rC = 2(1 — for the k^^ curvature 

transition, see (|1T]) and Section [9l 

The derivative of B is closely related to derivatives of S2; in sequence time A we find 

Ai?=/l^(^2-2) + 2(l-S2)^[i?-(S2-2)] wheniV2 = ^^^^^ 



dX 



d_ 



(E2 + 4) + 2(1 - [B - (S2 + 4)] when A^i = 



Integrating (|1U4|) thus amounts to integrating the two equations in pU5|) alternately: as 
long as A''2 = we use the first equation; whenever there is a 7jv2 or a TN2R2 transition the 
second equation is employed. For k £ N, let A^ G N be a sequence of times that represents 
the 'switches' between the two equations in (jl05p . i.e., A''2 = in (Afc,Afc+i) for even k 
and iVi = in (Afc,Afc+i) for odd k. Integration of (11051) and a somewhat cumbersome 
iteration yields 



BiX) = (S2-2)(A) + 6G(A) 



1 



+ G(A)[i?(0)-(S2-2)(0)] 



for A G {Xk,Xk+i) with k even (and a similar result for odd k). By the Leibniz criterion 
the sum Ylj=ii~^y ^i^^j)"^ converges as A; ^ cxd. Moreover, | X^j>fc(— l)-'G'(Aj)^^| < 
G(Aj+i)"^, hence G(A) ^^•^;j.(—l)-'G(Aj)~^ is bounded. Since also (S2 — 2)(A) is bounded, 
we can replace the sum with an infinite series and collect the remaining terms in a term 
labeled 6(A), 



B{X) = G{X) 



00 ^ 
ij(0) + 6^(-iy— — -(S2-2)(0) 
j=i 



+ b{X); 



here b{X) is bounded (and oscillatory). Consequently, generically, B{X) ~ const G{X) 
holds asymptotically, or, in different notation, Bi ~ const 



F Convergence of sums 

Consider a sequence of transitions Sq- on the oscillatory subset O; in our nomenclature, 
St denotes a generic sequence, i.e., a sequence that does not contain any double curvature 
transitions, see Section O therefore, the sequence Sr is associated with a Kasner sequence 
(Mi);gN) where / G N is the index that consecutively numbers the curvature transitions of 
Sr. 

In this appendix we first analyze the series 

^GT'\iAJ:,)\i, (106) 
I 

where Gi is the growth function and | {A'E,2)\i denotes the sum of the absolute changes in S2 
occurring for the frame transitions between the (/ — l)*'^ and the l^^^ curvature transition. 
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Series of the type (|106p arise in connection with the asymptotic suppression of the variable 
A''2, see Section [TO] and the discussion below, and in the context of the asymptotic freezing 
of the Hamiltonian variables rii and Vi, see Section [Hi In the following we prove that 
the series (I106p converges. To establish this result we analyze contributions of the sum in 
large and small curvature phases separately. 

First, we consider (|106p for a large curvature phase (which we assume to be associated 
with a sufficiently large value of r/„): I £ [k, If]. Obviously, |(AS2)|/ is uniformly bounded 
(by 2\/3), so that G;~"^|(AS2)|/ < const G^^ ■ Furthermore, in a large curvature phase, 
5, > Cp > 1, cf. ([62]), so that = g^^G^^ < Gy^G-^ for all I G [luk]. Hence, for 

a large curvature phase, each term in the sum can be bounded by a term of a geometric 
series. 

Second, we consider (jl06p for a small curvature phase; since small curvature phases domi- 
nate over large curvature phases in the probabilistic description of sequences, see Section [HI 
the subsequent considerations are of central signifiance. Recall that the prototype of a 
small curvature phase is the alternating sequence of T/Vi and Tfjg transitions in a neigh- 
borhood of the point T3, which is characteristic for billiard sequences Bq-- The following 
considerations are adapted to small curvature phases of this kind; however, the results hold 
for general small curvature phases with obvious minor modifications. Let I G [/in, lout] be 
a small curvature phase (where we assume > r/„ + 1 to obtain a phase consisting of at 
least two curvature transitions): 



l = hn + l l=kn + l j=kn 



2) I 



Gil E 



^ l + ui + uf \ [ 3(l + 2ni) 



^ 3(l + 2(^'"-fc)) ^ 3L(2^i°-L) 

kn 1 + ttin + (liin)2 'in 1 + ijin + (^in)2 

k=l 



'in ^ 



+ tii" + (ni"") 



where we have used (|60p for the growth factors gj and (|38p for |(AS2)|i. 

From dMI) we see that Gj,^^ < Gj,^^_-^^Gp^; furthermore, Gj2^^2 < Gi.^i'^ - + + 
u™ + (ti™)^), cf. ([6T]) . so that Gj"^ ^2 < ^i-^^p^- This leads to the important conclusion 
that the sum associated with a small curvature phase can be treated on an equal footing 
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with a single term in a large curvature phase, i.e. 



E^r'i(AS2)i/ 



. . . + (AS2)|u-iG^^i + (AS2)|/, 



'out-t" J- 



<™t (• • • + G^^i + Gi;:^ + + • • •) 



< const (. . . + G^_i_i + G^„^_iCp 1 + G^.^_iCp 1 + G^^^iC^ ^ + . . . 



) 



Therefore 



E^r^KASa)!/ < const E '^p^'' < °o 



« k 



which proves the asserted convergence of the series. 

For asymptotic sequences of transitions ASt the above considerations apply in the asymp- 
totic regime when ASq- shadows sequences on O and thus consists of approximate tran- 
sitions (7i)igp|. The convergence result also holds analogously for asymptotic sequences 
ASo- since the possible (approximate) Bianchi type VIq and VIIq orbits are associated 
with a very large increase in G, convergence in the sum is strengthened. 

The established convergence of ()106p has direct applications in our discussion of asymptotic 
freezing in Appendix [Hi The convergence result is also the basis for the proof of the 
assertion made in Section [10] that Eq. ([70]) simplifies to ([7T|) along an asymptotic sequence 
of transitions ASq-', this will been shown in the following: 

Integration of (|70]) . i.e., 



In the asymptotic regime, the integrals can be written as sums over the (approximate) 
transitions of the asymptotic sequence; the second integral contributes only for transitions 
involving We can make the estimate 



where the limits of the integral on the r.h.s. denote the r-time that elapses during the frame 
transitions between curvature transition number / and / -1- 1. In the asymptotic regime, 
the integrals of this type can be approximated by using exact transitions; in particular, 
there exists a uniform upper bound, which suffices to prove convergence in large curvature 
phases in complete analogy to the considerations above. 

The variable Ri is non-zero along Tr^ transitions; along other single transitions it vanishes. 
Note that multiple transitions involving Ri can be ignored in the present context, since 




(108) 



yields 





(109) 
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we focus on small curvature phases, where multiple transitions do not appear. For a Tr^ 
transitions we have — 2i?f = 9rS2, hence —2 /^'^^ Rfdr = (Il2)+ — {^2)-- Therefore, since 
1 2 J^'^^ Ri\ < I I/, the sum in pU9|) can be estimated by a sum of the type pU6|) : we 
thus obtain convergence and 

log {NiRl) = 2 j (2(1 - Y?) -T?- S2) dr + const' 

as r — > 00. Consequently, the term -4i?f G-i(iVii?|) in (fTOH]) can be neglected as r ^ 00; 
this leads to the simplified equation (j7ip and thus establishes the claim. 

To prove asymptotic freezing of ttx in Section [Hi we need to consider series of a slightly 
more general type than (|106p . Consider a sequence of transitions Sr and the associated 
series 

E^rv% (110) 
I 

where a > 0, c > (and for simplicity 2a — c > —1); note in this context that the series 

j:G-iogG.l, E^r^^, 

^ Ul ^ Ui 

see (|134p . can be estimated by a series of the type (|110p (with a = 1 — e, c = 1 — e). In 
the following we prove convergence of the series (jllOp . 

In analogy to (I107p . we obtain during a small curvature phase 



'out+l 'out+1 1 _i_ I 0,2 , „ 'out+1 / \2a-c 

y Gi%ui)-^ = Gi^ y ( (no-<(i+6)G'- y 

Z = «in + 1 «='in + l ' V y «=/in + l ^ ' 



where we have used an estimate of the type {l + u + u^) < (l + e)ii^ for some e > (which 
is small when iju is large). Further, by estimating the sum through the associated integral, 
i.e., 

y {uif"-" = y (u'^'-kf const / u^^-^du, 

we arrive at 

y Gj^inr < const G^;-^_-(^-)i- = const G^;(n^°)^-^ . (Ill) 

Modulo the large curvature phase terms whose sum converges straightforwardly, and which 
hence can be suppressed, the series (jllOp can be written as 

^out, J 

yGY^M-^ = y y GY^{u,r^, 

where [(/in)j) (^out)j] denotes the j^^ small curvature phase. Based on the estimate (jllip 
we therefore obtain 

yG^^im)-^ < const y (G(,„)r" (n';)'~^ , (112) 
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where is the initial value of u for the small curvature phase number j, which begins 
with transition hn,j- 

Regarded as a random variable, = is associated with the probability density w{>i:), 
see ([5T]) . Therefore, = does not possess a finite expectation value, cf. Section [HI 
however, (n™)"^"'^ = >c'^~^ does, since 

((n-)i-^) = r «;(x)x-idx~u;(0)i(r?-i)^~ - r^l'^ . (113) 
Jo c c 

It follows that the expectation value of the sum in pi2|) exists, i.e., 

\j=l / j=l 

Hence, since increases geometrically, i.e., > CpG(;.^)^, the limit n — > oo 

exists. We conclude that the series ()112p and thus (jllOp converges. 



G Hitting intervals stochastically 

Let us introduce a continuous random variable x £ [0, 1] with probability density 

w:[0,l]3 >i ^ w{k) G M , w(Ss) > . 

A simple example is a uniformly distributed random variable, i.e., 'w{>c) = 1. 

Now consider a sequence (5n)neN where (5„ G (0)1) and where 5^ — > (n ^ oo). 
The sequence (5n)neN generates a sequence of intervals [0, 6n] whose length decreases as 
n oo. 

The probability that x G [0, 5„] is given by 

P(xG [0,(5„]) = / w;(x)dx =: 5n ~ u)(0)5„ , (114) 

JO 

where the approximation holds for (5„ <C 1. (Evidently, the result is exact for a uniformly 
distributed random variable.) We denote the event x G [0, 5„] as a 'hit'; the associated 
probability P(x G [0,5„]) = 5„ is the 'hit probability'. 

Consider a series of n = 1, . . . , trials. The probability that no hits occur during N 
trials, i.e., x [0, 5^] Vn = 1, . . . , ('no hits'), is given by 

N 

P^(0) = n(i-^n); 

n=l 

analogously, the probability for 'one hit' is 

N N N ^ A ^ X 

p"a) = E«» n a-*.) = (na-«)Er^ = p"(»)ET^- 

n=l j=lj^n n=l i=l * i=l * 
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Finally, the probability for W hits' is the product 11^=1 ^n- 
As a simple example let 6n = en~^ with e <C 1. We obtain 

N N N 

logP^(O) = log(l -~^n)^Y.i- "^(0)'^" + Oi^n)) = -MO) E - + 0{e') ; 

n=l n=l n=l 

therefore logP^(O) — > — oo as — > oo, and P^(0) — > as A — > oo (independently of e). 
Similarly, 



P^(l) = P»^-^..exp -61. 

i=l ^ (>i \ 




^ 1 

i=l "I \ n=l / L n=l 

therefore P'^(l) — > as A — > oo (independently of e). The analogous result holds for 
the general case: one can show that P^(i) — > as A — > oo for all i £ N. We conclude 
that for an infinite series of trials the probability of getting a finite number i of hits is 
zero, P°°(z) = 0. We expect an infinite number of hits as A ^ oo (although the intervals 
between two hits are expected to become increasingly large). 

As a second simple example let 5„ = en~^ with e ^ 1. We obtain 



N N N 

logP^(O) = J2 log(l - ^~n) « ^ ( - W{0)5n + 0{Sl)) = -ew{0) ^ ^ + 0(6^) ; 



n=l n=l n=l 

in this case, logP^(O) — > — c as A ^ oo (with c > 0), hence P^(0) > as 

A — > oo. Analogously, one can show that (i) — > const > as A ^ oo for all i G N; in 
addition, P^(i) < P^(0)(c*/i!). Accordingly, for an infinite series of trials, the probability 
of getting a finite number i of hits is finite, i.e., P°°(z) > 0, where — > as z ^ oo. 
An important conclusion is that the expectation value of the number of hits is finite. 



(#hits) = ^iP^ 



< oo . 



1=1 



Note that the expectation value exists because falls off sufficiently rapidly as z ^ oo. 

Remark. For the above computations we have made use of the approximation 5n ~ w{0)5n, 
cf. pi4|) . Alternatively, we could have used an estimate of the type 6n < 1 1 it'll oo'^n- 



H Asymptotic constants of the motion 

The analysis of Damour et al. [13] indicates that the off-diagonal Iwasawa frame variables 
Hi, their conjugate momenta Vi, and the momentum vr;^ conjugate to A = log/?, see Sec- 
tionO are asymptotic constants of motion; i.e., these quantities converge to functions that 
only depend on the spatial variables (so-called asymptotic freezing). In this appendix we 
will derive these results from a dynamical systems perspective. 
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Asymptotic freezing of rii and Vi 

We begin by establishing asymptotic freezing of the quantities nj and Vi. In the conformal 
Hubble normalized variables these quantities are given by 

ni = — , n2 = — , ns = — (115a) 

-C/2 -C/3 -C/3 



Pi = -2FV5^3^ + n3p2, V2 = 2H^R2^, V3 = -2H^Ri^ 



(115b) 



where we have used the notation 



El :— Ei^ , E2 :— £'2^ , -£'3 E^' ; 

see Appendix lAl for details. Since ni = — ni, 712 = — n2 + nina, = — ns, cf. ([6]), 
asymptotic freezing of implies that also converge to constants. 

Consider the quantities nj. From the formulas for OqE^, see Appendix it is straight- 
forward to compute 

doni = 2R3^ , don2 = -2i?2^ + , dofis = 2Ri^ . 

Let us now consider a sequence of transitions (or a billiard sequence) that is approximated 
by an attractor sequence At in the asymptotic regime r — > 00. Recall that At is an 
infinite concatenation of TJvu ^-nd T/jg transitions on the billiard attractor C'g_4, see 
Section [TTJ As usual, for notational simplicity we drop the distinction between exact 
transitions and approximate transitions. To lowest order, the evolution equations for nj 
along At are given by 

S^ni = -2i?3^ , drn2 = -2niRi^ , drfis = -2Ri^ . (116) 

^3 ^3 

Let us first focus on the evolution of ni . We rewrite the evolution equation as 

drm = -2Rl , (117) 

and investigate the behavior of the quantity Ei/{E2R3). Along At we obtain 

therefore, Ei/[E2R3) behaves like the inverse of the auxiliary quantity A, see ()52p . Ac- 
cordingly, Ei/{E2R3) decreases like the inverse of the growth function, i.e., 

^ -G,7A(^),, (119) 



^E2R3)i \E2R3J0 

see ()56p . Inserting this result into (jll7p we are able to explicitly compute the change in 
ni for each transition: Ani is zero (in the lowest order approximation) for each Tr^ and 
each TjVi transition and 

(Ani)i = -(;^) . lASsli = -const 0^(1^ |AE2|i , (120) 



-79- 



H Asymptotic constants of the motion 



when % is a T^g frame transition; (AS2)i denotes the change in S2 during this transition; 
here we have used that Ei/{E2Rz) is constant during % and that dr^2 = 2i?|. 

To estabhsh asymptotic freezing of fii we must show that the summation of all (Ani)j is 
finite (where the sum is over all Tfjg frame transitions of At). Let X/jg = \s Tfjg }; 

the series 

^ (Ani), = const ^ G^-(i|AS2|i (121) 

can be rewritten as 

^G,-i|AE2b <oo, (122) 
I 

where / is the running index that consecutively numbers the curvature transitions in Aq-] 
|(AE2)|/ then denotes the sum of the absolute changes in S2 occurring for the frame 
transitions between the {I — 1)^^ and the /^^ curvature transition. The series (jl22p is 
proved to converge in Appendix (F] hence ()12ip converges. Finiteness of the sum (I12ip 
entails that ni converges to a constant value (which depends on the spatial variables), i.e., 
ni exhibits asymptotic freezing. 

The arguments used to establish asymptotic freezing of 77-2,^3 are analogous. Instead 
of Ei/{E2Rz) one can use the quantity E2/{E^Ri), which again exhibits a behavior of 
the type (I119p . In addition, the already established asymptotic constancy of ni can be 
employed in the equation for n2. 

Remark. Asymptotic freezing does not depend on the assumption that the solution con- 
verges to a attractor sequence (although this leads to simplifications in the computations). 
Analogously, one can show asymptotic freezing in the case that the asymptotic sequence 
converges to an arbitrary sequence of transitions on the billiard subset B (including double 
transitions). 

To establish asymptotic freezing of Vi we employ the equation (jllSp and the analogous 
equations for the quantities £'i/(£'3i?2) and E2/ {E'iRi): obviously, for the inverse quan- 
tities E2R3/E1, E3R2/E1, E3R1/E2, the evolution equations are of the form 

drA = 2{l-J:^)A. 

Furthermore we note that the evolution equation for the quantity A = 2H^ is 

5^A = -2(1 -S2)A, (123) 
see Appendix |X1 From these equations it is immediate that 

^rr^ = , 

where we have used that fis converges to a constant in the equation for Vi . It follows that 
Vi are constant in the asymptotic regime so that asymptotic freezing holds trivially. 
Remark. Based on the results of Appendix!^ we are able to compute the following ratios: 

R1R3A _ Vz{nzV2 + Vi) NiRl _ Ni{n3V2+Vi? EiA _ 1 

R2 ~ V2 ' iV2A2 ~ iV2 ' E2R3 ~ n3V2 + Vi ' 

where r.h. sides, and hence also the l.h. sides, are asymptotic constants of the motion. 
Since the quantity A satisfies (j96p . it behaves like the auxiliary quantity A (or like B in 
the generic case); we thus find — a posteriori — the asymptotic behavior of R2 and A''2, in 
consistency with the results of Section [TOl 
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Asymptotic freezing of tta 

The momentum Tr\ conjugate to the variable A is given by 



tta = A log 



(ei)(2-^^)(e2)(2-^^)(e3)(2-^3) 



see Appendix [Aj recall that Cq is defined by = HEa and that A = 2H^. The 
derivative of the variables Cq, reads 



dr^a — (1 + ) 



(124) 



cf. (f95]) in Appendix El Using these formulas and Eqs. ()77ap - ()77cp . and ( |77gP for 9t-Sq,, 
and (jl23p for 3t-A, we find that the derivative of vr;^ is given by 



A-'drTTx = log[(ei)^i {e2)''He3)'''] + 6(1 - S^) + 2Ri + 2i?^ , 



where 



Ai = -12(1 - S^) + 2i?i , A2 = 2{Rl - R 



A, 



-2Rl . 



(125a) 



(125b) 



In order to establish asymptotic freezing of ttx we consider a sequence of transitions that 
is approximated by an attractor sequence At in the asymptotic regime r — > 00 and 
investigate (|125p along this sequence of T/Vi, and T^g transitions. We begin by 

analyzing the behavior of certain functions in terms of the growth function G; we first 
note that 

A = 2H^ = Ao so that A^, = Aq G^^^^ . 

Prom (j77p and (jl24p we obtain the following equations: 



dr [(1 - S2)ef] = 4(1 - S^) [(1 - T?)e\] , 



62 



"i?I^' 


, dr 




^2 







2(1 - lY 



ei 



(126a) 
(126b) 



2(1 - 

hence the relations 

[(1 - S2)et] oc G^ 



describe the behavior of these functions along Ar in terms of growth function G. On the 
basis of these formulas Eq. (jl25p simplifies: For a Tr^ frame transition, (|125p reduces to 



oc , 


/?i^ 


OC G, 




OC G 




^2 




^1. 





a^TTA = A 



2i?ilog^ + 2i?i 

62 



2Ao G-^i?i 



log G + log i?3 + const 



(127) 



The fact that growth function G is constant along the T^g transition makes it possible to 
integrate (|127p . Since 9rS2 = 2i?3 we obtain that 



2 j RldT = AS2 , and 2 j Rj log R3 dr = AS2 (log AS2 - 1) 



(128) 
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where we have used that ii| = [S2 — (S2)-][(S2)+ — S2] for the second integral; as usual, 
AS2 = (S2)+ — (S2)-, where (S2)± denotes the values of 112 at the endpoints of the 
transition. Using (jl28p we find that 



Atta = Ao G"^ AS 



log G + log AS2 + const 



or, in index notation, where we assume that % is the Tr, transition under consideration: 



(A^A)i = Ao GJ^^^ (AS2)i - log G,(,) + log(AE2), + const 



(129) 



here, the constant is independent of i. If 7^ is a Tr^ transition instead of a T^^, the result 
is analogous. 



In the case of a T/Vj curvature transition, (jl25p reduces to 

d^Tix = 6A(l-S2)[-21ogei + l] = 6A0 G~^(l-S2) T-log G+i log(l-S2)+const] . (130) 

The integral of 2(1 — Y?) yields the logarithm of the growth factor g = C+/C- associated 
with the transition T^Vi, i.e.. 



2 y"(l - S2)(ir = j 



-=log^ 



^ogg . 



The integral of 2(1 - S^) log(l - S^) is more involved; using ()39bp we find 

<•(+ I 
'C- 



2y"(l - S2) log(l - S2)dr = (^log ^) log5 + y^"" ^ fog[(C - C-)(C+ - CM . (131) 



If the Kasner parameter u (= n_) of the T/Vi transition is large (i.e., if TJvi belongs to a 
small curvature phase), then p3ip simplifies by using the approximations 



u 



in particular 



_ c 



log[(C-C-)(C+-C)]dC 



12 

u 



log u 



log 6 — 1 



so that 



2 / (1 - S^) log(l - S2)(ir w - [-2 log M + log 12 - 2] 
J u 



for n » 1. Since G is not constant along the T/Vi orbit, Eq. (llSOp cannot be integrated 
explicitly. However, by making the estimate G > G_ , where G_ is the initial value of the 
growth function, we obtain 



IAttaI < 3AoG 



which reduces to 



-1 



log G_ log g + const log g + j - log(l - 



IA^aI <6AoG:^- 



log G_ + log u + const 
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in the case u ^ 1. In index notation, when we assume that the T/Vi transition under 
consideration is the i^^ transition of the sequence, i.e., % = T/Vj, we have 

I AttaIj < 6Ao Gj'^l^ -i- log + log n,(i) + const , (132) 

where the estimate holds in a small curvature phase (with sufficiently large ui). 

To establish asymptotic freezing of tt\ we must show that the summation of all {ATrx)i is 
finite (where the sum is over all transitions of At), i-e., 

^|A7rA|i<oo. (133) 

i 

Finiteness of this sum entails that ttx converges to a constant value (depending on the 
spatial variables), i.e., ttx exhibits asymptotic freezing. 

We merely state the main arguments of the proof and omit the computational details: 
in analogy to the considerations of Appendix [Fj we may focus on the behavior of (|133p 
in small curvature phases (i.e., the alternating phases of 7}jg and T/Vi transitions in the 
neighborhood of T3), since in large curvature phases, the treatment of the sum is simple. 
Since IAS2I ~ when n ^> 1, see (i38]l . the equations (I129p and ()132p are of the same 
type and lead to sums of the type 

E^r^logG.l and E^r^^, (134) 
I ' I ^ 

where / is the running index that consecutively numbers the curvature transitions in At- 
Finiteness of (jl34p is proved in Appendix[Fl hence finiteness of (jl33p ensues and we obtain 
asymptotic freezing of ttx- 



I AVTD singularities and the dynamical systems approach 

In this appendix we discuss the connection between work on AVTD singularities and the 
dynamical systems approach. There exist no rigorous results for inhomogeneous models 
with oscillatory singularities, but there are some for inhomogeneous models with so-called 
AVTD singularities. In this appendix we address how work in this latter field relates to 
the present approach. 

The term velocity dominated singularities was first introduced by Eardley et al. ^2] . Isen- 
berg, Moncrief and co-workers later used the approach of Eardley et al. to obtain rigorous 
results, see [32] and references in [16j . Let us here follow Isenberg and Moncrief |42j . 
They use the standard initial value problem as their starting point, the spatial metric 
and the extrinsic curvature being the dependent variables. They then consider the system 
of differential equations one obtains if one drop all spatial derivatives, except in the Co- 
dazzi constraints, and sets the three-curvature to zero. By doing so they obtain a velocity 
term dominated (VTD) ODE system of evolutionary equations constrained algebraically 
by the Gauss constraint; the solution of this system can subsequently be inserted into the 
Codazzi constraints yielding a VTD-solution. A spacetime is said to be asymptotically 
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velocity term dominated (AVTD) if it is a solution to Einstein's equations such that it, 
w.r.t. some appropriate norm, asymptotically approach a solution to the VXD equations, 
see p87, 88 in [32] for details. If furthermore the asymptotic limit yields a singularity this 
singularity is said to be AVTD. 

If one uses a conformal Hubble-normalized orthonormal frame, then the equations for the 
VTD ODE system are the Bianchi type I evolution equations and the Gauss constraint 
on the silent boundary — and the equations for and the Hubble variable H, with all 
spatial frame derivatives dropped and all variables set to zero except for (and possible 
matter variables), i.e., the lowest order perturbations of the evolution equations for 
and H w.r.t. the silent Bianchi type I subset; this therefore yields the generalized Bianchi 
type I solutions, which in the vacuum case are the generalized Kasner solutions. The 
VTD treatment of the Codazzi constraints corresponds to inserting these solutions into 
the lowest non-zero perturbation of the constraints w.r.t. the silent Bianchi type I subset 
(i.e., one sets all variables to zero except for E'^^S^^, and possible matter variables, in 
the constraints). 

In the vacuum case AVTD singularities occur for special spacetimes, notably ones with 
symmetries, as discussed in Section [131 In cases with one or two commuting spacelike 
Killing vectors, this naturally leads to an Iwasawa frame representation, and hence to 
a fairly simple correspondence between the present Iwasawa based conformally Hubble- 
normalized variables and variables that has been used in previous work in this area (thus 
the AVTD starting point is obtained by plugging in the Kasner circle K^ into the lowest 
order approximation for Ea^ and H, as done in Appendix [B]) . 

Arguably the most impressive result about AVTD singularities was obtained by Andersson 
and Rendall [22j for massless scalar fields and stiff perfect fluids. In this case the Bianchi 
type I equations on the silent boundary yields the generalized Jacobs solutions, and it 
turns out that an open subset of these solutions are stable and hence attract an open 
set of solutions. It should, however, be pointed out that the approach of Andersson and 
Rendall is naturally related to a Fermi propagated frame rather than an Iwasawa frame. 
Stiff perfect fluid spacetimes has also been investigated in terms of the UEWE dynamical 
systems approach by Coley and Lim [53j . 

We note that within the present framework AVTD singularities are more appropriately ge- 
ometrically described as singularities that are associated with (Hubble) asymptotic (Hub- 
ble) conformal spatial flatness. 

To study when singularities are not AVTD Moncrief and coworkers have introduced the 
"Method of Consistent Potentials" (MCP) [431 115] . and references in [iBj. In this case 
one inserts the VTD solutions (the generalized Kasner solutions in the vacuum case) into 
the full equations (or rather, into the Hamiltonian for the full equations); if all terms 
are decaying the spacetime is AVTD, if not it is conjectured to have local Mixmaster 
dynamics. We note that this corresponds to the linear perturbation of the Kasner circle 
K^ in the present context, as can be seen by inserting the generalized Kasner solutions 
into the relations given in Appendix [A] and comparing with the linearized K° result. Here, 
however, this result is seen in the context of the full conformally normalized state space 
picture. 

We note that the Hamiltonian approach using the 'dominant' Hamiltonian 7i suggests a 
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direct generalization of the VTD/AVTD approach, or equivalently, the bilhard subset and 
its lowest order perturbation into the physical state space (with the ODE solution, yielding 
a lowest order perturbation, inserted into the constraints), see Section [T2j An advantage 
of the present approach is that it naturally splits the problem into two parts and hence 
offers a possibility of achieving more modest goals than cracking the whole problem: it 
may be possible to prove some statements that are connected with the silent boundary 
only — it may not be necessary to also prove things about perturbations thereof. 
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